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Experimental Study of Cooling by Injection 
of a Fluid Through a Porous Material’ 


POL DUWEZ? anp H. L. WHEELER, JR. 
California Institute of Technology 


SUMMARY 


A method of coping with the high rates of heat transfer en- 
countered in the field of jet propulsion consists of making the 
parts to be cooled of a porous material and forcing the cooling 
fluid through the pores in a direction opposite to that of the heat 
flow. This method is referred to as ‘‘sweat-cooling.”” The study 
of sweat-cooling heat transfer discussed in this paper has the 
limited objective of defining the variables involved in the problem 
and of establishing the experimental téchniques. A full de- 
scription of these techniques is given. Extensive tables and 
graphs of the data obtained are also included. This experimental 
study of sweat-cooling has shown that either a liquid or a gas may 
be used efficiently as a cooling fluid. A theoretical expression is 
derived relating the surface temperature of the porous material, 
the weight flow rate of the coolant, and the temperature and 
weight flow rate of the main stream of hot gas. The validity of 
this expression is discussed on the basis of the experimental data 
that have been obtained. 


NOMENCLATURE 


B = ty/U 

Cp = specific heat of coolant 

6 = laminar-layer thickness 

D = diameter of sweat-cooled tube 

g = acceleration of gravity 

k = effective conductivity coefficient for coolant 

\ = friction factor for circular pipes 

M = mass flow for turbulent core per unit cross-sectional 
area 

m = mass flow of coolant per unit area of surface 

& = viscosity coefficient in turbulent core 


Presented at the Third National Flight Propulsion Meeting, 
LA.S., Cleveland, March 19, 1948. : 

* This work was done at the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under contracts with the Army 
Ordnance Department, Washington, D.C., and the Air Materiel 
Command, Army Air Force, Wright Field, Ohio. The authors 
wish to thank these agencies for the permission to publish the re- 
sults of this investigation. 

t Associate Professor of Mechanical Engineering and Chief of 
the Materials Section of Jet Propulsion Laboratory. 

} Research Engineer, Jet Propulsion Laboratory. 


viscosity coefficient of coolant 


Me = 

Pre = pclp/k = Prandtl Number of coolant 

Q = weight rate of flow of coolant per unit area 

q = heat transfer across unit area of protective layer 

Re = Reynolds Number of main stream computed from W 
and kinematic viscosity 

pe = density of coolant 

Ps = mass density of turbulent core 

7 = shear stress per unit area within the laminar layer 

7 = Shear stress at wall, determined by bulk characteristics 

T = average temperature of turbulent core 

To = temperature of gas entering coolant manifold 

T, = temperature of sweat-cooled wall 

Tz: = temperature of outside surface of porous specimen 

@ = (T— T) 

U = bulk or average velocity of main stream 

Uy = wall velocity 

V = maximum velocity of main stream measured in center 
of turbulent core 

v = velocity of coolant leaving porous wall 

W = weight rate of flow per unit area in main stream com- 
puted from V, 7, and the density of the gas (con- 
sidered to be air) 

¥ = velocity distortion factor 

§ = temperature distortion factor 


(I) IntTRopucTION 


T= HIGH TEMPERATURES encountered in the opera- 
tion of jet engines have imposed drastic require- 
ments upon the materials used in their construction. 
To meet these requirements, research may be directed 
toward finding alloys and ceramic materials that are in- 
herently refractory or toward materials capable of being 
efficiently cooled. 

Besides the conventional method of cooling by means 
of a liquid in contact with a material exposed to a high 
rate of heat transfer, a less orthodox method consists of 
making the part to be cooled of a porous material and 
forcing a cooling fluid through the pores. In this 


scheme the temperature of the coolant, which moves in 
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Fic. 1. Sweat-cooled test section. 
a direction opposite to that of the heat flow, increases 
gradually while passing through the porous material. 
On emerging from the wall the coolant forms a protec- 
tive layer on the surface exposed to heat transfer. This 
method, referred to as ‘‘sweat-cooling,’’ was proposed 
by the senior author at the Jet Propulsion Laboratory 
in September, 1944, and experimental work on the 
preparation of porous metals was started early in 1945. 
The sweat-cooling method opens new lines of re- 
search in the fields of heat transfer, flow of fluids through 
porous materials, and powder metallurgy. The study 
of the heat-transfer problem discussed in this paper has 
the limited objective of defining the different variables 
involved in the problem and establishing the proper 
experimental techniques. 


(II) DeEscrIPTION OF EXPERIMENTAL EQUIPMENT 


The present experimental investigation was limited 
to the case of a cylindrical duct made of a porous ma- 
terial through the walls of which coolant was injected. 
The test section, containing the porous duct 1 in. in 
diameter and 1'/; in. long, was designed so that it could 
be used with two different burners. In one burner an 
oxyhydrogen flame was used as the source of the stream 
of hot gas, the gas temperature being around 4,200° F. 
and the velocity about 300 ft. per sec. In the second 
burner a gasoline-air flame served as the source of hot 
gas. The gas temperature was only moderately high 
(around 2;000° F. maximum), but very high velocities 
corresponding to Mach Numbers approaching 1.0 


1948 


could be obtained, with resulting Reynolds Numbers 
up to 140,000. While the oxyhydrogen unit simulates 
conditions of temperatures and velocities present in the 
combustion chamber of a rocket motor, the gasoline-air 
burner operates approximately in the range of a gas 
turbine. 


(A) The Sweat-Cooled Cylindrical Duct 


The test section with the sweat-cooled cylindrical duct 
is shown in Fig. 1. The porous cylindrical specimen js 
clamped between the two sections of a stainless-stee] 
holder. Coolant is forced through the walls of the 
specimen from a small annular space between the 
holder and the specimen. The porous specimen jg 
thermally insulated from the holder by transite gaskets, 

For measuring the inside wall temperature of the 
specimen, a chromel-alumel thermocouple is installed 
in the following manner: A '/s-in. diameter hole js 
drilled approximately SO per cent of the way through 
the specimen wall. The remaining 20 per cent is 
drilled with a '/,-in. drill. A conical-shaped bead, ap- 
proximately '/i in. long and with a maximum diameter 
slightly greater than '/1¢ in., is formed on the ends of 
the chromel and alumel wires. This bead is generally 
cast in a carbon mold using silver solder. The bead is 
driven into the small hole in the porous wall until its 
tip is just flush with the inner wall. The wires coming 
from the back of the bead are separated and centered 
in the large hole. ‘The hole is then filled in with a ce- 
ment composed of magnesium oxide and sodium sili- 
cate. After the cement has been baked it is virtually 
impossible to remove the thermocouple installation from 
the wall. This type of installation allows for a maxi- 
mum of strength and thermal contact with the hot sur- 
face of the wall, while affording a minimum of leakage 
around the bead and a minimum of interference to the 
normal pattern of coolant flow through the specimen. 


The outside wall temperature was measured by a 
chromel-alumel thermocouple peened into the outer 
surface of the specimen. All thermocouple leads were 
brought out through the gasket between the two sec- 
tions of the specimen holder. The rate of flow of cool- 
ant was metered by means of a 60-cm. scale Fisher 
Porter Flowrator operated at a pressure of 60 Ibs. per 
sq.in. A series of floats was designed so that a wide 
range of flows could be covered. 


~ 8 Cc 
Fic. 2. Oxyhydrogen burner. 
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(B) The Oxygen-Hydrogen Burner 


The oxygen-hydrogen burner, shown in Fig. 2, con- 
sists essentially of a mixing chamber A, an injector B, 
and a combustion chamber C. The gases are mixed 
by passing them through a porous disc D. This ar- 
rangement produces a considerably hotter flame than 
the commercially available mixers. The injector B is 
the burner tip of a standard No. 15 Victor welding 
torch. The combustion chamber C, 4 in. long and 1 in. 
in diameter, is water-cooled. The water entering 
through the nipple £ cools a portion of the combustion 
chamber and the burner; this water is then discarded 
through nipple F. The water entering the forward 
portion of the chamber and circulating from G to H 
is metered, and its temperature rise is recorded. Hence, 
the heat transfer between the flame and a section of the 
combustion chamber preceding the sweat-cooled speci- 
men may be evaluated. This quantity is taken as an 
index of the conditiows of combustion and is referred 
to-as the potential heat transfer. The test section is 
bolted to the end of the combustion chamber. 

The rate of flow of hydrogen and oxygen through the 
burner was measured by calibrated orifice flowmeter. 
For all of the tests, the total gas flow rate amounted to 
approximately 6 X 10~* Ibs. per sec. at stoichiometric 
mixture ratio. The flame temperature was estimated 
by inserting various refractory materials of known 
melting point into the combustion chamber. By this 
rather crude means, the temperature of the gas stream 
was established as being between 4,000° and 4,500°F. 


(C) The Gasoline-Air Burner 


The experimental burner for the production of high- 
speed and moderately high-temperature gas is shown 
in Fig. 3. It consists of an injection chamber, a com- 
bustion chambey, a ceramic nozzle with a throat section 
1 in. in diameter, and a straightening tube 1 in. in 
diameter and 20 in. long. The axis of the chamber is 
inclined 30° from the horizontal to avoid accumulation 
of fuel. 

The injection chamber contains four tangential air- 
fuel inlets, arranged in pairs. It serves to initiate and 
control combustion and can also be considered as a pilot 
burner. The major portion of the air-fuel mixture is 
supplied, during normal operation, through the front 
inlets of the injector section. Oxidation patterns on 
the interior walls indicate that the gas in the upper end 
of the injection chamber is relatively stagnant, and 
combustion is reinitiated from this area if the flame in 
the combustion chamber is extinguished. 

The combustion chamber, which is flanged and bolted 
to the injection chamber, is primarily intended to 
furnish sufficient volume to complete combustion. 
Two air-fuel inlets are provided at the midpoint of the 
chamber. These inlets were seldom used for the injec- 


tion of fuel, because it was found that a supply of air- 
tuel mixture in this region tended to render the com- 
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Fic. 3. Gasoline-air burner. 


bustion unstable. They have been found of some value, 
however, when used simply as air injectors to cool the 
combustion gases. 

The lower end of the combustion chamber is filled 
with a nest of mild steel tubes 1 in. in diameter and ap- 
proximately 10 in. long. These tubes transform the 
spiral path of the flame into a more or less axial flow. 

At the exhaust end of the combustion chamber the 
gases are expanded through a fused silica nozzle with a 
throat section 1 in. in diameter. Between the nozzle 
and the specimen holder is a heavily lagged steel tube 
20 in. long. The purpose of this tube is to stabilize the 
velocity profile of the hot gas before the gas flows through 
the sweat-cooled section. 

Gasoline is supplied to the burner from a pressurized 
40-gal. tank, controlled by means of needle valves and 
metered through a Fisher Porter Flowrater. Air is 
supplied from two 25-hp. two-stage compressors, 
metered through an orifice and controlled by means of 
gate valves. 

Accurate determination of both the temperature and > 
velocity of the gas stream flowing through the 1-in. 
diameter test section was extremely difficult. In view 
of the preliminary character of the present study, it 
was not felt justified to sacrifice simplicity and rugged- 
ness of the metering apparatus to the accurate deter- 
mination of temperature and velocity. The velocity 
was ineasured only in the center of the stream by means 
of a water-cooled total head tube placed at the outlet 
of the test section. Temperature measurements were 
accomplished with an unshielded chromel-alumel ther- 
mocouple introduced radially into the stream through 
the wall of the straightening tube about 1 in. above the 
test section. A thermocouple recovery factor, obtained 
by comparison with the more delicate Pratt and Whit- 
ney type of probe described in reference 1, was intro- 
duced into the temperature reduction equation, but no 
attempt was made to introduce theoretical corrections 
for radiation losses. The quantities measured by 
thermocouple and total head tube was transformed into 
gas temperature and gas velocity (or Mach Number) 
by means of standard adiabatic equations. 
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(III) REsuLts oF EXPERIMENTS WITH THE OXYGeEy- 


| HypDROGEN BURNER 

| 

All of the experiments with the oxygen-hydrogen 
: { - +L | burner were performed with a gas temperature of ap- 
proximately +,200°F., a velocity of 300 ft. per sec., and 


| 
|: a Reynolds Number of the order of 1,500, as previously 
T mentioned. These conditions resulted in a potential 
| heat transfer of approximately 0.85 B.t.u. per sq.in, 
sec. to the conventionally cooled chamber upstream of 
. | the porous section. The experimental procedure con- 
| sisted of recording the surface temperature of the speci- 
Rete of Nitrogen Flow (10"* ioeqinse) en as a function of the weight rate of flow of coolant. 
R ; The results of a preliminary experiment on the tech- 
vm. 4. in exper: nique of cooling by water injected through a porous ma- 
terial exposed to a flow of gas at high temperature have 
shown (cf. reference 2) that there’is a critical value in 
the amount of coolant, above which the surface tem- 
perature of the porous material remains near or below 
- T the boiling point of water and bélow which the surface 
\ temperature increases rapidly with decreasing flow. 
\ 
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On the other hand, when a gas is used as a coolant, the 
curve of surface temperature vs. amount of gas coolant 
is continuous, since the coolant does not undergo any 
change of phase. 

The results of measurements with copper, nickel, and 
stainless-steel specimens are summarized in the curves 
of Figs. 4 and 5. In these graphs the surface tempera- 
ture is plotted against the flow of coolant. The meas- 
urements of surface temperature vs. rate of coolant flow 
for approximately the same value of potential heat 
P transfer were not always consistent. This is under- 
a 1.5 2.0 2.5 3.0 standable, since the same heat transfer to the cooling 

Rate of Hydrogen Flow (10°4 tb/sq in sec) coil may be obtained for more than one set of values of 
gas temperature and gas velocity, whereas it is prob- 
able that each of these variables taken independently 
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Fic. 5. Surface temperature vs. rate of hydrogen flow in experi- 
ments with the oxyhydrogen burner 


TABLE 1 


Results of Experiments with the Oxygen-Hydrogen Burner 


Flow of 
Coolant Heat Transferred to Coolant 


(10-4 Lbs. —Surface Temperature 7; (°F.)— (B.t.u. per Sq.In. Sec.).—-—__ Exit Velocity of Coolant (Ft. per Sec.) 
per Sq.In. Stainless Stainless Stainless 


Sec.) Copper Nickel steel Copper Nickel steel Copper Nickel steel 
Nitrogen Cooling 
6 900 970 0.122 0.132 . 2.79 2.94 
7 1,150 700 760 0.185 0.108 0.118 3.86 2.78 2.92 
8 910 570 620 0.165 0.097 0.107 3.75 2.82 2.96 
9 740 480 500 0.147 0.089 0.094 3.63 2.89 2.96 
10 630 420 430 0.136 0.084 0.087 3.66 3.01 3.05 
11 550 370 370 0.128 0.079 0.079 3.80 3.12 3.12 
Hydrogen 
1.4 900 740 910 0.389 0.313 0.393 9.07 8.00 9.07 
1.6 750 560 670 0.363 0.260 0.320 9.22 Pe ii 8.61 
1.8 630 440 510 0.336 0.220 0.261 9.34 a ae 8.23 
2.0 550 360 400 0.319 0.190 0.217 9.62 7.81 8.19 
2.2 490 300 . 320 0.306 0.164 0.179 9.95 7.96 8.17 
2.4 440 > 260 260 0.293 0.146 0.146 10.29 8.23 8.23 
2.6 400 230 225 0.282 0.132 0.126 10.65 8.54 8.42 
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has an effect on the conditions of heat transfer in the 
sweat-cooled section. As a consequence, there was 
considerable disagreement in the measurements per- 
formed in various tests with the same sweat-cooled sec- 
tion. 

In spite of this uncertainty in the conditions of heat 
transfer prevailing in the sweat-cooled section, the 
measurements brought out several interesting facts. 
As could be expected, the quantity of gas required for a 
given cooling effect is dependent upon the nature of the 
gas used as a coolant. The curves of Figs. 4 and 5 
show that, for a given wall temperature, the required 
weight rate of flow of nitrogen is approximately five 
times greater than that of hydrogen. 

This ratio of 5 between the hydrogen and nitrogen 
flows is not yet theoretically explained. Some per- 
tinent facts that may serve as a guide in establishing 
a theory explaining the difference between hydrogen and 
nitrogen cooling might be derived from this data. If it 
is assumed that the temperature of the coolant at the 
inside surface of the sweat-cooled tube is equal to that 
of the porous material (cf. reference 3), the quantity of 
heat g absorbed by the coolant before entering the main 
stream of gas is simply 


= — Tr) 


The velocity with which the coolant leaves the porous 
wall will be v = Q/p. This relation is based on the as- 
sumption that the pores are so close together that the 
individual streams of coolant coalesce shortly after 
leaving the porous wall. The values of g and v are 
given in Table 1. From this data, it can be seen that 
for a given surface temperature, the heat transferred 
to the coolant hydrogen is higher than for the coolant 
nitrogen. The reason for this behavior is probably 
related to the thermal conductivity of the two gases and 
the thickness of the boundary layer. 

A comparison of the results obtained with the three 
different metals indicates that heat transfer in sweat- 
cooling depends somewhat on the nature of the porous 
metal. For a given hot-gas stream and for a given 
amount of coolant, the wall temperature of the copper 
specimen was several hundred degrees higher than that 
of the nickel or stainless-steel specimens. This influ- 
ence of the nature of the porous material on the heat 
transfer in sweat-cooling has been evidenced also with 
the tests made with the gasoline-air burner. A quanti- 
tative explanation of this effect is still lacking. 


(IV) ResuLts oF EXPERIMENTS WITH THE GASOLINE- 
AIR BURNER 


The purpose of the measurements made with the 
gasoline-air burner was, to establish a relation between 
the surface temperature of the porous material and the 
Weight rate of coolant flow for different conditions of 
temperature and velocity .in the main stream of hot 
gas. 
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Fic. 6. Surface temperature vs. rate of nitrogen flow for a gas 
temperature of 1,300°F. in experiments with the gasoline-air 
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Fic. 7. Surface temperature vs. rate of nitrogen flow for a 
gas temperature of 1,500°F. in experiments with the gasoline-air 
burner. 


The four significant variables involved in this prob- 
lem are (1) the temperature T of the stream of gas 
flowing through the sweat-cooled tube; (2) the veloc- 
ity V of the same gas; (3) the weight rate of coolant 
flow Q; and (4) the surface temperature 7; of the 
sweat-cooled specimen. The procedure consisted in 
measuring the variation of the surface temperature 7; 
of the sweat-cooled specimen as a function of the rate of 
coolant flow Q for a series of values of velocity V and 
temperature T of the stream of hot gas. The standard 
values of T were 1,100°, 1,300°, 1,500°, 1,700°, and 
1,900°F. At each temperature from two to five differ- 
ent velocities were used. The velocity range was more 
or less dictated by the conditions of operation of the 
combustion chamber. At 1,300°F., for example, the 
velocities were ‘between 650 and 1,400 ft. per sec., 
whereas at 1,900°F. the range was from 1,300 to about 
1,800 ft. per sec. The outside surface temperature 72 


of the sweat-cooled specimen in contact with the in- 
coming coolant gas was also measured in all tests. 

A complete description of the results obtained with 
the four different porous materials is beyond the scope 
of this paper. 


Hence, only the experiments on stain- 
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Fic. 8. Surface temperature vs. rate of nitrogen flow for a 
gas temperature of 1,700°F. in experiments with the gasoline-air 
burner. 


less steel will be presented at length. The fundamental 
graphs expressing the results of these measurements 
show the dependence of surface temperature 7, of the 
specimen on the rate of coolant flow G for given gas 
temperatures and gas velocities (cf. Figs. 6 to 9, in- 
clusive). Curves have also been traced for the varia- 
tion of the temperature 7, as a function of the rate of 
coolant flow. These curves are not included in the 
present paper because they have essentially the same 
shape as those shown in Figs. 6 to 9. 

By comparing all the results obtained with one por- 
ous material, it can be shown that the curves of 7) vs. 
Q fall in a logical sequence, the curves being shifted to- 
ward a higher temperature when either the temperature 
or the velocity of the gas is increased. It is also inter- 
esting to note that the curves tend to be extremely flat 
and approach a horizontal asymptote when the rate 
of coolant flow is increased. 

A comparison of the curves of 7; vs. Q obtained for 
the same conditions of gas temperature and gas veloc- 
ity, but for the four different porous materials, indi- 
cates again that the 7, vs. Q relation depends on the 
nature of the porous sample. 

The results of the measurements, as given in the pre- 
ceding section, were reduced to obtain a more system- 
atic tabulation in which the gas temperature, gas 
velocity, and the rate of coolant flow would have the 
same sets of values for all the tests. The method em- 
ployed was as follows: From the curves of 7; vs. Q, 
presented in Figs. 6 to 9, inclusive, values of tempera- 
ture 7, were read corresponding to the rates of coolant 
flow of 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, and 1.6 X 10~* lbs. 
per sq.in.sec. These temperatures were plotted against 
the gas velocity V. By replotting the data in the form 
of 7; vs. V diagrams, accidental errors could be de- 
tected, since, if one of the original 7; vs. Q curves was 
slightly out of place, none of the points corresponding 
to that curve would fall on the continuous curves of the 
T; vs. V diagrams. Using the 7, vs. V diagrams, new 
values of 7, were then obtained for gas velocities of 
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600, 800, 1,000, 1,200, 1,400, 1,600, and 1,800 ft. per 
sec. A typical set of reduced data thus obtained is pre. 
sented in Tables 2 to 5, inclusive, for porous stainless. 
steel specimens. 
(V) INTERPRETATION OF RESULTS 

Several investigators have recently attempted to 
formulate a theory of heat transfer between a stream of 
gas and a porous wall through which a cooling gas is 
injected. Besides the theoretical treatments of Yuan‘ 
and Rannie,* the junior author has also proposed a 
semiempirical but simple method of correlating the ex- 
perimental results. In the following, the different 
theoretical treatments are briefly summarized and their 
application to the available experimental results is dis- 
cussed. 
In establishing an empirical equation for heat trans- 
fer, the following simplifying assumptions were made 
by H. L. Wheeler: (1) The flow within the tube con- 
sists of a turbulent core and a laminar protective layer; 
(2) the protective layer has reached steady-state condi- 
tions over the length of the tube; (3) values of specific 
heat, thermal conductivity, etc., are average values; 
and (4) radiation heat transfer to the wall is negligible. 
If 6 is the thickness of the laminar boundary layer ina 
conventional pipe, the two following equations can be 
written (cf. p. 46, reference 6): 


(A/8)psU? 


(1) 


TO 


If it is assumed that the shear force on the turbulent 
side of a line dividing the laminar and turbulent regions 
is equal to the viscous force on the other side, then 7 = 
7 and 


5 = 


Denoting the ratio u,,/U by 8 and replacing p,U by the 
mass flow M, the expression for 6 becomes 
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Fic. 9. Surface temperature vs. rate of nitrogen flow for a 
gas temperature of 1,900°F. in experiments with the gasoline-al 
burner. 
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TABLE 2 


‘Reduced Data for a Porous Stainless-Steel Specimen and a Gas Temperature of 1,300°F. 


10 (B.t 
-3 Lbs .t.u. per 
per Sq.In T; Sq.In. (10-8) Q (10-8) 
Sec.) (°F.) (°F.) Sec.) Ti — To T — To V W 
——V = 600 ft. per sec. — R. = 41,000 W = 0.095 Ib. per sq.in. sec. 
1.6 308 180 0.090 4.35 0.81 7.0 16.8 
Oe 334 200 0.088 3.80 0.79 6.3 14.7 
1.2 382 230 0.090 3.04 0.75 5.8 12.6 
1.0 434 270 0.088 2.45 0.71 5.1 10.5 
0.8 506 330 0.085 1.86 0.65 4.4 8.4 
0.6 594 410 0.075 1.37 0.58 3.6 6.3 
0.4 707 530 0.062 0.95 0.48 2.7 4.2 
——V = 800 ft. per see.—— R. = 55,000-—_——_—-~ W = 0.126 Ib. per sq.in. sec. 
1.6 387 210 0.122 2.97 0.75 5.8 12.7 
1.4 418 230 0.117 2.61 0.72 5.3 11.1 
1.2 467 265 0.109 2.27 0.68 4.8 9.5 
1.0 524 310 0.110 1.75 0.64 4.2 7.9 
0.8 593 380 0.102 1.38 0.58 3.6 6.3 
0.6 679 470 ‘ 0.088 1.04 0.51 2.9 4.8 
0.4 788 585 0.070 0.72 0.42 2.1 3.2 
— V = 1,000 ft. per see——. —R, = 69,000-—————~ W = 0.158 lb. per sq.in. sec. 
1.6 451 230 . 0.147 2.29 0.70 5.0 10.1 
1.4 483 260 0.140 2.03 0.67 4.5 8.9 
1.2 535 300 0.135 1.68 0.63 4.1 7.6 
1.0 595 350 0.127 1.37 0.57 3.6 6.3 
0.8 664 430 0.116 1.09 0.52 3.1 5.1 
0.6 748 520 0.088 0.83 0.45 2.5 3.8 
0.4 853 630 0.077 0.59 0.37 1.8 2.5 
——V = 1,200 ft. per sec. R, = 82,000-—————~ W = 0.189 Ib. per sq.in. sec. 
1.6 504 255 0.168 1.88 0.65 4.4 8.5 
1.4 539 285 0.159 1.66 0.62 4.0 e 
1.2 590 330 0.152 1.39 0.58 3.6 6.4 
1.0 650 390 0.141 1.14 0.53 3.2 5.3 
0.8 726 470 0.128 0.89 0.47 2.7 4.2 
0.6 808 560 0.107 0.68 0.40 2.2 3.2 
0.4 907 665 0.082 0.48 0.32 1.6 3.1 
——V = 1,400 ft. per see R. = 96,000. W = 0.221 Ib. per sq.in. sec. 
1.6 550 270 0.186 1.60 0.61 4.0 7.2 
1.4 590 310 0.177 1.39 0.58 3.6 6.3 
1.2 640 365 0.167 1.18 0.54 3.2 5.4 
1.0 700 425 0.154 0.97 0.49 2.8 4.5 
0.8 780 500 0.139 0.74 0.43 2.4 3.6 
0.6 860 585 0.116 0.56 0.36 1.9 2.7 
0.4 955 695 0.087 0.39 0.28 1.4 1.8 


5 = w/(A/88)M (2) 


Consider an element of fluid of unit area and thick- 
ness dy (y being the direction perpendicular to the wall 
of the tube). The viscous shear across dy will be 
udu/dy. If a fluid is injected through the porous wall, 
there will be superimposed on this factor a further stress 
pvy(du/dy) representing the rate at which momentum 
is gained across dy by a mass of coolant p,v flowing out- 
ward from the wall. The total shear stress then in the 
laminar layer will be 


= y(du/dy) + p.v,(du/dy) 


If du/dy may be represented by (u,./5), where y is a 
distortion factor which relates the velocity gradient at 
the wall to the average velocity gradient and is assumed 
to be a function of the ratio m/M, then 7 at a point y = 
will be 


tT = (Uy/5) + 


Taking this new value of shear in place of that given 
by Eq. (1), the thickness of the boundary layer 6 will be 


5 = — (m/p)} (3) 


Eq. (3) is similar to Eq. (2), except for an additional 
term representing the mass flow of coolant and a cor- 
rection for nonlinearity of the velocity profile. 

It is probable that values of \ and 8 do not have the 
same significance in a sweat-cooled tube as in a conven- 
tional pipe. In order to evaluate 6 on the basis of Eq. 
(3), experimental values of 6, ¥, and \ should be ob- 
tained for a sweat-cooled tube in which no heat transfer 
is occurring. 

The case of a sweat-cooled tube in which heat is 
transferred from the protective layer to the metal wall 
will be considered next. For a thermally conservative 
system, the heat transferred through the protective 
layer to the wall is equal to the increase in enthalpy 
of the gas passing through the porous wall, or 


q = &(k/6)(T — Ti) = me,(T, — To) (4) 


where é is a distortion factor relating the temperature 
gradient in the protective layer close to the wall to the 
average temperature gradient in the layer. The factor 
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Eq. (4) is correct, provided the coolant reaches the wall 
temperature. A justification for this assumption has 


been given in reference 3. 
From Eq. (4) it follows that the thickness of the 


boundary layer will be given by 

6 = (&k/me,) [(T — — To)] 
Substituting this relation in Eq. (3) and denoting the 
ratio of constants (u-c,)/k by 


T, — Pre 


Replacing the ratio of mass flows by the ratio of 
weight flows and the temperature ratio by 8, a more con- 
venient form is obtained: 


Reduced Data for a Porous Stainless-Steel Specimen and a Gas Temperature of 1,500°F. 

(10 B.t : 

S. (B.t.u. per 
per Sq.In. Ti Sq.In. v (10-3) Q (10-8) 

Sec.) ( °F.) Sec.) Ti — Ty T-— V W 

—V = 600 ft. per see.———. R. = 35,000-——————— W = 0.085 Ib. per sq.in. sec. 

126 338 185 0.102 4.50 0.82 ia 18.9 

1.4 371 205 0.101 3.88 0.79 6.6 16.5 

1.2 428 235 0.104 3.08 6.75 6.1 14.2 

1.0 487 285 0.101 2.49 O.71 5.4 11.8 

0.8 565 350 0.096 1.86 0.66 4.7 9.5 

0.6 670 450 0.087 1.41 0.58 3.9 yee 

0.4 813 590 0.073 0.94 0.49 2.9 4.7 

—V = 800 ft. per see———. Re = 47,000-—_———— W = 0.113 Ib. per sq.in. sec. 

1.6 414 220 0.132 3.25 0.77 6.0 14.2 

1.4 455 240 0.130 2.79 0.73 5.5 12.4 

1.2 513 275 0.129 2.28 0.70 5.0 10.6 

1.0 574 330 0.122 1.87 0.65 4.4 8.9 

0.8 660 410 0.115 1.45 0.59 3.8 7.5 

0.6 760 520 0.101 1.09 0.52 3.2 5.3 

0.4 886 660 0.080 0.76 0.45 2.3 3.5 
——V = 1,000 ft. per see — —R, = 58,000-——___-. W = 0.142 Ib. per sq.in. sec. 

1.6 476 240 0.157 2.59 0.72 5.1 11.3 

1.4 524 270 0.154 2.20 0.69 +.7 9.9 

1.2 583 310 0.150 1.82 0.65 4.3 8.5 

1.0 643 70 0.140 1.52 0.61 3.8 oe 

0.8 731 455 0.128 1.18 0.54 3.3 5.7 

0.6 828 565 0.110 0.90 0.47 2.6 4.3 

0.4 945 700 0.086 0.64 0.39 1.9 2.8 
——V = 1,200 ft. per sec.— R, = 70,000-—_———— IV = 0.170 Ib. per sq.in. sec. 

1.6 525 255 O: 10s 2.19 0.69 4.5 9.4 

1.4 580 290 0.174 1.84 0.65 4.2 8.2 

1.2 638 335 0.166 1.54 0.61 3.8 : 

1.0 700 400 0.154 1.29 0.56 3.3 6.0 

0.8 784 490 0.140 1.02 0.51 2.8 4.7 

0.6 878 600 0.117 0.78 0.44 2.3 3.5 

0.4 990 730 0.090 0.56 0.36 1.7% 2..4 
——V = 1,400 ft. per sec. — -Re = 81,000-—_—_—_—_-— W = 0.198 lb. per sq.in. sec. 

1.6 568 265 0.194 1.91 0.65 4.0 8.1 

1.4 627 305 0.190 1.60 0.61 3.7 el 

1.2 685 360 0.180 1.35 4.58 3.4 6.1 

1.0 750 430 0.166 1.12 0.53 3.0 o.1 

0.8 827 520 0.148 0.90 O47 2.5 4.0 

0.6 914 620 0.124 0.70 0.41 2.0 3.0 

0.4 1,020 750 0.093 0.51 0.34 * 1s 2.0 
——V = 1,600 ft. per sec——. ‘ R. = 93,000-— —_ W = 0.227 Ib. per sq.in. see. 

1.6 600 275 0.206 1.76 0.63 3.6 rt 

1.4 665 320 0.203 1.43 0.59 3.4 6.2 

1.2 725 385 0.192 1.20 0.55 3.0 S.a 

1.0 790 455 0.176 1.00 0.51 2.7 4.4 

0.8 860 540 0.155 0.82 0.45 2.3 3.5 

0.6 940 640 0.128 0.65 0.39 1.8 2.6 

0.4 1,040 765 0.095 0.48 0.32 1.3 1.8 

f, like y, is assumed to depend on the ratio m/M. 6 = Pr,/\[(E/SBY)(w/Q)] — (6/p)} (5) 


Since the terms on the right-hand side of Eq. (8) 
have been assumed to be functions of Q/W, Eq. (5) may 
be written 


6/Pr, = f(Q/W) (6) 


The validity of Eq. (6) can be discussed on the basis of 
experimental results. 

Values of 6/Pr, and Q/W, given in Tables 2 through 
5, have been plotted on the graph of Fig. 10. As evi- 
denced by this graph, all the experimental results ob- 
tained for a given porous material fall more or less on a 
single curve. The empirical relation (6) seems, there- 
fore, to have some value and is useful in reducing the ex- 
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TABLE 4 
Reduced Data for a Porous Stainless-Steel Specimen and a Gas Temperature of .1,700°F. 
q 
Sec.) (°F.) (°F.) Sec.) T — To T —Te 
-V = 600 ft. per see.—— — —R, = 30,000———————~ W = 0.077 Ib. per sq.in. sec. 
1.6 322 180 0.096 5.69 O85 7.8 20.8 
1.4 358 205 0.096 4.83 0.83 6.5 18.2 
1.2 417 240 0.100 3.81 0.79 6.0 15.6 
1.0 ‘487 290 0.101 2.98 O.75 5.4 13.0 
0.8 570 360 0.097 2.31 0.70 4.7 10.4 
0.6 688 470 0.089 1.66 0.62 3.9 7.8 
0.4 850 0.076 1.10 0.52 3.0 5.2 
——V = 800 ft. per see.——~ ——————-R, = 40,000-————_-. W = 0.103 tb. per sq.in. sec. 
1.6 415 205 0.133 3.84 0.79 6.0 15.5 
1.4 462 235 0.133 3.24 0.76 5.5 13.6 
1.2 525 275 0.132 2.64 0.72 5.0 11.6 
1.0 597 330 0.128 2.13 0.68 4.5 9.7 
0.8 691 410 0.121 1.65 0.62 3.9 7.8 
0.6 810 525 0.107 1.22 0.55 3.3 5.8 
0.4 949 0.186 0.86 0.46 2.4 3.9 
——V = 1,000 ft. per see-——~ —_—_—_——— -R, = 49,000-—-———_—~ W = 0.128 Ib. per sq.in. see. 
ai 493 235 0.164 2.92 0.75 5.2 12.5 
1.4 546 265 0.162 2.48 0.71 4.8 10.9 
‘2 613 310 0.159 2.04 0.67 4.4 9.4 
1.0 684 375 0.150 1.68 0.63 3.9 7.8 
0.8 788 460 0.140 1.29 0.56 3.4 6.2 
0.6 904 580 0. 121 0.97 0.49 2.8 
0.4 1,028 hapes 0.094 0.72 0.42 2.0 3.1 
——V = 1,200 ft. per see—— ————-R, = 59,000——_—_—_—~ W = 0.154 Ib. per sq.in. see. 
1.6 563 260 0.192 2.35 0.70 4.7 10.4 
1.4 618 295 0.187 2.01 0.67 4.3 9.1 
1.2 688 345 0.181 1.66 0.62 3.9 7.8 
1.0 761 420) 0.169 ‘1.38 0.58 3.5 6.5 
0.8 866 510 0.156 1.06 0.52 3.0 5.2 
0.6 979 630 0.132 0.80 0.44 2.5 3.9 
0.4 1,095 a 0.101 0.60 * 0.37 1.8 2.6 
——V = 1,400 ft. per see.—— ———_———-R, = 69,000-——-_—-_—-. W = 0.180 Ib. per sq.in. sec. 
1.6 628 285 0.217 1.96 0.66 4.3 8.9 
1.4 683 325 0.209 1.69 0.63 3.9 7.8 
1.2 753 380 0.200 1.41 0.59 3.3 : 6.7 
1.0 831 460 0.186 1.16 0.54 3.1 5.6 
0.8 934 560 0.169 0.90+ 0.48 2.7 4.4 
0.6 1,042 685 0.141 0.68 0.41 2.2 3.3 
0.4 1,156 Age 0.107 0.51 0.33 1.6 2.2 
——V = 1,600 ft. per sec—_— ——————-R.= 79,000— — W= 0.206 tb. per sq.in. sec 
1.6 689 310 0.242 1.66 0.62 3.9 7.8 
1.4 744 360 0). 230 1.44 0.59 3.6 6.8 
1.2 815 420) 04219 1.20 0.54 3.3 5.8 
1.0 896 505 0.202 0.99 0.49 2.9 4.9 
0.8 996 610 0.182 0.77 0.43 2.5 3.9 
0.6 1,096 740 0.149 0.59 0.37 2.0 2.9 
0.4 1,216 nes 0.43 0.30 1.4 1.9 


0.113 


perimental data to a single curve involving only two 
dimensionless parameters. 

A second theoretical study of the heat transfer in a 
sweat-cooled duct has been reported by Rannie.*’ This 
second study gives the temperature and velocity pro- 
files in the laminar protective layer and leads to a solu- 
tion of the heat transfer through the laminar protective 
layer [cf. Eq. (20), reference 5]. In this solution, how- 
ever, the thickness of the protective layer is a variable 
and is left undetermined. In order to obtain a compari- 


son with the experimental results, it was assumed in 
reference 5 that the boundary-layer thickness was the 
same as that for the flow through an ordinary pipe. 
This hypothesis may be acceptable if the velocity of 
the coolant normal to the wall is small in comparison 


with that of the main stream. For the experiments re- 
ported here, the ratio of these two velocities is smaller 
than 0.01 (cf. Tables 2 through 5). Thus, the final 
results given in reference 5 can be written 


1+ (4 - - 
Voy" 

This expression is Eq. (20) of reference 5, written in 
the nomenclature of the present paper. If it is as- 
sumed that y* = 5.6 and C = 0.046/Re°-*, then the 
quantity (JT — 7,)/(T — To) can be computed as a 


function of Q/W for different values of the Reynolds 
Number. The experimental values of (JT — 17})/- 
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TABLE 5 
Reduced Data for a Porous Staitiless-Steel Specimen and a Gas Temperature of 1,900°F. 
(10- S. .t.u. per : 
per $q.In. Ti T2 Sq.In. (10-3) Q (10-3) 
Sec.) (°F.) (°F.) Sec.) T; — To T —T V W ; 
V = 800 ft. per sec R, = 34,000 W = 0.094 Ib. per sq.in. sec, 
1.6 404 225 0.129 , 4.62 0.82 6.0 17.0 
1.4 470 260 0.135 3.67 0.78 5.6 14.9 
1.2 560 305 0.143 2.79 0.74 5.2 12.8 
1.0 650 370 0.141 2.19 0.69 4.8 10.6 
0.8 748 460 0.132 1.72 0.63 4.1 8.5 
0.6 900 600 0.120 1.22 0.55 3.5 6.4 
0.4 1,048 is 0.096 0.88 0.40 2.6 4.3 
——V = 1,000 ft. per see—— R, = 43,000———————~ W = 0.117 Ib. per sq.in. sec, 
1.6 507 250 0.169 3.26 0.77 5.3 13.7 
1.4 572 290 0.171 2.70 0.73 4.9 12.0 
1.2 657 340 0.172 2.15 0.68 4.6 10.3 
1.0 751 410 0.166 1.71 0.63 4.1 8.5 
0.8 851 515 0.153 1.36 0.58 3.6 6.8 
0.6 994 640 0.134 0.99 0.50 3.0 5.1 
0.4 1,144 0.105 0.71 0.41 2.2 3.4 
——V = 1,200 ft. per sec.-—— R, = 51,000-—————— W = 0.141 Ib. per sq.in. see. 
1°6 588 270 0.202 2.58 0.72 4.8 11.4 
1.4 654 320 0.199 2.17 0.69 4.4 9.9 
1.2 735 380 0.195 1.78 0.64 4.1 8.5 
1.0 832 455 0.186 1.42 0.59 3.7 (fee 
0.8 933 575 0.169 1.13 0.53 3.4 5.7 st 
0.6 1,070 690 0.145 0.84 0.45 2.6 4.3 _ th 
0.4 1,226 ae 0.113 0.60 0.37 1:9 2.8 al 
——V = 1,400 ft. per sec. R, = 60,000. W = 0.164 lb. per sq.in. sec. si 
1.6 655 295 0.228 2.16 0.69 4.4 9.8 
1.4 720 350 0.222 1.84 0.65 4.0 8.5 ec 
1.2 800 415 0.214 1.53 0.60 3.7 7.3 th 
1.0 900 500 0.203 1,22 0.55 3.3 6.1 
0.8 1,000 630 0.187 0.98 0.49 2.9 4.9 
0.6 1,135 730 0.157 0.73 0.42 2.3 3.7 
0.4 1,284 ete 0.119 0.51 0.34 Lo 2.4 
——V = 1,600 ft. per sec. R. = 69,000 -—_———Y~ W = 0.188 lb. per sq.in. sec. 
1.6 714 320 0.252 1.87 0.65 4.0 8.5 
1.4 775 375 ¢ 0.241 1.62 0.62 3.8 7.4 
1.2 858 450 0.231 1.34 0.57 3.4 6.4 
1.0 958 445 0.218 1.07 0.52 3.0 5.3 
0.8 1,060 690 0.194 0.86 0.46 2.6 4.3 
0.6 1,187 780 0.162 0.64 0.39 2.1 3.2 
——V = 1,800 ft. per see.— R, = 77,000————~ W = 0.211 Ib. per sq.in. sec. 
1.6 767 340 0.27. 1.65 0.62 3.7 7.6 
14 825 410 0.359 1.44 0.59 3.4 6.6 
1.2 912 480 0.248 1.19 0.54 3.1 5.7 
1.0 1,012 585 0.231 0.95 0.48 2.8 4.7 
0.8 1,113 740 0.205 0.76 0.438 2.4 3.8 
0.6 1,238 825 0.170 0.57 0.36 1.9 2.8 
‘ (T — Ty), shown in Tables 2 through 5, have been 
- plotted vs. the ratio Q/W in Fig. 11 for porous stainless 
de « steel. “On this figure, two theoretical curves corre- 
ole . sponding to a Reynolds Number of 30 X 10* and 100 
. 
X 10°, respectively, have also been traced. The ex- 
perimental results seem to fall on a single curve having 
34° the shape of the theoretical curves. The different 
j ms symbols corresponding to different gas temperatures 
3 — are too scattered to show the small dependency upon 
wes the Reynolds Number predicted by this theory. More- 
q over, the difference between experimental and theoreti- 
2 +2 wd cal values is, in general, appreciable. 
q _* Yuan‘ made a theoretical investigation of the flow of 
Foal a hot fluid over a sweat-cooled porous plate. He as- 
' c sumed that (1) both the fluid passing over the plate and 
q A the fluid forced through the plate are incompressible; 
7 (2) the viscosities of both fluids are constant; and (3) 
4 12 6 i) 4 no turbulent flow is involved in thé problem. The 
Ratio -2 (10°) momentum and energy equations for the boundary 
Fic. 10. Ratio (T — 71)/(T; — To) vs. ratio Q/W. layer were set up with the velocity of the injection as 


° 
= 
MAS), 
= 
: 


sec, 


ec. 


equation. 
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8 2 16 0 24 
Ratio (10°) 


Fic. 11. Ratio (T — 7;)/(T — To) vs. ratio Q/W. 


sumed to be uniform over the entire plate. In solving 


_ the laminar boundary-layer equation, the momentum 


and energy equations were reduced to an integral form 
similar to the Karman integral relation of the Prandtl 
The solutions were obtained by assuming 
the velocity and temperature profiles in the boundary 


layer as either a polynomial of the fourth degree or an 
exponential function. 

The solutions of the equations gave relationships be- 
tween velocity and temperature layer thicknesses on 
the one hand, and the distance, in the direction of flow, 
from the leading edge of the sweat-cooled plate on the 
other. Calculations were carried out for the tempera- 
ture profiles in the boundary layer at different Prandtl 
Numbers. These calculations were made possible by 
assuming a constant ratio of the thickness of the tem- 
perature layer to that of the velocity layer. This as- 
sumption was justified, since similar profiles were as- 
sumed for velocity and temperature distributions. 
The thicknesses of the temperature and velocity layers 
do not differ appreciably within the range of Reynolds 
Numbers encountered in this investigation. 


Knowing that the quantity of heat removed from 


_ the plate by the coolant is equal to the quantity of heat 


that is transferred to the plate through the protective 
layer, a relation between the wall temperature and the 
rate of coolant flow was established. The final result 
is a linear relation between the temperature ratio @ and 
the square of the ratio of coolant velocity v to average 
main stream velocity U’. 
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Ratio = (1073) 


Fic. 12. Ratio (T — 7,)/(T — To) vs. ratio v/U. 
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A 
50 100 2 
Reynolds Number 
Fic. 13. Nusselt Number ys. Reynolds Number for a 1'/.-in.. 


long duct and an 8-in. long duct. 


Although the theoretical study of Yuan is not exactly 
applicable to the conditions under which the experi- 
ments were carried out, it is interesting to compare the 
theoretical curve of @ vs. v/U with the experimental 
results. As shown in Fig. 12, the measurements ob- 
tained with a stainless-steel specimen are in excellent 
agreement with the theoretical curve computed for a 
Reynolds Number of 30,000. - This agreement may be 
just a coincidence, since, as indicated in reference 4, the 
relation between the temperature ratio 6 and v/ U is de- 
pendent on the Reynolds Number, whereas the experi- 


ments seem to indicate that the Reynolds Number ef- . 


fect is of secondary importance. In addition, as men- 
tioned previously, the results of tests made with porous 
copper, porous nickel, and a porous ceramic material 
did not duplicate the results obtained with porous stain- 
less steel and, if plotted on the graph of Fig. 12, would 
not fall on the theoretical curve. 

Normally, one would not expect the physical proper- 
ties of the porous material to influence the boundary- 
layer conditions. One possible explanation for the ef- 
fect of the material on the conditions of heat transfer 
is the existence of a longitudinal temperature gradient 
along the sweat-cooled wall. It is probable that the 
mixing of the cooling gas with the turbulent core pro- 
duces a progressive decrease in the rate of heat transfer 
along the sweat-cooled section. The result is a temper- 
ature gradient in the longitudinal direction. It there- 
fore seems probable that the thermal conductivity of 
the sweat-cooled material enters into the problem. On 
this basis, it is possible to explain that, while the the- 
ories predict the same coolant flow to obtain a given wall 
temperature, more coolant is actually required to cool 
copper than a ceramic material. 

In trying to correlate the experimental results with 
theoretical studies, the exact conditions in which tests 
were performed should be borne in mind. The rela- 


tively short length of the sweat-cooled section, for ex. 
ample, is an important factor. If a longer tube has 
been used, the surface temperature of the specimen 
would probably have shown a decrease in the down. 
stream direction. Because only one thermocouple was 
used on the inside surface of the 1'/»-in. long specimen, 
this gradient was not observed. Additional experj- 
ments with a longer sweat-cooled section will probably 
throw some light on this aspect of the problem. 


(VI) SUMMARY AND CONCLUSIONS 


The experimental study of heat transfer in a sweat- 
cooled duct has shown that: 


(1) The surface temperature of the porous material 
is a decreasing function of the coolant flow when both 
temperature and velocity of the main gas stream are 
constant. The surface temperature decreases slowly 
when the rate of coolant flow is increased bevond a cer- 
tain value. 

(2) The relation of surface temperature to rate of 


coolant flow is influenced by both the temperature and 


the velocity of the main gas stream. 

(3) A correlation of all the results measured with 
any porous material and cooling gas may be obtained 
by using either one of the ratios (T — 7)/(T, — 7) 
or (T — 7)/(T — Th») and relating it to a function 
Q/W or to v/U. 

(4) Only a fair agreement is obtained between the 
available theoretical studies and the experimental re- 
sults. It seems, however, that most of the difference 
may be due to the fact that the experimental conditions 
do not correspond to those assumed in the theories. 

(5) The experimental results indicate that under 
the conditions present in this investigation, the nature 
of the porous material has a marked effect on the aver- 
age heat transfer to the wall of the sweat-cooled speci- 
men. 


Appendix 


MEASUREMENT OF HEAT TRANSFER TO A CYLINDRICAL 
WATER-COOLED Duct 


The measurements of heat transfer to a convention- 
ally water-cooled duct were made with the purpose of 
obtaining a basis of comparison for the results obtained 
with the sweat-cooled section. Experiments were made 
with two ducts of different length, 11/2 and 8 in. long, 
respectively. Both ducts were 1 in. in diameter. 

The results of the heat transfer tests are summarized 
in a graph of Nusselt Number vs. Reynolds Number, 
presented in Fig. 13. The measurements obtained 
with the S-in. duct show a satisfactory agreement wiih 
the generally accepted relation 


(h/k)D = 0.02R,°-8 (A-1) 
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corresponding to the straight line shown in Fig. 13. 
The measurement obtained with a tube 11/2 in. long 
gave a Nusselt Number higher than predicted by Eq. 
(A-1). The results seem to be on a straight line, which 
would correspond to the relation 


(h/k)D = 0.95R.°:8 (A-2) 

An important result of the heat-transfer measure- 
ments with water-cooled ducts is that the conventional 
relation (A-1) does not apply to a tube 1'/2 in. long and 
1 in. in diameter. Hence, the experimentally deter- 
mined relation between Nusselt Number and Reynolds 
Number (curve 1 of Fig. 13) should be taken as a basis 
of comparison for the sweat-cooling results. 

The fact that the measurements of heat transfer 
made with a water-cooled section 8 in. long agree closely 
with the generally accepted Eq. A-1) is an indirect 
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proof that the measurements of velocity and tempera- 
ture of the main gas stream were reliable. 
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Letters to 


Dear Sir: 

In the May issue of the JouRNAL, F. A. Cleveland presented an 
excellent paper on ‘‘Afterburners for Turbojet Engines’’. 

He analysed afterburning in ducts with constant cross-sec- 
tional area and said that the fuel flow in this case can be selected 
in any arbitrary manner up to a maximum value being limited 
by the airflow available for the combustion [Eq. (16) in his 
paper]. 

This assertion may be correct in certain practical cases but 
does not apply to general conditions, since there exist other 
limits that sometimes can present themselves prior to the above- 
mentioned one. I should like to point out these limits here. 

If we disregard pressure losses due to friction and diffusion 
losses along the combustion chamber, static pressure decrease 
arises according to the energy equation 


—vdp = VdV/g (1) 


where v is the specific volume and V is the gas velocity. This 
results in corresponding total pressure decrease. Considering 
the basic premise A = const. (between stations 6 and 9) we can 
write Eq. (1), after a few transformations and integration, as 
follows: 


= A — yeMol Vo/Ve — A} (2) 
where the term 
= Ve?/gRts (3) 


characterizes the original condition of afterburning (at station 6) 
and where the velocity ratio Vs/Vs is a function of heat input 
Q (fuel flow) and of original condition. The relation Vs/Vs = 
4(Q) can be derived by simple means, but it has no practical 
significance, because a computation using a total heat-entropy 
diagram is more convenient. 

Eq. (2) shows that the pressure ratio p/p. can considerably 
decrease if the term (3) reaches high values. This can happen 
especially at take-off and in flight at sea level. Important de- 


crease in static pressure po with increasing heat input can possibly 
check the latter by the following two limits: (a) the static pres- 
sure py approaches the critical pressure; (b) the static pressure 
bs approaches the ambient static pressure po. 


the Editor 


The second case will occur instead of the first if the pressure 
difference ps — fo is comparatively small. If the fuel flow is 
selected higher than is compatible with one of the above-men- 
tioned limits, a decrease in airflow will result. 

As is to be seen from Eq. (2), the Mach Number J influences 
to a high degree the static pressure decrease, while the velocity 
Ve, which must be low enough to avoid flame blowout, is not the 
only criterion for static pressure decrease (disregarding here fric- 
tion losses, etc), since the value of the static temperature fs accord- 
ing to Eq. (3) must be also considered. In practical computa- 
tions, of course, losses due to friction and diffusion must be added 
to the aforesaid theoretical pressure decrease, and they will cause 
limit (a) or (b) to be reached even sooner. However, it seems 
worth while to compute f»/p. according to Eq. (2) in order to 
estimate real pressure py. 

The Mach Number , can be kept low by suitable diffuser de- 
sign. But first investigation should be made, considering the 
turbojet regulation proceeding, to find out those turbojet operat- 
ing conditions for which the turbine exit Mach Number 14 
reaches its highest value. 

When making up exact computations, there must also be taken 
in account the alteration of the original condition of afterburning 
(station 6) by the additional afterburning fuel flow. 


GUNTHER HELMECKE 
Cité Voisin, Bat.A 
Rue de Faulquiéres 
Decize/Niévre, France 


Dear Sir: 


The article on ‘“‘Buckling of Curved Sheets in Compression and 
Its Relation to the Secant Modulus,” by E. H. Schuette, JouRNAL 
OF THE AERONAUTICAL SCIENCES, Vol. 15, No. 1, January, 1948, 
is extremely interesting and worth while. Certain tests upon flat 
plate elements which were analyzed by the writer showed 
strengthening of such plates through curvature induced during the 
application of load, the certainty and degree of such strengthen- 
ing being related to the stiffening effect of the edge stiffeners. As 
such curvature did not exist initially transversely of the flat 
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plate elements, curved plate behavior had to be determined from 
the variation of curvature with load, failure occurring when the 
load became large ehough to permit buckling but only when the 
induced radius of curvature failed to keep pace with the require- 
ments for curved plate behavior at the stress existing. 

In view of the above observations, correlation of data obtained 
by Mr. Schuette might be improved, I believe, by measuring 
curvature close to the failure load and noting the relation of such 
variations from initial curvature to the R/T ratio of the speci- 
men. At least the relatively wide range of plotted points requires 
narrowing through definite allowance for such effects, ifany. The 
important data is the final, not the initial, curvature. 

To emphasize the point, it seems well to consider somewhat 
briefly the nature of the tests upon flat plate elements which led to 
such observations. 

Twenty-eight test sections were available. Nominal dimen- 
sions of these U-shaped, roof-deck type of beam sections were thé 
same except for two flat plate widths (and several thicknesses), 
the tests being parted into those involving 5.5-in. wide flat plates 
and those 10.5 in. wide. Each group was divisible in turn into 
two groups, one of which had the flanges held in fixed spacing 
laterally at several points spanwise, the other in which there was 
no such lateral restraint and flanges tended to spread laterally 
through a rotation of the flanges about the lines of junction of 
flanges and webs. Suppose R and U define restrained and unre- 
strained groups, respectively, and N and JW’ define the narrower 
and the wider web plates (in compression), respectively, so that 
there might be a trend from most edge restrained to least edge re- 
strained represented by RN, RW, UN, and UW. 

The above tests, be it noted, were part of a program carried out 
by Dr. George Winter at Cornell University for the American 
Iron and Steel Institute. 

It was found that the failure load might be calculated by means 
of the formula JJ, = f,S, to yield a value P, calculated for the 
failure load. J/, is the calculated failure moment, f, is the yield 
point of the material, and S, the section modulus of the reduced 
section. Reduction in aJl cases was based on the formula (6/t) = 
10,320/ a or (w/t), whichever is less, where w is the clear web 
width between flanges and 6 is the equivalent web width, ¢ being 
the thickness. This (b/t) formula corresponds to that for sup- 
ported edge plates. In the case of R sections, this formula gave 
substantially the lower limit of strength (two minor exceptions), 
while the upper limit was substantially given by the formula for 
clamped edge plates using 13,630 for the coefficient, there being 
four cases of exceeding the upper limit. The balance of cases were 
scattered in between. In the case of the U sections, all of the V 
and two of the W behaved as though the (b/t) formula in 10,320 
yielded average values, while the remaining two W sections be- 
haved as though unreduced. 

Now the basis for comparison will be the value of the ratio 
(P,/P,) of the ultimate test value to the calculated value. 

Out of 15 R sections, seven were RN. 
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RN.—With P, based on clamped edges, two had (p,/P,) of 
1.004 and 1.030. The latter value was in the range of (w/t) = 175, 
Below (w/t) = 100 (the stiffer sections), the ratio closely ap. 
proached 1.000. At least one section stiffened somewhat by 
curvature. 

RW.—Of these eight sections, three were strengthened by 
curvature, the ratios being 1.068, 0.992 (where 0.85 was expected), 
and 1.045. These sections were less stiff than RN. 

UN.—Four sections with P, calculated as supported. One 
section, with (p,/p-) = 1.042, was definitely strengthened by 
curvature. 

UW.—Four sections, all of which were effectively strengthened 
by curvature, with (P,/P,) (P, based on supported edges) taking 
the values 1.048, 1.017, 1.173, and 1.374. The last two, or 50 per 
cent of this most flexible group, behaved substantially as unre. 
duced. Basing P, on solid sections, the last two have (P,/P,) 
ratios of 0.966 and 1.176. 

To recapitulate, the number of sections strengthened by in. 
duced curvature were one, three, one, four, the last including two 
that were substantially solid. The numbers are about on the 
basis of relative flexibility, the more flexible, the more numerous, 
More exactly, 1/7, 3/8, 1/4, and 4/4 are the ratios of number 
strengthened to number in the group. 

These analyses were made by the writer because of the un. 
satisfactory character of the Sechler type of plot for the setting 
of allowable effective widths of flat plate elements as used by Dr, 
Winter. The points on the Sechler plot showed an excessive 
scatter, the range of which strongly indicated differences of be- 
havior of flat plates with character of edge stiffening. Further de- 
tails may be found in Proc. A.S.C.E. about February, 1946. The 
conclusion derived from the analyses was that there were marked 
differences in behavior with edge stiffening and that close esti- 
mates of the average failure load (most U group elements) or of 
the minimum ultimate strength (all R group elements) were 
possible with the simplest of effective width formulas derivable 
from buckling theory. (Further tests covered short columns with 
unstiffened flanges which behaved well when P, column values 
were based on an effective flange width of (b/t) 3,672/./f,). 

In view of the above analyses, I feel that it would be well 
worth while reporting radii of curvature of curved sheets close to 
failure and evaluating them with due regard to edge stiffening 
and (R/t) ratios in order to narrow the range of scatter. I feel 
sure that trends may develop which can be eliminated by closer 
analysis. 

The utility of these observations and analyses resides in the 
ability to derive simpler formulas; to calculate more exactly 
minimum strength values; to evaluate edge conditions; and last, 
but probably not least, the concept of a slight imposed initial 
curvature in flat elements to direct and ensure behavior as solid 
elements without any reduced effective width. 


EpwWARD ADAMS RICHARDSON 
Bethlehem Steel Company 


York 21, N.Y. 
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Hyperbolic Column Formulas for 


Magnesium-Alloy Extrusions 


E. H. SCHUETTE* 
The Dow Chemical Company 


SUMMARY 


Results are presented for tests of 89 Z-section columns of M, 
FS-1, J-1, O-1, and O-1HTA alloys. The tests were conducted 
in an effort to arrive at suitable column formulas for use in de- 
sign. In addition to prescribing such formulas, the following 
supplemeutary conclusions were reached: 

(1) The tangent modulus column curve does not give con- 
sistently good fits to test data for magnesium-alloy extrusions. 

(2) The end fixity achieved by a column tested flat-ended is 
probably influenced by the cross-sectional shape of the column. 

(3) Rounds and shapes follow qualitatively different column 
curves. 


INTRODUCTION 


— BEHAVIOR OF SIMPLE COLUMNS within the elastic 
range is well understood. For predicting column 
strength beyond the elastic range, a number of meth- 
ods,.both empirical and theoretical, have been pro- 
posed. For magnesium alloys, the problem of arriving 
at a suitable method has been complicated by a lack of 
adequate test data. A brief survey of column data 
for magnesium alloy extrusions gave indication not 
only that the data were insufficient for obtaining con- 
clusive results but that curves based on tests of rounds 
(comprising a large portion of the total data) would 
not be satisfactory for the design of shapes. As a con- 
sequence of that survey, a new test program, comprising 
89 specimens and covering M, FS-1, J-1, O-1, and O- 
1HTA alloys was undertaken. The results of that 
program are the subject of this report. 


SYMBOLS 
A = area of column cross section, sq.in. 
E = modulus of elasticity, ksi 
E, = tangent modulus of elasticity, ksi 
Fey = compressive yield stress, ksi 
K = empirical constant 
L = column length, in. 
L’ = column effective pin-ended length, in. . 
P = load, kips 
R = minimum radius of gyration of column cross section, in. 
¢ = fixity coefficient in the Euler column formula when 
given in the form, P/A = cr*E/(L/R)? 
n = empirical constant 


ksi = kips (thousands of pounds) per sq.in. 


Received February 12, 1948. 
* Division of Tests, Magnesium Laboratories. 


EXPERIMENTAL WORK 


Specimens and Method of Testing 


To avoid introducing an excessive number of vari- 
ables, a single Z-section shape, Dow die No. 4, was 
selected for all the alloys. The pieces chosen for test 
showed the following average dimensions: 


Flange width (overall): 1.50 in 
Web width (overall): 1.75 in 
Flange thickness: 0.182 in. 
Web thickness: 0.184 in 


Specimen lengths and areas are given in Tables 1-5. 
The ends of the specimens were milled flat, square, and 
parallel. The accuracy obtained by this procedure 
has not been evaluated, but the fact that it is generally 
as accurate as is warranted by the probable straightness 
of the columns is attested to by the consistent results 
and smooth curves obtained. In only two cases, num- 
bers 16 and 32, did it appear, upon placing the eolumns 
in the machine, that the ends were not properly squared 
with the column axis. Since the columns were milled in 
pairs, clamped together, this probably means that num- 
bers 15 and 31 were also below the desired level of ac- 
curacy. All four gave low points on the final plots. 
Moreover, the two specimens specifically noted (16 and 
32) in each case accounted for the low point of the pair. 

The specimens were tested flat-ended on the 30,000- 
Ib. range of a 300,000-Ib. capacity testing machine. 
The ram elevation was maintained essentially constant 
for all tests, and, for the ram setting, loads, and range 
used, the machine is believed to have given an accuracy 
of better than 1 per cent. 

The lower end of each specimen was located with the 
aid of concentric rings on the testing machine head. 
Checking of the vertical alignment was facilitated by 
the fact that a thin oil film on the upper head was 
squeezed out by the pressure of the specimen, leaving a 
clear mark and making it easy to determine whether 
subsequent specimens were bearing in the same place 
as the first. Testing was begun on the shortest speci- 
men, so that a minimum of misalignment was present 
when the mark was first produced. 

The testing machine was run at a rate of 0.3 in. per 
min. until near the maximum load, when the rate was 
reduced to 0.05 in. per min. up to failure. Only the 
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TABLE 1 
Column Test Results for M Alloy Z-Section 
Column F.y(L’/R), Fry(L’/R)?, er 
No. Fry, Ksi L, in. A, Sq.In. P, Kips P/A, Ksi (P/A)/Fy 1.691 XK L (Ksi)!/2 Ksi 
1 9.7 26.40 0.806 6.47 8.03 0.828 44.7 139.2 19,400 
2 9.7 26.40 0.804 7.09 8.81 0.909 44.7 139.2 19,400 
3 9.7 33.92 0.806 7.14 8.85 0.912 57.4 178.8 31,960 
4 9.7 33.93 0.806 6.99 8.67 0.894 57.4 178.8 31,960 
5 9.7 39.46 0.806 6.96 8.64 0.890 66.7 207 .8 43,100 
6 9.7 39.46 0.806 6.95 8.62 0.889 66.7 207.8 43,100 
7 9.7 44.80 0.806 6.89 7.31 0.754 75.7 236.0 55,700 
8 9.7 44.81 0.804 6.23 7.75 0.799 75.7 236.0 55,700 
9 9.7 53.32 0.806 5.07 6.29 0.649 90.1 280.6 78,700 
10 9.7 53.31 0.804 4.83 6.01 0.620 90.1 280.6 78,700 
11 9.7 63.32 0.806 4.20 5.21 0.537 107.0 333.2 110,900 
i2 9.8 63.30 0.807 4.14 5.13 0.524 107.0 335.0 112,200 
13 9.8 73.20 0.807 3.39 4.20 0.429 123.9 388.0 150,400 
14 9.8 73.21 0.807 3.41 4.23 0.432 123.9 388.0 150,400 
15 9.7 86.72 0.806 2.03 2.52 0.260 146.7 457 208,800 
16 9.7 86.7. 0.804 1.81 2.25 0.232 146.7 457 208,800 
7 9.5 106.93 0.809 1.57 1.94 0.204 180.9 557 310,000 
18 9.5 106.91 0.809 1.50 1.85 0.195 180.9 557 310,000 
Colt 
TABLE 2 N 
Column Test Results for FS-1 Alloy Z-Section i 
Column L'/R,  Fey(L'/R)3, 
No. Fry, Ksi L, In. A,Sq.In. P, Kips P/A,Ksi (P/A)/Fy 1691 XL (Ksi)?/2 Ksi 7 
19 13.6 22.19 0.819 11.79 14.41 1.061 37.5 138.3 19,140 7 
20 13.6 22.20 0.816 11.92 14.63 1.078 37.5 138.3 19,140 ‘ 
21 13.7 28.49 0.819 11.80 14.42 1.053 48.2 178.3 31,800 i 
22 13.5 28.49 0.817 11.19 13.70 1.015 48.2 7.4 31,400 8 
23 13.6 33.33 0.819 10.26 12.53 0.921 56.3 207.4 ,000 8 
24 13.6 « 33.32 0.816 11.16 13.69 1.007 56.3 207 .4 43,000 8 
25 13.6 37 .87 0.819 9.19 11.22 0.825 64.0 236.0 55,700 8 
26 13.6 87.85 0.816 9.25 11.33 0.834 64.0 236.0 55,700 8 
27 « 13.7 44.94 0.819 7.38 9.01 0.658 75.9 280.8 78,800 8 
28 13.6 44.93 0.819 7.06 8.62 0.634 75.9 279.8 78,100 8 
29 13.6 53.67 0.819 5.43 6.64 0.489 90.7 334.6 112,000 8 
30 13.6 53.67 0.819 5.33 6.51 0.479 90.7 334.6 112,000 8 
31 13.6 61.97 0.819 4.16 5.08 0.374 104.8 386.2 149,100 8 
32 13.7 61.97 0.819 3.93 4.80 0.350 104.8 388.0 150,400 9 
33 13.6 73.08 0.819 3.55 4.34 0.320 123.6 456 - 207,600 Seenee 
34 13.7 73.08 0.819 3.27 4.00 0.292 123.6 457 208,800 
35 13.5 89.69 0.817 2.35 2.88 0.214 151.8 557 310,000 
36 13.5 89.67 0.817 2.19 2.68 0.199 151.8 557 310,000 el 
or ot 
Meth 
Th 
mine 
TABLE 3 
Column Test Results for J-1 Alloy Z-Section 
Column L'/R, Fey(L'/R)*, 
No. Fy, Ksi L, In. A, Sq.In. P, Kips {P/A)/F,, 1.691 X (Ksi)!/2 Ksi 
37 19.1 18.71 0.788 15.37 19.52 1.022 31.7 138.5 19,190 Since 
38 19.1 18.71 0.787 15.38 19.56 1.024 31.7 138.5 19,190 wav | 
39 19.1 24.11 0.787 14.51 18.46 0.966 40.8 73.1. ~ 31,740 sah 
40 19.1 24.11 0.787 15.33 19.50 1.021 40.8 178.1 31,740 thoug 
41 19.2 28.07 0.792 13.89 17.54 0.914 47.5 208.6 43,500 each 
42 19.1 28.06 0.787 14.20 18.07 0.945 47.5 207 .4 43,000 
43 19.1 31.84 0.788 12.24 15.57 0.815 53.9 235.4 55,400 caleul 
44 19.2 31.84 0.792 11.95 15.10 0.786 53.9 236.0 55,700 
45 19.1 37.99 0.787 9.44 11.99 0.628 64.2 280.4 78,600 
46 19.1 37.99 0.787 9.51 12.09 0.633 64.2 280.4 78,600 
47 19.1 45.27 0.7. 7.438 9.44 0.494 76.5 334.2 111,700 
48 19.1 45.26 0.787 6.99 8.88 0.465 76.5 334.2 111,700 
49 19.1 52.23 0.78 6.09 7.73 0.405 88.4 386.0 149,000 
50 19.1 52.22 0.787 5.87 7.46 0.391 88.4 386.0 149,000 
51 19.1 61.61 0.788 4.52 5.74 0.300 104.2 456 207,600 
52 19.1 61.60 0.787 4.59 5.84 0.306 - 104.2 456 207,600 The 
53 19.2 75.16 0.792 3.21 4.06 0.212 127.2 557 310,000 Calcul 
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TABLE 4 


Column Test Results for O-1 Alloy Z Section 


Column L'/R,  Fey(L'/R)*, 
No. cy» Ksi L, In. A, Sq.In. P, Kips P/A, Ksi (P/A)/Fey 1.691 X L (Ksi)?/2 Ksi 
55 24.1 16.59 0.821 20.65 25.18 1.043 28.06 137.9 19,000 
56 24.1 16.59 0.821 20.67 25.20 1.045 28.06 137.9 19,000 
57 24.1 21.58 0.821 20.52 25.02 1.038 36.5 179.2 32,120 
58 24.1 21.57 0.821 19.60 23.90 0.991 36.5 179.2 32,120 
59 24.1 25.22 0.818 18.05 22.08 0.916 42.7 209.8 44,000 
60 24.1 25.22 0.811 18.32 22.62 0.939 42.7 209.8 44,000 
61 24.1 28.45 0.821 15.62 19.05 0.790 48.1 236.2 55,800 
62 24.1 28.45 0.821 15.90 19.39 0.804 48.1 236.2 55,800 
63 24.1 33.75 0.821 12.37 15.08 0.626 57.1 280.4 78,600 
64 24.1 33.76 0.821 12.64 15.41 0.640 57.1 280.4 78,600 
65 24.1 40.24 0.811 9.26 10.42 0.432 68.0 334.0 111,400 
66 24.1 40.24 0.811 9.23 11.38 0.472 68.0 334.0 111,400 
67 24.1 46.68 0.818 7.32 8.95 0.371 79.0 388.0 150,400 
68 24.1 46.69 0.811 7.27 8.96 0.372 79.0 388.0 150,400 
69 24.1 54.98 0.818 5.66 6.92 0.287 92.9 456 207,600 
7 24.1 54.97 0.811 5.66 6.98 0.295 92.9 456 207,600 
71 24.1 66.94 0.818 3.94 4.82 0.200 113.2 556 309,000 
7 24.2 66.94 0.813 4.07 5.01 0.207 113.2 557 310,000 

TABLE 5 
Column Test Results for O-1HTA Alloy Z-Section 

Column L’ /R, VFiy(L'/R),  Fey(L'/R)*, 
No. Fey, Ksi L, In. A, Sq.In. P,Kips P/A,Ksi (P/A)/Fe 1.691% L  (Ksi)?/2 Ksi 
7 37.5 13.28 0.796 28.66 36.04 0.961 22.46 137.5 18,900 
74 37.5 13.29 0.796 28.42 35.74 0.953 22.46 137.5 18,900 
75 37.5 17.17 0.789 29.05 36.84 0.983 29 .02 177.9 31,620 
7 37.5 17.16 0.789 27.75 35.20 0.939 29.02 177.9 31,620 
77 37.5 20.02 0.7 28.21 35.80 0.955 33.8 7.0 ,800 
7 37.5 20.02 0.789 28.75 36.44 0.972 33.8 207.0 42,800 
79 37.5 22.74 0.789 26.60 33.74 0.900 38.4 235.2 55,300 
80 37.5 22.74 0.789 27.58 34.98 0.932 38.4 235.2 55,300 
81 37.5 27.06 0.790 21.51 27.26 0.727 45.7 280.0 78,300 
82 37.5 27.05 0.796 20.92 26.30 0.701 45.7 280.0 78,300 
83 37.5 32.29 0.796 15.98 20.08 0.536 54.6 334.2 111,700 
84 37.5 32.38 0.796 15.73 19.79 0.528 54.6 334.2 111,700 
85 87.5 37.45 0.790 12.50 15.83 0.423 63.3 387.8 150,300 
86 37.5 37.45 0.790 12.33 15.62 0.417 63.3 387.8 150,300 
87 37.5 44.06 0.796 9.75 12.25 0.327 74.5 456 207,600 
88 37.5 44.06 0.796 9.69 12.17 0.325 74.5 456 207,600 
89 37.5 53.61 0.790 6.45 8.17 0.218 90.6 555 308,000 
90 37.5 53.63 0.790 6.14 7.77 0.207 90.6 555 308,000 


maximum load was recorded; no additional deflection 
or other data were taken. 


Method of Reducing Data 


The area of each length of extrusion used was deter- . 


mined from a single specimen as 


Weight (Ibs.) 
Length (in.) X Density (Ibs. per cu.in.) 


Since, for a given material, the areas determined in this 
way showed less than 2 per cent variation, it was not 
thought necessary to make a separate check of area for 
each individual specimen. The densities used in the 
calculation were as follows: 


M: 0.0635 Ibs. per cu.in. 
FS-1: 0.0645 
J-1: 0.0649 


O-1 or O-1IHTA: 0.0658 


The actual radius of gyration R was not calculated. 
Calculations for a simplified three-piece section, with- 


out radii, and combining the extreme values of recorded 
dimensions, gave the following values: 


Smallest flanges, smallest web: R = 0.3028 in. 
Smallest flanges, largest web: 0.3018 
Largest flanges, smallest web: 0.3060 
Largest flanges, largest web: 0.3074 


These values show a variation of less than 2 per cent. 
Since these maximum and minimum combinations were 
not all actually realized, the true variation was un- 
doubtedly even less than indicated. Consequently, 
the radius of gyration was assumed constant, and its 
actual value not determined. 

The correction for end fixity was obtained from the 
specimens in the elastic range, excluding numbers 15 
and 16, known to be inaccurate. The calculations are 
indicated in Table 6. The average value of (L’/R)/L 
obtained from these specimens was used in reducing all 
the data. With R assumed constant, an intermediate 
calculation for L/R is unnecessary. In order to check 
for any large error, however, the degree of fixity was 
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TABLE 6 
Calculation of (L’/R)/L 


Column (L'/R)/L, 
No. P/A, Ksi L,In. 1/In. 
13 4.20 73.20 123.7 1.690 
14 4.23 73.20 123.2 1.682 
17 1.94 106.93 181.9 1.701 
18 1.85 106.91 186.2 1.742 
35 2.88 89.69 149.2 1.663 
36 2.68 89.67 154.8 1.725 
53 4.06 75.16 125.8 1.673 
71 4.82 66.94 115.4 1.726 
72 5.01 66.94 113.2 1.692 
85 12.50 37.45 63.6 1.699 
86 12.33 37.45 64:0 1.710 
87 9.75 44.06 72.4 1.645 
88 9.69 44.06 72.6 1.650 
89 6.45 53.61 88.6 1.654 
90 6.14 53.63 90.8 1.693 
Average value of (L’/R)/L = 1.691 
ome 
rd LOGARITHMIC PLOT FOR O-l, & FS-i 
+ 
amMaG amcses 7 
20000-—1-+-30 000+} ++-200 000-+—+-300 000 
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calculated using an average of the previously listed ap. 
proximate values of R. The value obtained was ¢ = 
3.77, an entirely reasonable value for this type of test. 


ng. 


Method of Obtaining Column Formulas 


In order to include in the analysis data from other 
sources, and also to aid in obtaining general columy 
formulas, a set of generalized parameters was adopted, 
The failing stress was given in terms of the compressive 
yield stress—that is, as (P/A)/F.,. By dividing both 
sides of the Euler formula by F,,, 


(P/A)/F,, = E/F,,(L'/R)? 


the second parameter, F,,(L’/R)? or V F.,(L'/ R), is ob- 
tained. 

For possible help in arriving at suitable short-column 
formulas, the plots of (P/A)/F.y vs. F.,(L’/R)* were 
made on log-log paper (Figs. 1 to 3), and it was found 
that excellent fits could be obtained on such plots with 
straight lines, indicating formulas of the type 


(P/A)/Fey K[F,,(L’/R)?] 


Included in Figs. 1 to 3 are some additional column 
data on American Magnesium Corporation angle and 
I-beam sections, and an N.A.C.A. I-beam section, ob- 
tained from references 1 and 2. As given in the refer- 
ences, these data were reduced on the assumption of 
full fixity.. Based on the author’s experience, however, 
even though a very nearly uniform contact may be 
obtained between specimen and bearing plate, column 
imperfections and flexibility of the testing machine 
make a fixity coefficient of four unlikely of attainment. 
Consequently, for plotting in this report, the test data 
were adjusted for the following values of fixity coei- 
ficient c, calculated from points in the elastic range: 


Am. Mag. L 3.97 
Am. Mag. I 3.65 
N.A.C.A. I 3.79 


The N.A.C.A. tests were on Dowmetal O-1HTA. 
The American Magnesium Corporation tests were on 
alloy AM-C5S8S, equivalent to O-1, and AM-C58S-T5, 
essentially equivalent to O-1HTA. 

Probably the straight lines shown in Figs. 1 to 3 are 
not the best fits that could be made to each individual 
set of data. They were selected with a view to obtain- 
ing formulas as much alike as possible for the several 
alloys and with exponents readily handled in slide-rule 
calculations. On the whole, little adjustment was 
necessary to achieve these ends; the extreme consist- 
ency of the several test series plotted in Figs. 1 to 31s, 
in fact, surprising. 

Because log-log plots are often misleading as to at 
curacy of curve fitting, the test data and curves ait 
replotted on rectangular coordinates in Figs. 4 to /. 
Here the abscissa has been changed to ~/ F.,(L’/R) 
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material whose compressive yield stresses vary from 
1, ob- | 302 ksi (Am. Mag. L) to 37.5 ksi (Dow Z). In Fig.5, 
reler- the corresponding spread is from 21.2 ksi (Am. Mag. I) 
on of | to 24.1 ksi (Dow Z), and in Fig. 6, from 13.6 ksi (FS-1) |*- 
nies to 19.1 ksi (J-1). The curves in Figs. 5 and 6, more- |" 


7 be over, are identical. 
pn The recommended column formulas, obtained from - 
the plots just presented, are as follows: 
ment. . @ — 
- data M alloy: | | | 
2——+ t 
coel- 
180V F, | 
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has been proposed as a method for predicting short- 
column strength. In a number of instances, this 
method has given good results for aluminum alloys. 
A brief study was therefore made.to determine its pos- 
sible applicability to extruded magnesium alloys. 

In Fig. 8 are shown four tangent-modulus column 
curves for O-1HTA—a curve obtained from the typical 
O-1HTA stress-strain curve given in reference 3 and 
three curves obtained from the stress-strain curves for 
the column material of references 1 and 2. Complete 
stress-strain curves for the material tested in this inves- 
tigation were not available for obtaining tangent- 
modulus values. It is immediately noted that, while 
the data plotted in Fig. 4 are extremely consistent, 
there is considerable scatter in the curves of Fig. 8. 
. Detailed comparison shows also that in only one case 
(Am. Mag. L) is there satisfactory agreement between 
the tangent-modulus column curve and the test data. 

There are two possibilities to consider here. Either 
the differences in Fig. 8 are due to difficulties in obtain- 
ing accurate tangent-modulus values from stress-strain 
curves, or materials possessing different tangent moduli 
are capable of giving column test results of the con- 
sistency evidenced by Fig. 4. It is, however, imma- 
terial which explanation is the more nearly correct 
one; in either case, it appears that the tangent-modu- 
lus method must be considered unsatisfactory because of 
the much better fit obtained by the proposed hyper- 
bolic formulas. 

Consideration was also given to formulas of the 
Rankine type: 


P/A = F.,/[1 + R(L’/R)"] 


If such a formula is to produce a good fit, a plot of 
[F.,/(P/A)].— 1 vs. L’/R on log-log paper should show 
a straight line. A few trials of such a plot were made, 
and they definitely did not show straight lines. Conse- 
quently, because of the relative simplicity of the pro- 
posed hyperbolic formulas, which did give good fits, 
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study of the Rankine-type formula was not pursued 
further. 

With regard to such empirical curves as the straight 
line and the Johnson parabola, a brief study of the 
Figs. 4 to 7 will show that the general shape of the curve 
is such that it could not be well represented by either g 
straight line with a horizontal cutoff or an inverted 
parabola. 


Minimum Property Column Curves 


Because of the general interest in curves for material 
of minimum properties, it is desirable to construct 
graphs based on the proposed hyperbolic formulas, 
incorporating minimum values of compressive yield 
stress. Such minimum values have been prescribed 
for three of the materials tested—namely, O-1HTA, 
O-1, and J-1. The column curves for these materials, 
with minimum properties, are given in Fig. 9 in the 
conventional form of P/A vs. L’/R. 


DISCUSSION 


Effect of Shape on Column Fixity 


In the section on* ‘Method of Obtaining Column 
Formulas,’’ different values were ascribed to the fixity 
coefficient for the two different sections tested by the 
American Magnesium Corporation. This procedure 
might well cause the reader some concern, in view of the 
fact that for the present investigation, a constant value 
was taken for all five materials. However, in this case, 
it is felt that the end justifies the means. The pro- 
cedure gave considerably better fits than would have 
been obtained if an average value of c had been used for 
both sections. It would appear, moreover, that the 
proposed hyperbolic curves could still be justified if 
the American Magnesium data were left out of con- 
sideration. If, then, the use of fairly reasonable as- 
sumptions makes these latter data agree closely with 
the rest, it is well to proceed on the basis of those as- 
sumptions. 

From the above reasoning, one may conclude that 
the fixity obtained in a flat-end column test is very 
likely influenced by the shape of the section tested. 
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FORMULAS FOR MAGNESIUM-ALLOY EXTRUSIONS 


Comparison with Column Data on Rounds 


In the ‘‘Introduction,’’ it was stated that a survey of 
previously obtained data gave indication that column 
data on rounds would be unsatisfactory for the design 
of shapes. The present data now enable a further 
check of that situation. In Fig. 10, therefore, the 
curves obtained herein for shapes are plotted along with 
applicable test data on rounds taken from reference 4. 
It is evident, even though the data on rounds are in- 
sufficient for accurately prescribing the curves, that 
these data show lower slopes than were obtained for 
shapes; the curves leave the horizontal cutoff at 
lower values of F,,(L’/R)? and intersect the Euler 
curve at higher values of (P/A)/ F.y. 


CONCLUSIONS 


(1) The following formulas give excellent fits to the 
test data in the short-column range: 


Malloy: 
180V Fy 
P/A = TR (Max. P/A = 0.90Fy) 
FS-1, J-1, and O-1 alloys: 
2,900V Fey 
= ———— (Max. P/A = F, 
P/A (L’/R)}5 ( ax / cw) 
0-1HTA alloy: 
3,300V F., 
= ———— (Max. P/A = 0.96F 
P/A (Max. P/ ev) 


All stresses are used in kips per sq. in. These formulas 
are recommended for design. 
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= FIGURE 10 
= COMPARISON OF SHAPES AND ROUNDS 


(2) The tangent-modulus column curve does not 
give consistently good fits to test data for magnesium- 
alloy extrusions. 

(3) The end fixity achieved by a column tested 
flat-ended is probably influenced by the cross-sectional 
shape of the column. 

(4) Rounds and shapes follow qualitatively different 
column curves. 
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Letters to 


Dear Sir: 

I wish to congratulate J. F. McBrearty for his interesting 
paper “A Critical Study of Aircraft Landing Gears,” which 
appeared in May, 1948, issue of the JouURNAL OF THE AERO- 
NAUTICAL SCIENCES. 

During World War II, I have had the opportunity of working 
as test engineer with Olaer Comp. in France, maker of hydraulic 
installations for aircraft. I was attached to the drop-test 
department of that company, and the machine in use was similar 
to one shown in the picture in the McBrearty paper, but of smaller 
size. Test results were recorded on moving film. The as- 
sembly of landing gear and shock absorber was dropped onto a 
platform with piezoquartz unit connected to the recording ap- 
paratus. Already at that time some of the aspects mentioned 
and summarized in Mr. McBrearty’s paper were observed. 

Shock absorbers in use were of the hydropneumatic type pat- 
ented by Olaer Comp. The amount of oil air pressure, and valve 
settings were adjusted independently to suit different types of 
aireraft, their load, landing speed, etc. Of course, each specific 


the Editor 


group of shock absorbers had a maximum and minimum value of 
permissible static load (gross weight of the airplane). 

Some of the results are of particular interest to the investi- 
gators of drag failures in the main landing gear initiated by com- 
pression forces in elements designed for tension stresses, as in the 
example mentioned by Mr. McBrearty on page 269 under No. 
(1). Drop tests of the landing gears produce a graph of oscilla- 
tory character mainly dependent on the damping characteristics 
of the shock absorber itself; however, the elastic properties of the 
whole structure and of-the tires have a considerable influence upon 
the final result. For each particular gear assembly with a partic- 
ular shock absorber there were differences in the character of the 
resonance vibrations for different loads and heights of free drop. 

To be sure, these vertical vibrations do not correspond exactly 
to forces and accelerations present during a normal airplane 
landing. It was suggested that ¢ moving platform be constructed 
under the drop-test machine in order to provide the conditions 
similar to those that would occur in real landings. Of course, 
there are many difficulties, but none of them is of insurmountable 
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character at the present time. Electronic devices can be easily 
used for proper recording of the resulting forces, deflections, and 
accelerations in any plane (mainly vertical and horizontal). 

These vertical oscillations in some peculiar landing conditions 
may result in resonance-actuated horizontal oscillations in the 
drag plane, causing the total instantaneous forces to exceed the 
permissible values. From my past experience as a pilot and test- 
flight engineer I know cases where normal landing with great 
forward velocity and without much sinking caused greater 
drag forces on the landing gear than landing with the same air- 
plane and the same load with a high sinking speed. The only 
explanation of this is the combination of the vertical and hori- 
zontal forces together with oscillations from the vertical shock 
absorber which excite critical response at the natural frequency 
of the gear assembly, increasing the oscillatory and spin-up drag 


forces. 


Jerzy ILLASzEwicz, P. ENG. 
Victoriaville, P.Q., Canada 


‘Dear Sir: 

Certain calculations and conclusions presented by Alan Pope 
in his paper ‘“‘On Airfoil Theory and Experiment” (JOURNAL OF 
AERONAUTICAL SCIENCES, vol. 15, No 7, pp. 407-410, July, 1948) 
concerning the theoretical aerodynamic centers of airfoils with 
thickness, are in error. 

The formula for M., the moment about the center of the airfoil’s 
“parent” circle is stated correctly (page 407), but this center is 
not the midpoint of the airfoil, as the author seems to believe. 
In fact, the center is located at a distance (2 + e)c forward of the 
trailing edge, while the mid-chord point is at a distance ¢/2 = 
2c (1 + e?) forward of the trailing edge. 

Thus, the correct formula for calculating A.C., is as follows: 


Mae = Me — Lit — (2+ — (a.c.)t] 


etc., where LZ denotes the lift and ¢ the chord. 
The result, which can easily be confirmed by subtracting 
[(2 + €)c — t/2]/t from Professor Pope’s result, is 


1 1 


or, retaining only O(e?), 


a.c. = (1/4) + (€?/2) 


Thus, von Karman and Burgers? are correct in stating that the 
aerodynamic center lies at the quarter-chord, to the order of 
approximation represented by O(e). Incidentally, it seems ad. 
visable to retain, at most, e? terms, since the expression used by 
Professor Pope for ¢[i.e., £ = 4c(1 + €?)] is an approximation and 
is correct only to this order. 

The numerical values of a.c. according to the corrected formula 
are considerably different from those given by the author; for 
example, 


ac. = 


d/t = 0.09, e = 0.069, a.c. = 0.252 
d/t = 0.15, e = 0.115, a.c. = 0.256 


Obviously, Fig. 2 of the paper requires correction. 

Thus, it is seen that the location of the theoretical aerody- 
namic center (inviscid-fluid approximation) at the quarter. 
chord is not, as Professor Pope has stated, a ‘‘thin-airfoil approxi- 
mation” and is not in the same category as the rough value 2z 
for the lift-curve slope. 
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and no proofs will be sent to the authors. 
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Announcing a New Department . . . Readers’ Gorum 


It has been suggested by several members of the I.A.S. Editorial Committee that the Journal open 
a special department to fill the need for the more immediate communication of essential results of research 
in the greatly expanded and accelerated fields of the aeronautical sciences. Therefore, it has been 
decided to establish a new department, the “Readers’ Forum,”’ which will include letters to the editor 
for the quick publication of (1) brief reports of important findings in aeronautical sciences and (2) dis- 


(1) All letters should be accompanied by a suggested short title, giving some indication to the reader of the subject 


(2) These letters should not exceed 800 words in length. 
(3) Publication will be completed 6 to 8 weeks after receipt of the material. 
(4) The Editorial Committee will not hold itself responsible for the opinions expressed by the correspondents, 


It is hoped that the new department may be started with the January, 1949, issue of the Journal, 
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Deep Ring Analysis 


KUANG-SHENG CHING* 
Bureau of Aircraft Industry, China 


SUMMARY 


A paper with the title of “Calculation of the Stresses in An- 
nular Frames’’ by H. Fahlbusch and W. Wegner! was published 
in the November, 1941, issue of the Journal of the Royal Aero- 
nautical Society. It deals with the depth effect insofar as the neu- 
tral axis of the ring is offset from the skin line, but no considera- 
tion is given to curvature (the distribution of stress due to bend- 
ing being as for a straight beam) or to the deformation due to 
direct load and shear force. 

In the present paper the curved-beam theory and also the 
direct and shear deformations have been taken into considera- 
tion. The result obtained agrees with the above-mentioned 
paper in the region of small ratio of depth to outside radius but 


* differs considerably when the depth-radius ratio is large. 


INTRODUCTION 


tke INSTALLATION of a jet engine in the fuselage. 


usually requires the introduction of a deep ring 
structure that clears the jet pipe and acts as main 
support for the empennage structure. 

Curves are available in the Royal Aeronautical 
Society Data Sheets on Stressed Skin Structures for 
the calculation of circular rings. However, as men- 
tioned in those sheets, these curves are reasonably ac- 
curate provided the ratio of the depth of the ring to the 
outside radius is less than one-tenth. The curved- 
beam theory has to be used to calculate the bending 
stress.2 Therefore it is obvious that we have to con- 


sider the depth effect on the stressing of a thick ring - 


which forms important structural components in a 
single jet engine aircraft. 


SYMBOLS 


= Young’s modulus 

= shear modulus 

outside radius 

= radius to centroidal axis 

= cross-sectional area of the ring 

= cross-section area of the ring to take shear 

= distance between centroidal axis and neutral axis 

= shear flow in supporting skin 

= shear distribution factor across the section in the calcu- 
lation of shear energy? 

K = a(AE/A,G), a constant 

U = strain energy 

M = bending moment, positive producing tension in concave 

side 

V = shear, positive producing positive moment as @ increases 

H = direct force positive producing tension 

M, V,H = applied bending moment, shear, and direct loads 

Mp, Vz, Hp = bending moment, shear, direct load at @ = 0° 


ese 
I 
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Fic. 1. 


Fie. 3. 


The strain energy of a curved beam under combined 
loading for a circular ring is 


y = Me | Het _ Mol, 
Jo \2AEeR,  2A,G. 2AE AER, 


(1) 
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To find Mp, Vr, and Hp, we have 


0U/OMp = 0, dU/@Vp=0, dU/OHp = 0 
i.€., 
JI dM, oM,/R. * R,aM,/ * * ou, (2) 
/10My_ oH, 1 oV, 
In Fig. 1, by statics, fo = : 
Ro So'alRo — R. cos (0 — y)]dy (5) 1+ ~ WE 
(11) 
V, = Vgcos@ — Hysin@ + Ro fog sin (6 —y)dy (6) 14kx- Wx.) 
H, = Hz cos 0 + Vp sin — RoJo’g cos (6 —y)dp (7) 1+ Kwx —wx) 


Two basic cases of loading are considered here. 
Case 1 is a circular ring with an applied moment at 6 = 
0°, and Case 2 is a ring with a vertical shear at @ = 0°. 
Both of these applied load systems are assumed to be 
balanced by the shear from the supporting skin. In the 
second case the shear distribution in the skin is as- 
sumed to be in the form of flexural bending, which is 
nearly true in the case of a deep ring because the ring 
is so much more rigid than the supporting skin. 

In case (1): 


g = —M/2R°r (8) 
In case (2): 


q = (V/Rom) sin (9) 


For case (1), substituting Eqs. (5), (6), (7), and (8) 
into Eqs. (2), (3), and (4), it can be found that 


= 
Ve = —(38M/2rR.)fi 
Ap = 0 
where 
1 1— WX 
In the above equation 
W=R/R, = e/Ro 


Now considering case (2) we substitute Eqs. (5), (6), 
(7), and (9) into Eqs. (2), (3), and (4) and get 


Mr (3 VR,/ 4a) fe 
Ve = 
Hp = (3V/4n)fs 


where 


Substituting these results into Eqs. (5), (6), and (7) 


for case (1) we get: 


= (4 | 


M M thin ring 
1 
-1+ (a - cos0] (13) 
M M thin ring 
+ 2] sin 6 | 
where 


M\ 1 
= —[(x — 6) — 2sing] 
thin ring 2r 


M 
x [1 + 2 cos 6] 
= [2 sin 6] 
and for case (2) 
-2[3 (fe + fa) 


3 + cos 


W (2 
(14) 
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ty 


where 


ring Qr 2 (x 6) sin @+ 1 


1 
(7) x — 6) cos — 5 sin 


(@) = aE cos 6 + (x — @) sin 


If we put X = 0 and K = 0 in Eqs. (13) and (14), 
the results coincide exactly with those in the paper by 
Fahlbusch and Wegner. 

The variables in f,, fe, and fs are X, W, and K, and all 
of them depend on the section properties. For a 
typical I-beam (or box beam) section as shown in Fig. 
3, which the author believes can represent an ordinary 
ring section for heavy load, we find 


{o.1388/ 
1 — 0.9258 
[ozs log, (1 — B) (1 — 0.0758) ” 
0.03 log. 9.9958 

W = 2/(2 — 8) 

~ AG = 18° 
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Here 8 = h/Ro, and we assume that the shear dis- 
tribution factor is equal to one that is nearly true for 
a beam with thin web and that A, is equal to the web 
area. E/G is assumed to be 3 for duralumin or steel 
sheet 

The values of X, fi, fo, and fs have been calculated 
and plotted in Figs. 4 and 5. A curve of X for a solid 
rectangular section is also plotted in Fig. 4 for compari- 
son. 

The values of fi, fo, and f; deduced from Fahlbusch 
and Wegner? are also plotted in Fig. 5; we see that f; 
varies only slightly with a range of 6 for both results 
but f; varies considerably. All the values of fi, fe, and 

fs vary greatly with a large range of B—e.g., from 0.3 
to 0.5. In this graph the x-axis represents the theory of 
thin ring as in ordinary textbooks. 

The values of J/,, Vs, and H, have been calculated 
based on the sectional properties of Fig. 3. The re- 
sults are plotted in Figs. 6, 7, and 8. It is interesting 
to compare these results with those of Fahlbusch and 
Wegner. The variation due to the depth effect in 
M,/VR., V,/V,and H,/V is comparatively small, but as 
in the case of M,/M, V,R,/M, and‘H,R,./M which in 
the paper referred to, do not differ from the thin ring 
theory, however in the present paper, the variation is 
of considerable magnitude. It is obvious that this 
variation is to be attributed to the inclusion of the 
shear and direct load energy and the curved beam effect. 

Therefore the author believes that, in the stressing of 
the deep ring at the end of the fuselage with a jet pipe 
passing through it, the depth effect should be con- 
sidered, otherwise the variation due to depth effect 
might be considerable as demonstrated in Figs. 6, 7, 
and 8. 


Fic. 8. 


REFERENCES 


1 Fahlbusch, H., and Wegner, W., Calculation of the Stresses in 
Annular Frames, J.R.Ae.S., November, 1941. 
? Timoshenko, S., Strength of Material, Vol. II. 


. 


Notice to Contributors 


Manuscripts submitted for publication must be double- or triple-spaced originals. There should be wide margins on 
both sides of the sheets and triple spacing around formulas to allow for the marking of directions for the printer. 


Original drawings (jet black India ink on white paper or tracing cloth) should accompany the manuscript. (Blueprints 
are not acceptable.) Photographs must be on glossy white paper. 


Only the simplest formulas should be typewritten; all others should be carefully written in pen and ink. The difference 
between capital and lower-case letters and Greek symbols should be clearly distinguished. 


See the Inside Back Cover of this issue for further details. 


: 
| 
| \ 
| Me 
} 
| Me f 
| / j | Ir 
| | | | 
co 
| 
| | ° | SP 
| | L | th 
| | | | 4 mi 
| | 
op 
pe 
me 
the 
Ge: rec 
ste 
sel 
od: 
| 
| tu 
| in 
| M 
| te 
| act 
| ta 
| 
= 
| 
tor, 


| 


in 


Application of the Performance Operator 
to Aircraft Automatic Control 


ROBERT C. SEAMANS, Jr.,* BENJAMIN G. BROMBERG, L. E. PAYNE? 
Massachusetts Institute of Technotogy, McDonnell Aircraft Corporation, and Jackson & Moreland 


SUMMARY 


This paper contains a discussion of techniques that provide ef- 
fective means for analyzing the performance of automatically 


controlled aircraft. These techniques are based on the perform- 


ance operator—a concept developed by Dr. C. S. Draper and his - 


associates in the Instrumentation Laboratory of Massachusetts 
Institute of Technology. The performance operator method has 
already been successfully applied in the complete analysis of 
complex Naval fire-control systems. 

Until recently, the fundamental aeronautical design problems 
were largely concerned with fluid mechanics, structures, power 
plants, and static and dynamic stability as applied to aircraft 
controlled by human pilots. The increased emphasis on high- 
speed, automatically controlled flight has forced the designer out 
of these well-explored fields into comparatively uncharted realms 
of air-borne system dynamics. In these regions, the performance 
operator provides a straightforward procedure for understanding 
the necessary physical blending of the sensing and control equip- 
ment with the air frame in a manner to provide satisfactory flight 
maneuvers. 

In this paper, the performance operator is described in general, 
and methods are introduced for determining specific performance 
operators for system components and for complete systems. The 
performance operator of an airplane can be determined experi- 
mentally from information obtained when sinusoidally oscillating 
the control surfaces in flight. It is shown that a simpler method, 
requiring considerably less flight time, involves the application of 
step displacements of the control surfaces. Flight-test data pre- 
sented in the paper, obtained with the U.S. Air Force A-26 
Medium Bomber, indicate good correlation between the two meth- 
ods. Finally, it is shown that the performance operator makes it 
possible to design satisfactory longitudinal control equipment 
for the A-26 airplane. 


INTRODUCTION 


i bes INSTRUMENTATION LABORATORY, of the Division 
of Industrial Cooperation of Massachusetts Insti- 
tute of Technology, working under a contract initiated 
in May, 1946, by the Armament Laboratory of the Air 
Materiel Command, U.S. Air Force, is investigating 
the performance of air-borne tracking systems. Sys- 
tems of this type include radar for directional sensing, 
computers for intelligence, servos or human pilots to 
actuate the controls, and aircraft to carry the useful 
load. Air-borne tracking systems are primarily mili- 
tary in nature, but, as shown in this paper, there are 
many commercial applications. 


Presented at the Instruments Session, Sixteenth Annual Meet- 
ing, ILA.S., New York, January 26-29, 1948. 

* Assistant Professor of Aeronautical Engineering and Direc- 
tor, Tracking Control Project, Instrument Laboratory. 

t Chief, Guidance and Control Section, Guided Missile De- 
partment. 

t Technical Publications Staff. 


In designing a multiloop air-borne system, it is first 
necessary to determine the static and dynamic charac- 
teristics of each component. In order to test the servo 
control units for adaptability to military tracking sys- 
tems, arrangements were made with the Minneapolis- 
Honeywell Regulator Company for the loan of an auto- 
matic pilot. These tests have been carried out both 
in flight and in the laboratory. A fixture was used in 
the laboratory to simulate thé torques occurring about 
the hinge lines of the several control surfaces. 

The Cornell Aeronautical Laboratory pioneered in 
dynamic flight tests and, hence, was logically equipped 
to assist in the flight program. During tests conducted 
by this group, the response of the aircraft was deter- 
mined separately for applied displacements of the ele- 
vator, rudder, and ailerons. In each case, sinusoidal 
displacements and those of a step character were em- 
ployed in order to ensure both steady-state and transient 
data. The Instrumentation Laboratory developed 
special instruments for the test, including a three-axis 
rate indicator, a three-axis accelerometer, and sensing 
elements for angle of attack and angle of yaw. Both 
longitudinal and lateral flight tests have been satisfac- 
torily completed, and the aircraft has been transferred 
recently to the Instrumentation Laboratory for instal- 
lation and test of necessary tracking equipment. 

The method selected to present the longitudinal and 
lateral flight test results is important, since it must ful- 
fill two necessary functions. The method must give 
the engineer and the designer a clear mental picture of 
the physical phenomena associated with the airplane’s 
static and dynamic performances. It is equally im- 
portant that the method yield information that can be 
used readily and economically in establishing the sta- 
bility and the response time of the overall system. For 
example, by techniques developed at Cornell Aeronau- 
tical Laboratory, it is possible to evaluate the important 
stability derivatives from flight tests. Until recently, 
the dynamic stability derivatives had never been ex- 
perimentally determined for the complete aircraft; 
therefore, this result provides the aircraft designer with 
a new tool of inestimable significance. When these 
stability derivatives are inserted in the airplane equa- 
tions of motion, the airplane output for any control 
displacement can be determined. However, the control 
displacements in a multiloop system are functions of 
quantities involving all the components in an auto- 
matically controlled aircraft. In order to obtain the 
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Fic. 1. Functional diagram of the operating component, d. 


response of the complete system, it may be necessary 
to solve as many as 20 simultaneous differential equa- 
tions. Without the aid of an adequate simulator, this 
“classical approach”’ entails a prohibitive outlay of ef- 
fort. 

A technique has been developed by the Instrumenta- 
tion Laboratory that permits convenient determination 
of the response of compound systems. The technique 
is based on a mathematical concept, called the per- 
formance operator. This‘operator expresses the rela- 
tionship between the output of a system and the actuat- 
ing input, or forcing function, responsible for the out- 
put. In use, the performance operator employs con- 
ventional mathematical procedures in such a manner 
that complex functions involving derivatives, vectors, 
and nonlinearities can be segregated and divided into 
understandable effects. 


(1) THE PERFORMANCE OPERATOR 


The performance operator may be broadly defined as 
a concept designed to express in readily identifiable 
form the performance characteristics of an operating 
system or subsystem. ‘Performance characteristics” 
are understood to be those characteristics of a system 
which alter its output with respect to its input or in- 
puts. 

A simple operating component, designated d, is 
shown in the functional diagram of Fig. 1. Impressed 
on the system is a forcing function, or actuating input 
quantity, g... The action of the operating component 
on the actuating input gives rise to an output quantity, 
Gout. Depending upon the nature of the component, 
the output may differ from the actuating input in di- 
mension, magnitudé, phase, and integral or differential 
order. In this sense, the component can be said to 
operate on the actuating input to produce the output. 


! 
; OPERATING - OPERATING H 
| COMPONENT COMPONENT | 
Gin = Aina, d, Moutya, = Aina, dp F out 
ours, = outa, = 
! 


q out = out) Fin 


Fic. 2. Functional diagram of the system, d. 


_input quantity, 
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In terms of the performance operator this may be ex- 
pressed : 


The brackets ¢ } designate the term so enclosed as an 
operator; the form ¢P} indicates that the term is a per- 
formance operator relating the magnitude of the output 
to the magnitude of the actuating input; the subscript 
d indicates the designated component; and the syb- 
scripts gj, and ou: indicate by their order and form the 
particular performance operator intended. 

In the operating component of Fig. 1, an additional 

factor has been introduced in the form of a modifying 
A modifying input is dis- 
tinguished from an actuating input by noting that it 
can have no effect upon the output unless an actuating 
input is present; whereas, an actuating input can af- 
fect the output in the absence of a modifying input. 
The effect of a modifying input is to change the per- 
formance operator. 

The location of the center of gravity of an airplane 
is an input that modifies the longitudinal response of 
the airplane to applied elevator displacements. In 
this example, the elevator displacement is the actuating 
input quantity. There are several output quantities 
that may be of interest—namely, rate of pitch, normal 
acceleration, and angle of attack. For this reason, at 
least three performance operators are needed in speci- 
fying the longitudinal characteristics of a particular 
aircraft. Each one of these performance operators is 
modified by the center-of-gravity position. 

In practice, several operating components are usually 
found combined into a more complex system. The 
scheme of concepts, definitions, and notation chosen to 
describe the behavior of the system must be effective in 
permitting immediate, unambiguous specification of 
the performance characteristics of any selected operat- 
ing component of the system. This makes it necessary 
for the notation associated with each performance oper- 
ator to include complete information on the identity of 
the operating component involved, its output, and its 
actuating input. In each case, in addition to the per- 
formance operators, it is expected that a functional dia- 
gram (including a block for each operating component) 
is available for reference. From this diagram, it is 
possible to identify the modifying inputs that affect 
the performance operator associated with a given oper- 
ating component. 

The performance operator of a system comprising 
any number of operating components—whether con- 
nected in series, parallel, or with feedback returns— 
can be related to the component operators by carrying 
out the straightforward algebraic processes indicated 
on the pertinent functional diagram. Fig. 2 shows two 
operating components, d; and d2, connected in series to 
form an open-end system, d. In this figure the per- 
formance operators for the individual components are 
specified, and the simple relationships between com- 
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ponent inputs and outputs are noted. However, ref- 
erence to the functional diagram indicates that the in- 
put to operating component d; is also the input to sys- 


tem d, that the output from operating component d; is - 


the input to operating component d2, and that the out- 
put from operating component d2 is the output from 
system d. It then becomes evident that 


$Phatgindoud = (2) 


It has already been stated that the performance char- 
acteristics of a system are understood to be those that 
alter its output with respect to its input. When the 
input and the output are scalar quantities, the operator 
is also scalar and serves only to account for changes in 
magnitude. However, when the input and the output 
are vector quantities having both magnitude and direc- 


tion, the operator not only must relate the magnitude of: 


the output to the magnitude of the input but must in- 
troduce a space orientation of the output vector with 
respect to the input vector. When this is the case, the 
vector nature of the performance operator is designated 
by the addition of a bar over the symbol—i.e., 


This equation is an extension of Eq. (1) in the sense that 
it permits a study of systems that experience derivative 
and integral effects and, in addition, are involved with 
space-dimensional geometry. The open-parentheses 
appearing as a subscript to the performance operator 
symbol in this equation signifies the designated compo- 
nent or system. 

For the mathematical treatment of certain vectors, 
it is expedient to refer to the complex plane. This 
plane is familiar to most engineers as that defined by an 
axis of reals and an axis of imaginaries. All real num- 
bers fall on the axis of reals, and numbers that are 
the product of real numbers and +/—1 (which repre- 
sents a unit vector along the axis of imaginaries and is 


‘ usually represented by the symbol j) fall on the axis of 


imaginaries. The complex plane derives its title from 
this dual nature of numbers therein. 

In general, a quantity in the complex plane is com- 
posed of both a real and an imaginary number, and its 
magnitude is a function of the size of each. Because its 
graphical representation may fall anywhere between 
the axis of reals and the axis of imaginaries, the quantity 
also has direction. From this point of view, a quantity 
occurring in the complex plane can be considered as a 
vector and is so classed in this paper. Since this type 


of vector is always contained within a plane, certain 
common vector operations—for example, the vector 
cross product—which involve three-dimensional space 
are never required. However, despite limitations of 
this kind, the vector remains an extremely useful con- 
cept to the aeronautical engineer in dealing with com- 
plex numbers. 
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Fic. 3. Diagram showing the derivation of the vector perform- 
ance operator. 


(Il) THe Vector PERFORMANCE OPERATOR 


In developing the concept of a vector performance 
operator, it is necessary to expand the brief description 
of complex plane quantities contained in Part (I). In 
the complex plane a quantity is composed of both real 
and imaginary numbers and, hence, can be represented 
by a vector extending from the origin of the axis system 
to a particular coordinate in the plane. This quantity 
has magnitude with respect to the origin and direction 
as measured from the real axis, thus exhibiting its vector 
characteristics. An additional characteristic assigned 
to the vector quantities used in this paper is rotation 
about the origin at constant rates. The time for one 
complete revolution (27 rad.) of a vector about the 
origin is called the forcing period, 7,; the rate in 
revolutions per second is the reciprocal of the forcing 
period and is called the forcing frequency, m, At any 
instant, the angular position of a vector quantity of 
this character depends upon the forcing frequency and 
upon the elapsed time, ¢, for in this time it will have ro- 
tated through n,t parts of a complete revolution. It is 
more convenient to express the position of the vector 
quantity in radians, in order to avoid the continual use 
of the conversion factor, 27. For this purpose, a new 
rate is defined called the angular forcing frequency, 
wy—that is, 


ao, = = (4) 
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The generalized forcing function of concern in this 
paper is the input quantity, g,,. Vectorially, this 
quantity has the direction w/t and a magnitude, or am- 
plitude, @¢n)2, as shown in Fig. 3. This vector input 
quantity is equal to the sum of its real and imaginary 
components, or 


Jin = Qcin)a COS wt + JQ SIN wt 


[COS wt + j sin w t | 


(5) 


The common practice, when projecting a vector onto 
X, Y, Z cartesian coordinates, is to make use of unit 
vectors that, by definition, have a magnitude of unity 
and are directed along the respective Y, Y, Z axes. In 
this manner a vector can be expressed in terms of unit 
vectors and the projected magnitudes of the respective 
vectors along the axes in question. This is usually ex- 
pressed as 


(6) 


The vector input quantity occurring in the complex 
plane is similarly expressed, as shown by Eq. (5). In 
this equation the unit vector along the imaginary axis is 
denoted by the symbol j; the unit vector along the 
real axis is implied, although for the sake of brevity it is 
not actually included. 

The merit of the vector quantity is more fully under- 
stood when Eq. (5) is expanded by the use of a Mac- 
laurin series, for then 


din = 1 + + 


A = iAx + jAy kAz 


dnd 1 + (jut) + 2! 3! 


(7) 


In this"equation the vector input quantity is expressed 
concisely as the product of its amplitude, g(n)2, and the 
exponential, e/*', The latter unit vector is rotated 
wt rad. from the realaxis. With this exponential form, 
the operations of differentiation and integration require 
only the multiplication or division of the vector quan- 
tity by jw, A comparison of Eqs. (5) and (7) shows 
that, when equals 1/2 rad. 


= 


Gin = = (8) 
From this, it can be seen that 
j= (9) 


Hence, multiplication by jw, changes the magnitude of 
the vector function in direct proportion to the value 
of w, and advances the vector 7/2 rad., or 90°, in a 
positive (counterclockwise) direction. Similarly, divi- 
sion by jw, changes the magnitude inversely with w, 
and retards the vector by 90°. 
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* vector performance operator of Eq. (3). 


The response of a component or system to an applied 
vector input yields a vector output quantity that, for 
steady-state conditions, rotates at the same rate or 

* frequency as the input. In general, however, the out- 
put is not in phase with the input but is advanced or 
retarded by a phase angle. The output vector quantity, 
with amplitude Q(our, is shown in Fig. 3 leading the 
input by the angle (PA) )[.¢.4.); therefore, it can be 
written 

+ (PA) | 


Gout Q(out)al@ (10) 
It is desirable to consider the operation on the vector 
input that is required to obtain the related output re- 
sponse. This operation, which affects both the magni- 
tude and the phase of the output, is specified by the 
This equation 
can be combined with Eqs. (7) and (10) to give 
j(PA) 


Eq. (11) brings out a further advantage in the use of 
this type of vector function—that is, it yields a perform- 
ance operator that is independent of time. However, 
the vector performance operator for a_ physically 
realizable component or system 7s a function of the 
forcing frequency. 

The relationship between the amplitude of the output 
and the amplitude of the input is simply expressed by 
the amplitude ratio, (AR) 

(AR) = (12) 
It is convenient to consider the amplitude ratio as the 
product of two factors: (1) a reference amplitude ratio 
(AR) 0 faingoul(ret)» Which is independent of the forcing 
frequency but which carries the dimensions of the am- 
plitude ratio; (2) a dynamic amplitude ratio 
(AR) ofeingoul(dyn)» Which depends upon the forcing 
frequency but which is dimensionless, In equation . 
form, 


(AR) teindout] = (AR) (ret) (AR) 0) fein gout) (13) 


The selection of the reference amplitude ratio de- 
pends upon the particular component or system under- 
going specification. With some components and sys- 
tems at zero frequency, when conditions are static, a 
given change in the level of the input yields a zero 
change in the output; with others, a proportional 
change is introduced; with the remaining, an infinite 
change is produced. An example of the latter is 4 
servomotor that, when excited, continues rotating until . 
it either reaches a limit stop or is de-energized. An- 
other example is an airplane with the elevator displace- 
ment as the input quantity and the pitch angle as the 
output quantity. The limit in this case is a function 
of air speed and gravitational effects. 
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When considering a component or system exhibiting 
either zero Or infinite response under static conditions, 
the choice of the reference amplitude ratio is arbitrary 
and should be selected for optimum convenience. 
However, when variations in the input produce propor- 
tional changes in the output, the component or system 
has a sensitivity, So fendouls equal to the ratio of the 
static output to the input. Whenever it exists, this 
sensitivity serves as the reference amplitude ratio. 
In this circumstance and when the frequency is zero, 
the dynamic amplitude ratio is equal to unity. In 
summarizing the developments of this section, it can 
be generally stated that a 


IPA) (Ky 
= (AR) 


(AR) (ret) X 


j(Pa) 
(AR) feingourl (dy 


§ Pho 


(14) 


and, when a finite response exists under static condi- 
tions, 


j(PA) 


(15) 


The sensitivity carries an algebraic sign (plus or minus) 
so that for static conditions the phase angle of Eq. (15) 
isequal to zero. It should be noted that the reference 
amplitude ratio of Eq. (14) cannot always be defined 
in this same manner, because the phase angle at zero 
forcing frequency may be either 90° or 270°. 

The vector quantities described in this paper are 
mathematical abstractions that help to establish a new 
approach to dynamic analysis. In addition, the vector 
performance operator can be used directly to obtain 
the response of a linear component or system to all 
periodic functions. A first step in obtaining this result 
is accomplished, as described in Part (III), by relating 
the periodic functions to vector quantities. 


(III) QUANTITIES 


(a) General 


A series of vector quantities can be used to represent 
any periodic real quantity, gp. The series consists of 
one nonrotating vector along the axis of reals and an 
infinite number of paired vectors—the two vectors of 
each pair rotating in opposite directions at angular 
frequencies that are integral multiples of the funda- 
mental forcing frequency wy. Stated mathematically, 
the series is 


For complete accuracy, all integral values of the 
number ”, including zero, must be used in Eq. (16). 
The series usually converges rapidly, so that values of 
" greater than 20 are seldom necessary to obtain prac- 
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Fic. 4. Relationship between simple harmonic function and 
equivalent vector quantities. 


tical engineering accuracy. The coefficient K, of 

the exponential time function in Eq. (16) has vector 

properties and can be evaluated by first multiplying 

each side of the equation by e~'"““ and then inte- 

grating the resulting equation over a time interval of 

one fundamental forcing period, T,. This gives 

Tyi/2 
K, = (1/Tp) (17) 
—Tp/2 

Two types of periodic quantities, the simple harmonic 

and the square wave, are of particular interest. The 

first type is useful when analyzing the steady-state 

response; the second, when analyzing the transient 
response. 


(b) Simple Harmonic Function 


The simple harmonic function is useful because it 
lends itself conveniently both to experimental test and 
to mathematical analysis. There are many methods 
in the laboratory for producing harmonic motion. 
Ball-and-disc integrators, Scotch yokes, and bell cranks 
are a few of the mechanisms that can be employed to 
obtain this type of displacement. Harmonic voltages 
can be generated by means of alternators, tuning forks, 
or oscillators. By selecting the proper time reference, 
all of the resulting quantities can be written: 

Qsh = Ja COS (18) 


The amplitude of the quantity g, is the maximum de- 


_ viation of the quantity from its mean value, which in 


this case is zero. The harmonic quantity completes a 
full cycle every time the angle w,t is advanced 360° 
or 2rrad. It can be seen that the time for a complete 
cycle, .called the forcing period T,, is related to the 
angular forcing frequency wy, by the same equation 
[Eq. (4)] that relates comparable parameters used to 
describe the vector input and output quantities. 
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QUANTITY 
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sy = SQUARE WAVE 
QUANTITY 


T, 


WHEN 


Thy 
| 


<t<(>-t), 
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EQUATIONS 16 AND 17 YIELD 
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Fic. 5. Relationship between square wave function and infinite 
series of vector quantities. 


The simple harmonic quantity g,, can be expanded 
into the vector series of Eq. (16) by substituting Eq. 
(18) in Eq. (17) and solving for the coefficient K,. 
When the frequency of the fundamental vector quantity 
wy in the series is set equal to the frequency of the 
harmonic motion wy. 


Ry = = (19) 
and 
R, = KR. = Ky =... = Rin = Ki-y = ... = 0 
(20) 
so that, by the use of Eq. (16), 
qn = + /2)* (21) 


Eq. (21) states that the simple harmonic function 
equals half the sum of two conjugate unit vectors, one 
rotating positively and the other negatively. These 
vectors are represented graphically, along with an alter- 
nate proof for the relationship of Eq. (21), in Fig. 4. 


(c) Square Wave Function 


The square wave has been used extensively by engi- 
neers for testing amplifiers and other electronic circuits. 
When the period of the square wave is long in compari- 
son with the transient response time of the circuit, the 
square wave response during each half cycle is similar 
to that of a transient—in effect, giving the latter re- 
sponse in a repetitive form that can be observed on an 
oscilloscope. 


* For certain functions, especially those that are nonsymmetri- 
cal with respect to the zero time axis, it is more convenient to use 
a sine function than a cosine function. The sine function lags 
the cosine function by +/2 rad., so that, when gens = ga Sin wyt, 
the use of Eqs. (9) and (21) shows that 


gen = gal (e947 — 


In Fig. 5, a generalized square wave function is shown 
pulsating between zero and a level designated as g, 
because the function usually has this value at zero time. 
The square wave function or quantity, 950, has a cyclic 
period, T,;, and is symmetrical in the respect that half 
the period is devoted to the zero level and half to the 
level go. In Fig. 5, the square wave function is leading 
the zero reference time by an amount f, or, stated math- 
ematically, 


> Qsw = Qo 
Ts Tn 
> Qsw = O 


The use of these conditions in Eq. (17) permits the 
evaluation of the mth square wave coefficient—i.e., 


1 (Ti/2) — ty 


Tn 


K n(se) = Qoe “dt 


(qo/2jnx)(1 — (99) 
This equation, when combined with Eq. (16), gives the 
following expression for the square wave quantity: 


(23) 

The square wave function is related in Eq. (23) toa 
series of vectors; consequently, when the input to a 
system can be represented by such a function, the out- 
put of the system can be obtained by application of the 
appropriate vector performance operator to each of the 
terms in the series. The amplitudes of the exponentials 
in Eq. (23) vary inversely with the order of the integer 
n and, in addition, depend upon whether 1 is zero, an 
odd integer, or an even integer; for when 


n = 0, (1/2jnr)(1 — e~ = 1/, 
(24) 
1 
n = odd integer, ——(1 — e*") = — (2%) 


(1 — =0 (26) 


n = even integer, 


With these three relationships, Eq. (23) can be con- 
verted to the form 


A discussion of the response of a system to a square 
wave input is contained in Part (IV). In Part (V) 
techniques are presented for establishing the vector 
performance operator of a component or system by 
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either mathematical or experimental means. It is also system, once displaced, never again reaches equilibrium 
shown in Part (V) that the simplest experimental but rather yields an error with progressively increasing 
method involves the application of long-period square magnitude. 


waves, which become step functions in the limit. This Stability is a necessary criterion for dynamic per- 
method is based on the theoretical considerations of formance, but equally important are the degree of sta- 
Eq. (27). . bility and the response time. The degree of stability is 
a measure of the remoteness of a system from the con- 

(IV) RESPONSE TO A SQUARE WAVE INPUT dition of instability; it is associated with the number of 


It may rightly be asked why stress is placed on the oscillations required to damp out a given transient 
response of a system to a square wave input. Inorder resulting froma step input. The response time is that 
to answer this question it is needful to consider what required for the transient to become reduced to a 
constitutes satisfactory performance under the actual stated fraction of its initial value. 
operating conditions of typical automatically controlled The automatic pilot experiences step inputs not only 
aircraft. For example, in most multiengined aircraft when the airplane encounters a gust but also when the 
used by the air lines, an automatic pilot relieves the pilot switches from manual to automatic controls and 
human pilot during the more routine phases of a flight. when he retrims the aircraft, with the automatic equip- 
The purpose of the automatic pilot is to hold the pitch ment engaged. The degree of stability and the re- 
and heading initiated by the pilot. When an airplane sponse time are important with each of these step in- 
encounters a gust while in flight, the nose may be forced puts. The classical theoretical determination of the 


-up or down; in either case the actual pitch of the air- response to a step input is laborious and is seldom used. 


plane becomes altered from the pitch value set by the However, the degree of stability and the response time 
pilot. In a satisfactory system, the elevator is then can readily be evaluated by the use of the square wave 
actuated by the servo controls to restore the pitch angle input, which, from Eq. (23) can be written 

to the pitch setting. A system that responds in this 

manner is stable as opposed to an unstable system, ‘> (38) 
which augments an initial transient error. An unstable n= —o 2jnn 


The response of a component or system to a square wave input can be determined (1) by operating individually 
on each of the vectors in the series of Eq. (28) and (2) by superimposing the results, giving 


n= @ 1 


The ability to superimpose the vector responses to each of the terms in the series implies the existence of linear equa- 
tions for the component or system. This same limitation is placed on most mathematical approaches, excepting 
step-by-step integration processes. To date, this technique has only been used for components and systems that 
have static sensitivities. In such cases, a vector performance operator can be replaced by the equivalent form of 
Eq. (15); therefore, 


does * (AR) (aya) Att + (PA) ] 
‘out. 1 —e e? Of¢ingout} 30 
© Gout] 2. ( ) ( ) 


However, according to the definition of static sensitivity given in Part (II), the dynamic amplitude ratio is unity 
and the phase angle is zero at zero frequency. Therefore, when Eqs. (24), (25), and (26) are substituted in Eq. 
(30), 


Govt + S0 > (AR) (472) me, + (PA)oteingoutt] (31) 
2 = 1,3,... jn 
—3,... 
2S AR yn). 
1,3,... 


The response of a component or system to a square series of Eq. (32). This series contains only odd har- 
wave forcing function can be obtained once the sensitiv- monics, so that, when the dynamic amplitude ratio is 
ity, dynamic amplitude ratio, and phase angle have equal to unity, the 19th harmonic, or tenth term, in the 
been determined by the summation of the Fourier series is '/;th, or approximately 5 per cent, of the. fun- 
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BROKEN LINE CURVES REPRESENT SQUARE WAVE INPUT AND THE SUBSEQUENT RESPONSE 


SOLID LINE CURVES REPRESENT STEP FUNCTION INPUT AND THE SUBSEQUENT RESPONSE 


WHEN + >> RT, THE SQUARE WAVE RESPONSE IS THE SAME AS THE STEP FUNCTION 


T 
RESPONSE FOR THE TIME INTERVAL O < ¢ < a 


Fic. 6. Comparison of the response of a system to a square 
wave input with the response to a step function input. 


damental. Usually, the amplitude ratio is consider- 
ably less than unity at the higher frequencies, so that 
the elimination of the 19th and higher harmonics incurs 
an‘inaccuracy of less than 2 per cent. Several methods 
for summing the Fourier series are discussed in Appen- 
dix B: These include a description of an electromechan- 
ical synthesizer that produces a graphical plot by os- 
cillographic means, according to the relationship of Eq. 
(32). In this synthesizer the phase angle and the am- 
plitude ratio of the first 17 harmonics (nine terms) can 
be introduced by the rotation of selsyn stators and the 
attenuation of the corresponding selsyn armature out- 
puts. 

In a stable component or system the transient in- 
accuracy resulting from art abrupt change in the input 
decreases and approaches zero asymptotically. Conse- 
quently, an infinite time is required to eliminate the in- 
accuracy: for all idealized components or systems that 
do not have the limitations of Coulomb friction and 
electrical hash level common to all practical units. 
With an eye to such considerations, the response time 
of a component or system is arbitrarily defined as the 
time for the transient inaccuracy to fall below, never 
again to exceed, 5 per cent of its initial value. 

In Fig. 6, the response to a step input is compared 
with the response to a square wave input. The half 
period of the square wave is selected enough greater 
than the response time so that at the termination of 
each half cycle the response is substantially quiescent. 
Therefore, at the beginning of each new half cycle, the 
response is nearly constant— a condition that also holds 
when time equals zero. For this reason, the compon- 
ent or system responds to the square wave input during 
the first half cycle as though a step input were applied, 
provided the half period is sufficiently longer than the 
response time. The validity of this heuristic reasoning 
is enhanced by the additional proof given in Appendix 
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A and is corroborated by numerous examples success- 
fully handled in this manner. 

The fact that the square wave response is nearly 
equal to the step response during the first half cycle of 
the square wave permits its use when determining both 

* the response time and the degree of stability. A low 
degree of stability in a system implies such nearness to 
the unstable region that a slight variation of system 
parameters might lead to instability. In comparison, 
a system with a high degree of stability can be con- 
structed with low tolerances on the system parameters 
without danger of instability. A measure of the degree 
of stability is the number of oscillations required to 
damp an initial transient. More definitely, the ratio of 
a transient output peak, Q(our)(m+ 1, to its preceding 
peak, Qout)(m), indicates the degree of stability, for, 
when this ratio equals unity, the system is on the bor- 
derline between instability and stability and, when it 
approaches zero, the degree of stability becomes large, 
The stability coefficient, SC, is a quantitative measure 
of the degree of stability and is related to the peak ratio 
by the expression 


or, in percentage, 
SC = {1 [Qcout)(m lav} 100 (34) 


From Eq. (34), it can be seen that the stability coef- 
ficient is negative for unstable systems, is zero for 
borderline systems, and is 100 per cent for systems 
reaching a steady-state value “dead beat’”’ with no over- 
shoot. For best results in evaluating the stability 
coefficient, several peak ratios should be averaged, as 
is indicated in Eqs. (33) and (34). 

Before obtaining the response of a system to a square 
wave input, it is necessary to obtain its vector per- 
formance operator. This entails a knowledge of the 
vector performance operators of each of the system 
components. Several methods for evaluating the com- 
ponent vector performance operators are discussed in 
Part (V). 


(V) MetHops FoR DETERMINING COMPONENT PER- 
FORMANCE OPERATORS 


(a) By Use of Performance Equations 


The performance operator of a component can be ob- 
tained either by computation or by direct measurement. 
The mathematical approach presupposes that the per- 


formance equations describing the characteristics of ° 


the component are available in accurate form. The 
performance equations offer the only methdd available 
when a component is in the design stage or when the ef- 
fects of anticipated design changes of a completed unit 
are being studied. 

The performance equations of a component are de- 
rived by summation of the pertinent physical quanti- 
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ties—e.g., forces, torques, currents, voltages. The re- | = (37) 
sulting equations representing the physical phenomena 
can then be expressed in the following general form, 
plus an additional number of performance equations in 
direct proportion to the number of output quantities: 


The vector outputs are related to the vector input by 
the use of the vector performance operator defined in 
Eq. (3); therefore, 


$P $0 (38) 


F +...) + (35) ang 
= re (36) when expressions for the input and output quantities 
and their derivatives are substituted in the perform- 


Although the number of performance equations de- ance equations [Eqs. (35) and (36)], each term con- 
pends in direct proportion upon the number of output tains the kernel, Gin)tse’“’. This kernel can be fac- 
quantities, which in this case are the dependent vari- tored from the equations, leaving only the vector per- 
ables, there is no limit to the number of input quantities. formance operators, the parameters of the component 
However, only one input is shown in Eqs. (35) and _ and the forcing frequency. It is possible to solve the 


.(36) because, with linear equations, the effects of the equations simultaneously for the vector performance 


various inputs can be superimposed and can be con- operators, one of which can be written as the quotient of 
sidered separately, according to Eq. (7), in the form determinants, or 


lL... 


(A R) O gi (4 0) 


Consequently, the amplitude and phase relationships of the vector performance operator are dependent upon the 
forcing frequency and the component parameters, according to Eq. (40). 

The accuracy involved in computing the performance of a component is a function of the accuracy to which 
the values of the component parameters are known. The accuracy in the case of the airplane is admittedly ap- 
proximate, because there is considerable latitude in the interpretation of certain of the stability derivatives. 
Several of these derivatives, such as dC,/da, dC»,/da, and dC,/dég, can be determined to a satisfactory degree of 
accuracy from wind-tunnel tests on large scale models, but other derivatives, such as the damping in pitch, dC,,/dg, 
require special tests and in most cases are not determined experimentally. Some of the derivatives are seriously 
affected by changes in the propeller thrust, but at present there is not available sufficient generalized data to com- 
pute with accuracy the effects of these power changes. 

In contrast to the statements in the preceding paragraph, the techniques introduced in this paper do not require 
knowledge of the individual stability derivatives when the performance operator is obtained from flight test. 
Measurements for this purpose can be made at any desired altitude, attitude, and speed within the limits of the 
aircraft, and such effects as Mach Number and propeller normal force variation with pitch angle are automatically 
included. The effects of lag in the build-up of aerodynamic forces and torques caused by accelerations that re- 
quire rapid changes in the flow pattern around the airplane are also included. 


(b) By Experimental Application of a Simple Harmonic 
The vector performance operator of a component can be evaluated experimentally by measuring the amplitude 
and phase of the response to a simple harmonic input of known amplitude and frequency. The input can be 


represented mathematically as. 


Jin = Qcinjq COS (41) 


or, by Eq. (21), 
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The response of the component can be determined by operating independently on the two conjugate vectors of 


Eq. (42), giving 


Qin)a just Qin)a 


2 


(AR) otemeoul {etter + (PA) olgingoutt] + } 


(AR) COS + (PA) } (43) 


As indicated by Eq. (43), the response to a simple har- 
monic input is itself a simple harmonic quantity, os- 
cillating at the same frequency as the forcing frequency. 
However, the response (the output), is not in phase with 
the input, and its maximum deviation from zero, called 
the amplitude, equals the product of the amplitude 
ratio and the input amplitude, so that 


Jout = Q(out)a cos { wt + (PA) } (44) 


The measurements that must be made to determine 
the amplitude ratio and the phase angle are shown in 
Fig. 7. A of Eq. (43) with Eq. (44) indi- 
cates that 


(A R) Ol¢ingoul = in)a (45) 


Therefore, the amplitude ratio can be obtained im- 
mediately by dividing the output amplitude by the in- 
put amplitude. Actually, a datum, or zero reference, 
line is not always available, in which case it is more con- 
venient to take the quotient of the double amplitudes. 
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t 
(PA) 360 (IN DEGREES) 


Use of a sinusoidal input for the experimental deter- 


Fic. 7. 
’ mination of the vector performance operator. 


In Fig. 7 the output response of the component is 
leading the input quantity by a time ¢,. The fractional 
lead equals the ratio of the time lead to the forcing 


period, but, since the forcing: period is equivalent to’ 


360°, the phase angle becomes 
(PA) = (4z/Ts) 360 (46) 


in degrees. 

The technique employed in carrying out these tests 
consists in establishing steady-state conditions by al- 
lowing sufficient time for the transient to disappear and 
then in recording by oscillographic or other means 
both the input and output quantities. The frequency 
of oscillation is changed in predetermined amounts 
until the desired frequency spectrum is obtained. For 
components of a commercial automatic control, the 
desired frequency range is from 0 to 10 rad. per sec. 
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Fic. 8. Diagram showing method used to obtain vector per- 
formance operator from step function data. 
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Military applications sometimes require response data 
at higher frequencies. 

The sinusoidal rate-of-pitch response for a U.S. Air 
Force A-26 medium bomber flying at an indicated air 
speed of 180 m.p.h. at 10,000 ft. altitude is presented 
in Fig. 10. In determining the associated vector per- 
formance operators the center of gravity is located at 
26 per cent of the mean aerodynamic chord. The 
maximum amplitude ratio occurs at a forcing frequency 
of 2.3 rad. per sec., which represents a forcing period of 
2.7 sec. This period corresponds to the short-period 
oscillation of the airplane. The long-period oscillation 
resulting from dynamic changes in air speed is greater 
by one order of magnitude. 

A year’s work is required to obtain the level flight 
jongitudinal and lateral performance operators by 
means of sinusoidal inputs when effects of air speed and 
center-of-gravity location are thoroughly investigated. 
For economic reasons, this period of time might be pro- 


_hibitive. In addition, for rapidly changing flight condi- 


tions resulting from turns, climbs, or dives, the time 
required to establish steady-state conditions is too long 
for the sinusoidal method to be effective. The step 
input has great value for these configurations, as well 


as for those in level flight, for the data is then collected 
in seconds rather than in hours. ; 


(c) By Experimental Application of a Step Input 


The response to a step input can be presented as a 
transient curve departing from its initial value at the 
instant the step input is applied. A typical transient 
response is illustrated in Fig. 8. In order to obtain the 
vector performance operator, this curve is approxi- 
mated by the series of steps indicated in the figure. 
For perfect matching of the actual response, infinitesi- 
mal steps are required, but for engineering purposes a 
choice of 15 to 20 steps gives satisfactory results. The 
time ordinate of the transient response, up to a point 
sufficiently beyond the response time to obtain the nec- 
essary accuracy, must be divided into even time incre- 
ments, At—the number of increments depending upon 
the number of steps required. It is necessary to meas- 
ure the change in output from the transient response 
during each of the time increments. These changes are 
denoted by the symbols 


dou ou 
af |. [2 “|. ete. 
Q(in)o 


The incremental steps can be expressed as functions of time by the use of the square wave quantity of Eq. (27). 
The square wave conforms to a step only for a time equal to the half period of the square wave, but no upper limit 


need be placed on this time duration. 
sponse, and for this reason 


Qcinyo in)o Tn = jn 


However, the sum of these steps is made to approximate the output re- 


af Jout 


1 


The time lead of Eq. (27) is replaced in Eq. (47) by a time delay. The time delay is selected so that most steps 


bracket the response curve that is to be matched. 


In this manner the area under the discontinuous approximating 


curve is kept as nearly as possible equal to the area under the actual response curve. 
However, the square wave response is expressed in Eq. (31) as a summation of vectors in a form similar to that 


of Eq. (47). When the period of the square wave response equals the period of the incremental steps, it is possi- 
ble to equate these two expressions and to solve for the vector performance operator, which becomes 
Soteatoul = — + — + --- (48) 
and 
(in)0 (in)O0 


Eq. (48) shows that the input-output sensitivity can 
be determined by the direct addition of the incremental 
steps. However, this sum equals the ratio of the 
steady-state output to the input, which can be meas- 
ured directly; hence, this equation serves only as an 
arithmetical check on the incremental values. The 


amplitude ratio and the phase angle can be evaluated 
by the vector addition indicated in Eq. (49). This ad- 
dition, which can be made either graphically or alge- 
braically, can be carried out for any forcing frequency 
wy. This forcing frequency replaces the product of the 
order of m and the fundamental frequency appearing in 
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Fic. 10. Comparison of rate of pitch vector performance oper- 
ator of airplane obtained from elevator step displacement with 
operator obtained from sinusoidal oscillation of elevator. 


Eqs. (27) and (47) because, in this case, both of these 
quantities are arbitrary—the only restriction being 
that the fundamental forcing frequency be small. 

The rate-of-pitch response of the airplane to a step 
displacement of the elevator is shown in Fig. 9. This 
response is found for the same modifying quantities of 
altitude, air speed, and center-of-gravity position as the 
sinusoidal response data already presented. For this 
reason it is possible to compare the vector performance 
operator obtained by sinusoidal input displacements 
with that obtained by a step displacement.: This 
comparison is demonstrated in Fig. 10, in which the 
circles represent individual sinusoidal measurements 
and the solid curve is a result of the step response. The 
agreement of the two sets of data is remarkably good, 
demonstrating the practicality of both techniques and 
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the fact that the airplane can be considered to be a linear 
component. 


(VI) ANALYsis OF MULTILOOP SYSTEMS 


Parts (I) through (V) furnish comprehensive back- 
ground for an analysis of multiloop systems and sug- 
gest the basic steps to be followed in designing a spe- 
cific system, These steps are summarized as follows: 

(1) Specify the input and the output, and consider 
the various systems that might give the desired rela- 
tionship between the two. Construct a functional dia- 
gram for each of the systems decided upon, and de- 
scribe each of the physical quantities connecting the 
components contained in the systems (the components 
must be physically realizable—for example, circuit, 
motor, shaft, airplane, or instrument). Define the 
relationship between the input and the output of each 
component in terms of the vector performance operator, 

(2) Express the system performance operator in 
terms of the component operators. This step is shown 
in Fig. 11 for a typical system. 

(3) Determine the values of the component vector 
performance operators by the methods introduced in 
Part (V). When the components are available, base 
the values on experimental test. When the equipment 
is in the design stage or otherwise unavailable, obtain 
the values from computations based on mathematical 
reasoning. 

(4) Effect a preliminary design by the use of ideal- 
ized components for which the less important dynamic 
effects are neglected (for example, the performance 
operator of an idealized derivative circuit is based on 
the assumption that perfect differentiation occurs at all 
forcing frequencies). The components should be ideal- 
ized to the extent that the overall system can be repre- 
sented by a third or fourth order equation; then, the 
methods of Routh can be used to ensure positive roots 
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and a stable system. A more complete description of 
the idealized performance operator is given in Appendix 
C. 
(5) Check the stability of the system, with all pos- 
sible dynamic effects included, by the use of polar plots 
of the vector performance operator. These plots are 
commonly called Nyquist diagrams. According to 
Nyquist criteria, a system is stable when the loci of all 
the open loop performance operators remain at the right 
of the minus unity coordinate on the complex plane 
and may be stable even though one of several loci passes 
to the left of the coordinate. The relationship between 
the input and the output of a typical single-loop system 
isshown in Fig. 11. The open-loop characteristic of the 
system can be obtained by measuring the variations 
in feedback for given correction signals, with the com- 
parator removed. The vector performance operator 
governing the action of this open loop is given by the 
product of the operators of the two components—e.g., 


_ the driving component and the feedback component. 


(6) Evaluate the closed-loop performance operator 
of the system by combining the component operators 
vectorially. Present the system performance operator 
as a function of forcing frequency in terms of amplitude 
ratio and phase angle. A further check on the dy- 
namic response of the system can be obtained by obser- 
vation of the performance operator presented in this 
form. For example, it can be shown that systems be- 
come unstable for which the phase angle becomes lead- 
ing as the frequency approaches an infinite value. 
Furthermore, the degree of stability of a system can be 
estimated qualitatively by evaluation of the amplitude 
ratio. When this ratio exceeds 1.3, the stability coef- 
ficient is too low for most practical systems. 

(7) Determine the square wave output of the system 
by introducing the closed-loop performance operator of 
(6) into Eq. (32). Quantitative values of the response 
time and the stability coefficient are then obtained 
from this square wave output by the method described 
in Part (IV). Care must be exercised when selecting the 
period of the square wave, for, with too short a period 
in comparison with the response time, steady-state 
conditions are not reached at the completion of each 
half cycle. In this case, the system does not have zero 
initial conditions at the start of each new cycle, and 
the pattern of the response is arbitrarily affected. 
Conversely, with too long a period, the fundamental 
forcing frequency is low, and an unnecessarily large 
number of terms is required for the Fourier series re- 
sponse to converge to reasonable limits. 

The distortion of the square wave response involves 
the so-called Gibbs’ phenomenon. Gibbs discovered 
that the summing of a Fourier series causes the graphi- 
cal representation of the sum of the series to peak at 
points of discontinuity or at points where abrupt 
changes occur. The phenomenon may be observed 
in the square wave response by a sudden deviation 
from steady-state conditions just prior to the end of a 


half cycle. This fluctuation may be minimized by ac- 
counting for more terms in the Fourier series; how- 
ever, even when the fluctuation amounts to as much as 
10 per cent of the steady-state level, the remaining sec- 
tions of the square wave response are practically un- 
affected. 

More serious distortion of the square wave response 
occurs when the amplitude ratio is inconsistent with 
the phase angle. When values of the amplitude ratio 
are known for all frequencies, the phase angle of a 
physically realizable stable system is uniquely defined. 
There can be several reasons for a discrepancy between 
the amplitude ratio and the phase angle. There is the 
possibility of a mistake in calculation, and there is the 
possibility of incurring an inaccuracy when measuring 
the component vector performance operators during 
experimental tests. Finally, there is a chance that an 
important frequency effect, such as a resonance peak, is 
being neglected because the upper frequency limit of the 
summation is too low. This distortion, as with the 
Gibbs’ phenomenon, gives rise to a deflection of the 
output response after the steady-state condition is 
reached but before the beginning of the next cycle. 


(VII) DeEsIGN or A Typical LONGITUDINAL CONTROL 
SYSTEM 


Aircraft control equipment must be sufficiently 
adaptable so that a given type unit will give satisfac- 
tory service in several different kinds of airplanes. With 
this provision, a manufacturer can supply a relatively 
economical product for both commercial and military 
use. However, the desired motion of the control sur- 
faces in correcting for a divergence of the aircraft from 
its intended flight direction depends upon the particular 
airplane and the use to which it is placed. Two system 
adjustments are provided with the equipment: one is 
called the proportional control, which governs the con- 
trol surface displacement for a given error signal; the 
other is called the derivative control, which governs the 
control surface displacement for a given rate of change 
of the signal. Asa specific case, the effect is considered 
here of the proportional and derivative adjustments for 
a typical longitudinal control system that includes a 
Minneapolis-Honeywell elevator servo and an A-26 
medium bomber. 

When in cross-country flight, an airplane is usually 
called upon to fly for protracted periods at constant al- 
titude. When checking his vertical flight direction the 
pilot utilizes such instruments as the altimeter and the 
rate of climb meter and must periodically retrim the air- 
plane to account for changes in angle of attack resulting 
from the consumption of fuel. When in level flight, 
the airplane must be pitched above the horizontal by 
an amount equal to the angle of attack; consequently, 
trimming the airplane requires a setting of the pitch 
angle, Oe. The difference between the pitch setting, 
Oset, and the correct airplane pitch angle for the in- 
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@ITCH ANGLE SETTING INTRODUCED (1)0, TOTAL PITCH INACCURACY OF 
/ AIRPLANE 
PLOT, Ov 


\ PITCH INACCURACY 
AIRPLANE PITCH ANGLE INDICATED 
BY VERTICAL GyRO, Da, (rye TOM “URACY OF VERTICAL 


Pitch oF ainecane, O, Xa, SRPLANE NC -LIFT LINE 


x 
X REFERENCE AMIS SET BY PILOT 
XX... VERTICAL GYRO HORIZONTAL 

V9" REFERENCE AXIS 


HORIZONTAL AXIS CONTAINED IN 
PLANE OF SYMMETRY 


CORRECT POSITION OF NO-LIFT 

CORRECT PITCH ANGLE LINE 

OF AIRPLANE, Osicorr) 


ALL LINES CONTAINED 'N PLANE OF SYVMETRY OF AIRPLANE 
ANGLES MEASURED PITCH OF AIRPLANE (Y-AXIS) 
ACCORDING TO NACA STANDARD CONVENTIONS 


EFINITION 
S PLOTS INACCURACY IN SETTING PITCH ANGLE 
set 


~ Paicorr) 


THEN 

(1)6, + (DO... + (T)6, 
WHERE 

CONSEQUENTLY 

INDICATED PITCH ANGLE CORRECTION 


= -(T},, = 


Fic. 12. Geometrical factors involved when controlling the 
pitch angle of an airplane. 


tended flight, @s corr), is the pilot’s setting inaccuracy 
(I) 6,-¢—that is, 


(Z) set = — Favcorr.) (50) 


A vertical gyro is used to provide a horizontal refer- 
ence within the aircraft, from which the pitch angle of 
the aircraft can be continuously measured. The ver- 
tical gyro indication of pitch angle, 4;, is inaccurate by 
an angle (J), ‘when the gyro is out of plumb with the 
earth, in which case 


(Z) = O41 — 


The angle between the horizontal reference established 
by the vertical gyro and the horizontal reference set 
by the pilot is the indication of pitch inaccuracy, (J) 4:. 
This indication is supplied to the servo controls to ini- 
tiate the required correction. The total inaccuracy 
in pitch, (I)@4, is equal to the sum of the inaccuracies 
introduced by the pilot, the vertical gyro, and the servo 
control. As seen from Fig. 12, 


(1) 64 = 64 — 94 (corr.) 


(51) 


The inaccuracy of the pilot setting is a random func. 
tion beyond the scope of all but psychologists, and the 
inaccuracy of the vertical gyro is small except when the 
airplane maneuver involves large horizontal accelera- 
tions; therefore, for the purposes of this paper, the 
total inaccuracy is assumed to be equal to the indicated 
inaccuracy. Fig. 13 isa functional diagram of the longj- 
tudinal control system. In this diagram, the compara- 
tor, in which the vertical gyro indication is matched 
with the pitch angle setting, is shown as a differential, 
The comparator has the dual function of indicating the 
inaccuracy and of transforming the inaccuracy into a 
suitable correction signal, (C)@4;, where 

= = — Oat (53) 

In the longitudinal control system the indicated 
pitch correction is the only input to the elevator com- 
puter. The function of the elevator computer is to 
provide the necessary proportional and derivative sig- 
nals for stabilization and control. These signals are 
transmitted to the elevator servo as the elevator com- 
puter output voltage. The elevator servo consists of 
an appropriate power amplifier, a motor, gear train, 
and output shaft. A followup on the output shaft 
generates a voltage that is matched with’ the input 
voltage from the elevator computer. Cable wrapped 
around a drum fastened to the servo output shaft is 
connected at its ends to the elevator horns. In this 
manner servo shaft displacements are transferred to the 
elevator, which governs the pitch of the aircraft. The 
input to the longitudinal control system can be con- 
sidered to be the pitch setting, and the output, the pitch 
angle, so that the system vector performance operator 
becomes 


$Pbucs) = = (AR)acsetP 


oy 


With the performance operator of the complete longi- 
tudinal control system established, the next step is to 
determine the performance operators of the several 


= + (1) (52) components. 
PITCH COMPARATOR 

ELEVATOR ELEVATOR PITCH 

PITCH INDICATED COMPUTER SERVO ANGLE 

PITCH ANGLE ANGLE / PITCH ANGLE output SHAFT ELEVATOR oF 
PILOT SETTING gf CORRECTION VOLTAGE ELEVATOR | DISPLACEMENT DISPLACEMENT AIRPLANE 

s 
INDICATED PITCH 
ANGLE OF AIRPLANE | VERTICAL ] 


Fic. 13. Functional diagram of aircraft longitudinal control system. 
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(a) Elevator Computer 


The elevator computer develops an output voltage, 
V,, equal to the sum of two voltages—one proportional 
to the indicated pitch correction (C)@4; and the other 
proportional to the rate of change of the pitch correc- 
tion—so that 


Ve = Seo(C) + Scal(C) (55) 


and 


$P}, V./(C)6a1 = Sep + 


2 
cp 


The elevator computer is designed to permit individual 
adjustment of the proportional sensitivity, S.,, and 
the derivative sensitivity, S.4. The magnitudes of 
these sensitivities determine the dynamic response of 
the longitudinal control system. 


Elevator Servo 


The elevator servo selected for the longitudinal con- 
trol system is part of the C-la Minneapolis-Honeywell 
autopilot. The elevator servo is alternately driving and 
braking when introducing an angular displacement, an 
action called pecking. The pecking action is sufficiently 
fast that the unit can be considered a linear compo- 
nent, and its performance can be expressed as 


{P}, = = = SAR) (57) 

An increase in the rate of pecking throttles the re- 
sponse of the elevator servo and increases its effective 
damping. The throttling can be adjusted within 
limits. In Fig. 14, which shows plots of the elevator 
servo amplitude and phase angle as functions of the 
computer voltage forcing frequency, the throttling is 
setata maximum. It can be seen that with maximum 
throttling the C-la elevator servo acts as though it 
were critically damped and had a natural frequency of 
approximately 7 rad. per sec. 


(c) Elevator Rigging 


The elevator is coupled to the servo output shaft 
with a stranded steel cable, which stretches in propor- 
tion to the reactive elevator load. The source of this 
load can be considered to be entirely aerodynamic, 
since the inertia of the elevator is sufficiently small in 
comparison with the elasticity of the cable so that the 
elevator natural frequency is greater than 60 rad. per 
sec. The aerodynamic load, called the hinge moment, 
restrains the elevator motion so that its displacement, 
for a given servo displacement, is always less when the 
aircraft is in the air than when it is on the ground. 
Therefore, 


(58) 
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Fic. 14. Vector performance operator of C-la elevator servo 
manufactured by Minneapolis-Honeywell. 


This equation brings out the important fact that the 
sensitivity of the elevator rigging, S,, is mot constant; 
it is a function not only of the rigging dimensions but 
also of the elevator surface area and the indicated air 
speed. 


(d) Airplane 


The experimental means used to determine the vector 
performance operator of the airplane is discussed in 
Part (V). The most economical timesaving method 
involves the measurement of the rate of pitch following 
a step displacement of the elevator. Note that this 
method cannot be used directly for pitch angle evalua- 
tion, because this quantity increases linearly with time 
following an elevator step input and does not reach 


‘the necessary constant steady-state value required 


when using the Fourier series methods of Fig. 8. How- 
ever, the pitch angle is the integral of the pitch rate; 
therefore, 


6 
oz 
Consequently, it can be seen by inspection 
(AR) a (50041 (AR) atsx41/s (60) 


and 


(PA) = — (4/2) (61) 


AIRCRAFT AUTOMATIC CONTROL || 
| 
1 the | 
1 the | | 
lera- 
the 
ated 
ara- 
ched 
i 
to a («le tan (56) 
-20 = 
(53) 
ated 
s to 
 Sig- 
ts of 
haft 
oped 
ft is 
this | 
) the 
The | 
con- 
i 
itch 
ator 
(54) | 
| 
is to 7 
eral 


550 JOURNAL OF THE AERONAUTICAL SCIENCES-SEPTEMBER, 1948 


It is shown by Eqs. (60) and (61) that the pitch 
angle performance operator of the airplane can be ob- 
tained from the pitch rate operator of Fig. 10 by divid- 
ing the pitch rate amplitude ratio by the forcing fre- 
quency to obtain the pitch amplitude ratio and by in- 
creasing the pitch rate phase angle by 7/2 rad. or 90° 
to obtain the pitch phase angle. The resulting data, 
originally derived from the pitch response to a step dis- 
placement of the elevator, are presented in Fig. 15. 


Bat 


SEP 


cp 


(e) Vertical Gyro 


The vertical gyro is essentially a long-period pendy. 
lum and is unaffected by most airplane maneuvers ex. 
cepting procedure turns lasting 1 or more min. Since 
this discussion is concerned with effects having a much 
shorter period, the vertical gyro output can be assumed 
equal to the input, in which case 


$ Pd = = 1 (62) 


The open-loop vector performance operator of the longitudinal control system equals the product of the vector 
performance operators of the computer, servo, rigging, airplane, and vertical gyro. 
of these five components are expressed by Eqs. (56), (57), (58), (59), and (62), respectively, from which 


The performance operators 


open-loop vector performance operator (transfer function) 


Soa\? j [wy Sed/Scp] S[(PA)s+(PA) 41626 
= 1+ («, e \ (AR) (63) 


to be designed 


In Eq. (63), the parameters associated with the ele- 
vator servo performance and the airplane appear 
within the braces marked ‘‘previously specified.” They 
are therefore not subject to the will of the designer. 
The only parameters remaining within his control ap- 
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Fic. 15. Angle of pitch vector performance operator for A-26C 
No. 898 at 180 m.p.h., IAS, with c.g. at 26 per cent M.A.C. 


previously specified 


pear within the alternate braces. They include the 
loop gain, which is equal in magnitude to the product of 
the computer; servo and rigging sensitivities, S,,5,S,; 
and the ratio of the elevator computer derivative sensi- 
tivity to its proportional sensitivity, S.¢/S-», which has 
the dimension of time. The influence of these param- 
eters can be judged by referring to the open-loop 
vector performance loci of the longitudinal control sys- 
tem contained in Fig. 16. These loci are evaluated 
from Eq. (63) and show that the system is stable when 
there is no derivative control and when there is a gain 
of one.* When the gain is increased to three, the locus 
falls on the left of the minus unity coordinate, indicat- 
ing an unstable control system. With the introduction 
of a computer derivative signal, both loci are rotated 
positively about the origin, thereby increasing the 
stability of the system. 

The longitudinal control system can be more accur- 
ately evaluated by using the relationships of Eq. (54) 
to compute its closed-loop vector performance operator 
in the form of amplitude ratio and phase angle. The 
results presented in Fig. 17 give an easily understand- 
able physical representation of the frequency spectrum; 
in each case it can be seen that the static performance is 
ideal, having an amplitude ratio of unity at zero fre- 
quency. However, the dynamic performances of the 
four cases presented differ widely. The resonance 
peaks occur at a higher frequency when the gain is in- 
creased, since this increase in proportional feedback ef- 

* The product of the sensitivites, S->S,S,, must actually havea 


negative value because, by N.A.C.A. convention, a positive ele- 
vator displacement gives a negative change in pitch angle. 
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Fic. 16. Open loop vector performance loci of aircraft longi- 
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Fic. 17. Vector performance operator of aircraft longitudinal 


control system. 


fectively stiffens the system. The increase in gain de- 
creases the stability, as evidenced by the height of the 
resonance peaks and by the leading phase angle ob- 
tained at forcing frequencies greater than 3 rad. per 
sec. The use of derivative feedback decreases the 
height of the resonance peaks and thus gives evidence 
of a greater degree of stability. 

The square wave response, shown in Fig. 18, is ob- 
tained by combining the data of Fig. 17 in the relation- 
ship of Eq. (82). Note that the response with no 
derivative and a gain of three has peak values that in- 
crease exponentially with time. This response does 
not reach equilibrium during each half cycle and, hence, 
fails to commence each new cycle with initial conditions 
equal to zero. 

The response time, which is the time for the transient 
to decay to 5 per cent of its initial value, can be evalu- 
ated directly from the response curve; the stability 
coefficient can be determined by the use of Eq. (34). 
These results are listed in Table 1. 

The system with a proportional factor of one and a 
derivative of 0.3 sec. clearly has the best characteris- 
tics, since the response time is the shortest and the de- 
gree of stability is the highest. It is interesting to note 
the appearance in this case of an exponential in addition 
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TABLE 1 


Effect of Proportional and Derivative Control on Response Time 
and Stability Coefficient 


Derivative Response Stability 
Proportional Factor, Time, Coefficient, 
Factor, ca/ Sep 
Sep Sede Sec. Sec. Per Cent 
-1 0 6.2 40 
—-1 0.3 4.0 65 
0.3" 5.9 26 


to the damped sinusoid. A still further increase in 
derivative control or a decrease in proportional control 
increases the characteristic time of this exponential, 
thereby yielding a longer response time. 


CONCLUSIONS 


(1) The use of the performance operator and the 
square wave response presents a method of analysis of 
systems that have (a) an unrestricted number of de- 
grees of freedom because the differential equations of 
motion do not have to be solved and (b) any number or 
any degree of complexity of loops. 

(2) The method permits the direct use of test re- 
sults, when they are available, without the need of re- 
ducing the test data to parameters defined by an as- 
sumed differential equation. 

(3) The method allows each component of a system 
to be considered separately and the action of the loop 
to be obtained by a simple combination of the indi- 
vidual parameters. 

(4) The powerful tools of the servomechanism engi- 
neer—for example, the Nyquist criterion—can be ap- 
plied without recourse to arbitrary and empirical rules 
with regard to the transient. 

(5) In the experimental determination of the per- 
formance operator of the airplane, effects such as aero- 
dynamic lag, compressibility phenomena, and power 
effects are included which cannot now be computed 
with accuracy. 


Appendix A 


USE OF THE FOURIER INTEGRAL FOR DETERMINING THE 
STEP-FUNCTION RESPONSE 


It is generally accepted that the study of the transient 
response of an automatically controlled system to a 
step function is the most effective means for its 
proper design and adjustment. However, the methods 
hitherto available have been difficult to use with all 
except simple systems having a limited number of de- 
grees of freedom. The vector performance operator in 
a comparatively simple manner provides the steady- 
state sinusoidal response of any system. Hence, in 
studying automatic control systems, a suitable means 
for converting the steady-state sinusoidal response to 
the step-function response would be a practical tool. 


1948 


The mathematical bridge between the steady-state 
sinusoidal response and the transient motion is the 
Fourier integral. Except for the most simple systems, 
however, its use in a formal analytical solution of 
transients is extremely difficult because the necessary 
definite integrals are not easily obtainable. Never. 
theless, the existence of the Fourier integral and the 
Fourier transform provides the fundamental idea that 
any transient input motion can be expressed as a con. 
tinuous function of frequencies and that the ratio of 
any frequency component in the input wave to the cor- 
responding frequency component in the output wave js 
determined by the steady-state parameters of the sys. 
tem at that frequency. This idea is not new and js 
probably the most basic principle in electric communica- 
tions theory, even though the analytical use of this 
principle has not been widespread because of the inte- 
gration difficulties mentioned. 

The developments of Parts (I) and (II) furnish a 
general method for establishing a step-function analy- 
sis, based on the Fourier integral. From the definition 
of the performance operator [Eq. (1)], 


out = (64) 


By extending the reasoning of Part (IV), which carries 
out the operation of Eq. (64) for the square wave input 
function, the step-function response can be obtained 
by transforming the step input function into the fre- 
quency domain. The theory of the Fourier transform 
permits the expression of any arbitrary time function 
as a continuous function of frequency, as follows: 


Gin = (1/24) (65) 
and 
qin = (66) 


In these equations, gin is the Fourier transform of the 
input, representing the amplitude and phase of a single 
harmonic component entirely analogous to the complex 
coefficient of the Fourier series. In the latter, the 
Fourier coefficients are defined only for the integer 
values of frequency and have no definition between 
integers; a plot of these values against frequency 
would consist of ordinates erected at each integer, such 
plots being called step line spectra. The interval of 
frequency is of differential order in the Fourier trans- 
formation representation, and gi. can be plotted as a 
continuous function of wy. 
For a step-function input, where 
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Gin = [from Eq. (66)] (67) 
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The integral of Eq. (67) is not convergent; therefore, 
, mathematical artifice must be employed, as follows: 


This equation is equivalent to Eq. (67). By this means 
the response of a system characterized by the vector 

ormance operator ¢P} can be obtained by use of 
Eq. (66), as follows: 


(69) 


=— 


‘In general, whether the lim,+» can be taken before 


integration or after is dependent upon the nature of 
§P}. However, in all practical cases, the limit can be 
taken prior to the integration of Eq. (70). Hence, 


* (71) 
Qc in)o Wy 
The vector performance operator ¢P} has both ampli- 
tude and phase, as shown in Part (II). Hence, 


ou 1 AR i(wrt+ 


ino Wy 


Expanding Eq. (72), 


AR 
cos (ws + + 


ao 


A 
j sin (wt + PA)dw, (73) 
J 


The first integral on the right-hand side of Eq. (73) 
has a singularity point at w, = 0 and therefore requires 
special treatment. Guillemin* has carried out the com- 
plex integration of a similar integral which can be ap- 
plied to the first integral of Eq. (73). Using the knowl- 
edge that for all practicalsystems, when w, = 0, PA = 0, 
and AR = S, the first ategeat can be evaluated as 
follows: 

Dai — cos (wt + ) yes (74) 


* Communication Networks, Vol. II, pp. 471 and 472; John 
Wiley and Sons, Inc., 1935. 


In the second integral of Eq. (73), the j's cancel, and, 
because w,and PA are odd functions and AR is an even 
function, the integral can be written as follows: 


sin (wet ~ PA) )\dwy = 


l AR 
— sin (wi + PA)dw, (75) 
T Wr 

0 


Therefore, Eq. (73) becomes 


out 1 AR 
= —sin + PA)dw, (76) 
Wy 


Qcinyo 2 


The response of a system to a step function can be 
evaluated by either analytical or graphical use of Eq. 
(76). The computational work in solving Eq. (76) ap- 
pears to be rather tedious, whereas the analytical 
solution may be difficult, especially when the vec- 
tor performance operator is obtained only in curve 
form. 

A numerical procedure for evaluating Eq. (76) can be 
obtained by changing the integration to a summation. 
As in most numerical schemes, the infinitesimal value 
of dw, is replaced by a small but finite value Aw,. 
Thus 


1 AR 
Bsn (wt + PA)dw, = 
Wy 


a=@ 


Tn=1 Wyn 


(wynt + PA) Aw, 


Let 
Aw; = 2w Bin = 


where wy, is a reference frequency. Then 


1 AR 
— sin (wp + PA)dw,; = 
Wy 


tbo 


sin + PA) (77) 
Wn=1,3,5... Nn 
or, finally, the response of a system to a step function 
can be approximated by 

Sets? 

2 wWn=1,3,5... 

Eq. (78) is precisely the same as that derived in Part 
(IV) for the response of a system to a square wave. 
The wy, thus takes on a physical meaning— it is the 
angular velocity of the basic harmonic of the square 
wave. This result might have been anticipated, in 
that the Fourier integral is the limit of the Fourier 
series as the fundamental frequency w,, is reduced to an 
infinitesimal. 
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Appendix B 


TECHNIQUES FOR SUMMATION OF FOURIER SERIES 


The performance of an automatically controlled air- 
craft when subjected to input disturbances is important 
because these disturbances (mathematically, step func- 
tions) are repeatedly applied to the aircraft in the form 
of gusts, pilot trim adjustments, and small excitations 
generated within the automatic equipment because of 
its minor imperfections. A system similar to an air- 
plane has two fundamental properties that govern its 
response to a step function—namely, its degree of sta- 


Jout S0 taindout] + 2S 0) tain gout] 
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bility and its speed of response. These properties cap 
be specified quantitatively in terms of the stability 
coefficient and the response time parameters defined 
in Part (IV). In designing a system, values of the 
stability coefficient and the response time should be ob. 
tained, so that the system can be compared with other 
systems constructed for the same task. 

It is shown in Part (IV), and confirmed in Appendix 
A, that these parameters can be determined from the 
response of a system to a square wave input of suf. 
ficiently long period. Furthermore, the square waye 
response can be related to the vector performance oper. 
ator by use of the equation 


sin [nwpt + (PA) (79) 


n = 1,3,5.. 


2 us 


The success of this approach hinges on the convenience 
of the method used to effect the summation of the 
Fourier series of Eq. (79). 

In one method of summation, values for ¢ are sub- 
stituted into Eq. (79), using sufficient terms for the de- 
sired accuracy. The summation performed in this 


manner can be arranged in a tabular form suitable for 
This 


computation by calculating machine operators. 
method is tedious, but can be carried out to any degree 
of accuracy; consequently, it is occasionally used to 
check results obtained by other techniques. 

The square wave response can be obtained more 
rapidly by using graphical methods for summing the 
Fourier series. One possible graphical scheme, de- 
scribed by Bromberg,* involves the use of ‘‘transient 
analysis charts.’’ In any graphical scheme each sinu- 
soidal term is represented as a vector quantity in the 
complex plane with magnitude 


[Gin Gout] (A R) O [ain @out} (dy 9) 


and with angle displacement with respect to the axis of 
reals 


[nwpt + (PA) 0 


The position of each vector can be quickly plotted for 
the desired times by means of the transient analysis 
chart or, when this is not available, by means of com- 
passes and protractor. However, the imaginary parts 
of the vectors are equal to the sine functions of Eq. 
(79); therefore, by means of a T-square and dividers, 
the vectors can be projected onto the axis of imaginaries 
and the indicated summation can be effected. 

The square wave response can be determined most 
conveniently by means of an electromechanical synthe- 


* Bromberg, B. G., Graphical A pplication of Fourier Series to 
Automatic Tracking Systems, Thesis submitted for Sc.D., M.I.T. 
1947. 
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sizer, one of which has been constructed by the Instruy- 
mentation Laboratory, M.I.T. In this device each 
harmonic is represented by a synchro as shown in 
Fig. 19. The synchro rotor is driven by a gear train 
so constructed that the speed ratio of any synchro with 
respect to the speed of the fundamental synchro equals 
the particular harmonic number. Nine synchros are 
available in the synthesizer, which makes it possible to 
take into account the 17 odd harmonics. The angular 
position of each synchro housing can be set to give the 
proper phase for the particular harmonic, and the out- 
put from each synchro can be attenuated by means ofa 
potentiometer to give the proper magnitude. The re- 
sulting output voltages are transmitted to Consolidated 
Engineering Corporation recording equipment, where 
the signals are amplified, rectified, and summed. This 
summation is impressed on a galvanometer in order to 
obtain a curve of the square wave response on oscillo- 
graph paper. 


AMPLIFIERS AND 
RECTIFIERS 


RECORDING 
OSCILLOGRAPH 


PHASE ANGLE ADJUSTMENTS 
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Square wave synthesizer with associated recording 


Fic. 19. 
equipment. 
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Appendix C 


THE IDEAL PERFORMANCE OPERATOR 


When analyzing compound systems, it is necessary 
in some cases and expedient in others to idealize certain 
of the components included in the system. In a dy- 
namic analysis the effect of Coulomb friction and elec- 
trical hash on mechanical and electrical components is 
yncertain and is usually neglected. This idealization 
js necessary and justified because, for all practical and 
usable components, these effects are minimized by care- 
ful design. Derivative and proportional circuits are 
provided in the elevator.computer described in Part 
(VII) to stabilize and stiffen the longitudinal control 
system. Most derivative circuits provide exact differ- 
entiation of their inputs over only a limited frequency 
band. On the basis of the analysis of Part (VII), the 
electronic engineer can be supplied with frequency 
specifications for the design of the elevator computer 


‘that will justify the assumption of perfect differentia- 


tion. 
It is sometimes necessary to idealize components to 


an extent known not to be realizable in actual practice 
in order to simplify the performance equations of the 
system for a preliminary analysis, as mentioned in (4) 
of Part (VI). The simplification of the performance 
equations should at least permit generalized conditions 
for stability to be established by the methods of Routh 
and, if possible, should permit a complete solution of 
the performance equations. When a component is 
knowingly idealized for these purposes, an additional 
subscript is affixed to the performance symbols, in 
which case Eq. (14) becomes 


faingout) “deal (AR) 0) (80) 


For the simplified analysis to be significant, the dif- 
ference between the actual performance operator and 
the idealized performance operator must be relatively 
small; consequently, when including the effects ne- 
glected in the idealization, it is convenient to express 
this difference as the alteration of the idealized per- 
formance operator necessary to obtain the actual per- 
formance operator. Changes in amplitude can be rep- 
resented by a performance ratio, PR, and changes in 
phase angle can be represented as a deviation, D, or 


(PR) (AR) 0 faingout] 
qingouty] 
(A R) [ginQout] deal) 


(81) 


and 


(D)(PA) ()taing (PA) ()faingout) (PA) (ident) 
(82) 
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Eqs. (81) and (82) permit the vector performance opera- 
tor of a component to be written in the form 


$0 


D j(D) ‘inJou 


In this manner the actual vecter performance operator 
is related to the ideal vector performance operator, and 
it becomes a simple matter to include effects originally 
neglected in the analysis. 
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SUMMARY 


The great increase in activity in the aircraft gas-turbine field in 
recent years has focused attention on problems associated with 
the attainment of very high efficiency of the various turbine com- 
ponents. The expansion or turbine component has had its share 
of attention from many research and development engineers. 
Since design work is progressing at such a rapid rate, it is felt 
that it would be helpful for all concerned to pause at various 
points in the progress to summarize present knowledge, evaluate 
the progress thus far obtained, and enumerate the problems still 
confronting us. In this paper, the authors have tried to present 
their views of such a summary. 


INTRODUCTION 


F  onaeesy THE AIRCRAFT GAS TURBINE is an out- 
growth of the stationary gas turbine and steam 
turbine, it has peculiarities that make it a different 
and more sensitive animal. It differs from stationary 
turbine power plants mainly by the fact that extremely 
high cycle efficiencies must be obtained with a much 
lighter weight power plant than the stationary turbine 
types. This requirement has a pronounced effect upon 
the expansion unit of the turbine, because it means 
that the same amount of work must be accomplished 
with a fewer number of stages. For example, a two- 
stage aircraft turbine may operate with a pressure 
ratio of 5 to 1 or more, while a stationary turbine oper- 
ating with the same pressure ratio may use eight or 
more stages. These fewer stages, must therefore, work 
a lot harder, which means that velocities or Mach 
Numbers, turning angles, tip speeds, and stresses must 
be much greater. Tip speeds of the order of 1,300 
ft. per sec. are used as compared to 700 to 800 ft. per 
sec. for a stationary gas turbine. Velocities slightly 
over sonic in the first stator and near sonic in the re- 
maining nozzles are required. Axial leaving velocities 
are usually above 600 ft. per sec. for the aircraft turbine, 
_whereas they are usually below this figure for the sta- 
tionary power plant. Rotor turning angles of the 
order of 135° are sometimes used in aircraft turbines, 
whereas those of stationary turbines are generally 
nearer 100°. 

Blades are usually of the reaction type rather than 
the impulse type, such as generally used in a stationary 
turbine. 


Received March 8, 1948. Based on a paper titled ‘“‘Aero- 
Thermodynamic Aspects of Aircraft Turbine Expansion Unit 
Design,” which was presented at the National Aircraft Propulsion 
Meeting, I.A.S., Cleveland, March 28, 1947. 
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The higher tip speeds, stresses, and lightweight 
construction usually require greater radial clearances 
than those obtained in a stationary power plant. This 
tends to increase leakage losses. The requirement for 
extremely high work per stage is usually so much 
greater compared to the tip speed allowed that the 
U/C is lower than that obtained in stationary turbine 
practice, and this means that the relative Mach Num- 
bers entering the rotor blades are considerably higher, 
With such high Mach Numbers and the greater turning 
angles required, the design of the blade shape and pass- 
age for high efficiency is much more difficult. 

The greater work per stage or fewer stages to do the 

same work means that the blade heights tend to in- 
crease at a faster rate than those for a stationary tur- 
bine, which introduces problems associated with flow 
over a rapid divergence of the inner and outer contours 
of the annulus. In order to keep this divergence from 
being too extreme, a tendency in aircraft turbine prac- 
tice is to increase the axial velocity from the first to the 
last stage. Thus, the axial velocity leaving the last 
stage tends to be higher than it would be in a stationary 
turbine. These factors tend to make the last rotor 
blade of the aircraft gas turbine shorter than the blade 
for a comparable job in the stationary turbine, but, be- 
cause it operates at a higher tip speed and does more 
work, the inlet temperature to it is greater. Therefore 
centrifugal, vibration, and bending stresses can be more 
critical. 

Another characteristic peculiar to the aircraft gas 
turbine which makes it more difficult to design is the 
fact that it must be efficient under much more varied 
operating conditions than the stationary turbine. For 
instance, even at constant r.p.m. the aircraft turbine 
must operate under a great number of different alti- 
tudes and air-speed flight conditions. These changes 
in temperature, pressure, air speed, etc., result in a 
considerable variation in pressure ratio across the ex- 
pansion unit at substantially the same engine r.p.m. 
Thus, a wide variation in Mach Number and in angle of 
attack on the blade is encountered, particularly when 
the number of stages has been held to a minimum. 

As newer and better high-temperature materials 
are developed, it becomes possible to improve the gas 
turbine cycle efficiency by increased combustion cham- 
ber outlet temperature. This tends to require aif- of 


liquid-cooled rotor blades, which introduce compromises 
The requirement for 


in blade shape and construction. 
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low weight dictates short chord length of blades, which 
introduces Reynolds Number problems as well as dif- 
{culties in design of blade shapes. 

All these differences peculiar to the aircraft turbine 
mean that knowledge gained from stationary turbine 
practice is only the beginning and that a great many 
problems remain to be solved in order to develop an air- 
craft turbine to its highest efficiency. 


DESIGN PROCEDURE 


J) Design Conditions and Assumptions 


Before beginning the design of the turbine expansion 
component, the following factors must be known: tur- 
bine expansion unit inlet temperature and pressure, 
weight flow, discharge static pressure, desired shaft 
output and turbine efficiency, maximum permissible 
wheel speed and maximum diameter, and r.p.m. range. 
From this preliminary information the’ designer can 
then determine the general dimensions and characteris- 
tics of the machine. 

The number of stages should be kept as small as 
possible because of weight considerations but should 
not be so small as to affect markedly the efficiency. 
Pressure ratios per stage as high as 3:1 and 4:1 are 
now used without seriously lowering efficiency. 


ID Velocity Diagram 

Having decided upon the number of turbine stages, 
the percentage of work done by each stage must be 
selected. In general, the first stages are designed for a 
greater amount of work than the later stages for the 
following reasons: 

(a) Higher temperatures in the first stages permit 
higher velocities (more work) without exceeding Mach 
Number limitations. 

(b) Aircraft turbines operate under widely varying 
pressure ratios. Since the incidence angles and Mach 
Numbers of the later stages are more affected by overall 
pressure ratio changes than the initial stages, the less 
work done in the later stages, the less will be the effect 
on overall power output under changed conditions. 

(c) Work done in any stage is, in the final analysis, 
afunction or rotor turning angle, and this angle should 
be kept smaller in the last stage to provide a low leaving 
velocity in a direction as close to axial as possible. 

Because of these factors, the Mach Number for each 
nozzle is made less than that for the preceding nozzle. 

In the design of most machines, the assumption is 
made that total energy of the fluid is constant, root to 
tip. With this assumption, several designs with respect 
to basic rotational airflow may be used, among which 
are: 

(a) Free vortex, in which tangential components of 
absolute velocity vary inversely with radius from root 
to tip and axial velocity is constant at all radii. 

(b) Constant mass flow per unit area at all points 
from root to tip (oVaxiai = Constant). In this type of 


diagram the maximum Mach Numbers may be reduced 
for the same work. 

(c) Solid rotation in which tangential components 
of absolute velocity vary directly as the radius from 
root to tip. In this type of diagram the axial velocity 
varies greatly from root to tip. 

With the exception of the free vortex type, the calcu- 
lation of a velocity diagram involves considerable trial 
and error, since the axial velocity is not constant root 
to tip. A radial pressure gradient is necessary in all 
types of diagram to balance the centrifugal forces aris- 
ing from the absolute whirl velocities. For equal work 
root to tip, the vele-ity vectors at one radius, coupled 
with the centrifugal pressure gradient, determine the 
density and absolute velocity at any other radius. 
Since the absolute whirl component at any radius de- 
pends on the centrifugal pressure gradient, the axial 
velocity at any radius may be determined. All these 
velocities at both stator and rotor exit must be balanced 
so that the stage work root to tip is equal. Free vortex, 
though not necessarily the most efficient type of design, 
permits simple geometrical determination of all veloci- 
ties at any radius, having established a particular veloc- 
ity at a particular radius. The free vortex type of de- 
sign has, therefore, been chosen as a basis for discussing 
design procedure. 

Having decided upon the type of design, the vector 
diagram angles and velocities are calculated as follows: 

(a) The approximate first stator nozzle area and 
blade height are determined by means of a conventional 
nozzle formula,using the given weight flow, inlet pres- 
sure and temperature conditions, pressure ratio across 
the nozzle, and a stator exit nozzle angle. Maximum 
diameter ratios of the order of 0.85 are now in use. 
From the approximate first stator blade height, root 
and tip design radii that will be within the physical 
blade height dimensions are next selected. 

(b) Since high pressure drop through the first stator 
is desirable for previously mentioned reasons, a stator 
root velocity close to sonic is assumed. The use of 
highly supersonic velocity, over 1.2 Mach Number, is 
undesirable, since a purely convergent nozzle is in- 
efficient above this Mach Number and a convergent- 
divergent nozzle above this Mach Number has too 
critical performance under variations in operating con- 
ditions. 

In order to determine Mach Number, the static en- 
thalpy, /72,, must of course be known and is determined 
from equation 


He, = H, — (V2/2gJ) 


where //, is the stator entering total enthalpy, B.t.u. per 
Ib., and V is the actual stator discharge velocity, ft. per 
sec. The value of H2, thus calculated is valid for any 
nozzle efficiency, since any energy not converted into 
velocity must, from the energy conservation law, still 
exist in the form of heat (temperature). Having H.,, the 
static temperature 7», may be ascertained from tables, 
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such as Thermodynamic Properties of Air by Keenan 
and Kaye. 

The discharge pressure required for generation of the 
desired velocity is a function of the nozzle efficiency 
and may be determined as follows: 


Horg = Hy, — (V?/2gJE) 


where [27x is the theoretical discharge enthalpy and £ 
is the assumed nozzle efficiency, which, for careful de- 
sign, is close to 92 per cent. The pressure correspond- 
ing to this theoretical enthalpy is then determined from 
air tables or, more laboriously, by application of the 
adiabatic gas law. 

(c) The work in each stage per pound of fluid is di- 
rectly proportional to LU’ AU where 


U 
AU 


rotor speed in ft. per sec at any diameter 

algebraic difference of tangential components 
of the rotor absolute entering and leaving 
velocities at the diameter in question 


The sum of U and AU(U + AU) forms the base of 
the velocity diagram at any diameter. (See velocity 
diagram of Fig. 1.) This base is also equal to the sum 

_ of the tangential component of the stator leaving veloc- 
ity and the tangential component of the rotor relative 
leaving velocity. In the design calculations these two 
tangential velocity components are so proportioned that 
the rule of decreasing Mach Number for successive noz- 
zles is maintained. At the same time, the axial velocity 
of the rotor is maintained equal to, or slightly greater 
than, the axial velocity in the stator. These factors de- 
termine the axial velocities in the stator’and rotor and 
the angles leaving the stator and rotor and should not 
be allowed to result in a rotor relative inlet Mach Num- 
ber greater than 0.70. 

(d) The tip diagram is determined from the root 
diagram for a free vortex machine by use of the fact 
that axial velocity root to tip is a constant and Vr = 
constant, where 


V. = whirl or tangential component of absolute 
velocity 
r = radius 


(e) The static enthalpy and specific volume leaving 
the stator were found in (b) above. The static enthalpy 
and specific volume entering the rotor are the same as 
that leaving the stator, since the process of moving onto 
the wheel has not changed them. 

(f) The static enthalpy leaving the rotor is deter- 
mined from the fact that the difference in static en- 
thalpy from entering to leaving conditions is equal to 
the energy difference represented by the squares of the 
relative velocities. As before, the static temperature is 
determined on an adiabatic basis directly from the 
static enthalpy and air tables, while the pressure is 
dependent upon the nozzle efficiency. 

(g) Knowing the specific volume and axial velocity 
at root and tip sections and the weight flow required, 
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the annular area and blade height required at the stator 
and rotor exit is calculated by means of the continuity 
equation W/A*= V/v (weight flow per unit area equals 
ratio of velocity to specific volume). In the same man- 
ner, using the same specific volumes but the relative 
exit velocities, the nozzle throat area required is cal- 
culated. 

(h) The same process is followed in successive 
stages with velocities leaving the nozzles being assumed 
on the basis of decreasing Mach Number in each succes- 
sive nozzle. Also the axial velocities are allowed to in- 
crease in order to keep the increase in blade height as 
small as possible. 


(i) In the last stage of the machine an attempt is 
made to proportion the stator and rotor enthalpy drops 
so that the absolute discharge velocity is as nearly axial 
as possible, since the whirl component of this velocity 
is considered useless even in a jet machine. 


(j) Fig. 1 shows the vector diagram of a free vortex- 
type turbine such as might be used to supply compres- 
sor power for a turbojet machine. For those interested, 
detailed calculations for this machine appear in the Ap- 
pendix. It will be seen that the discharge velocity of 
the machine is not axial. The design, however, is 
based upon sea-level conditions, and the increase in 
overall pressure ratio resulting from operation at speci- 
fied altitude is such that this velocity shifts to a nearly 
axial position. 


(UI) Blade Profile Shape 


Having determined the required blade angles, root 
and tip, for both stator and rotor, the profile shapes 
that will most efficiently operate at these angles must be 
developed. As mentioned before, because of requite- 
ments that dictate high turning angles and Mach Num- 
ber limits for stator and rotor exit and rotor inlet veloci- 
ties most blades are of the reaction rather than the im- 
pulse type. Therefore, the remarks to follow apply 
more specifically to the reaction type blade. 


The chord length of the blade must be sufficient to 
provide a Reynolds Number greater than 100,000 based 
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on fluid conditions at the blade passage exit. In addi- 
tion, high turning-angle blades require correspondingly 
high chord lengths in order to provide a sufficiently gen- 
erous radius of curvature for the channel walls. The 
so-called radius ratio of the curvature should be main- 
tained so that the local Mach Number over the convex 
surface of the blade does not exceed the Mach Number 
at the blade passage exit, in order to prevent local dif- 
fusion or adverse pressure gradient, with losses and 
possibly separation. 


The mean velocity in the channel should be constant 
or steadily increasing in order to avoid the possibility 
of diffusion in the passage. It is sometimes preferable 
to allow the air to diffuse ahead of the airfoil. This 
may be accomplished by a change in blade shape and 
use of a positive angle of incidence. 


The solidity or chord pitch ratio has an optimum 
value, which appears to be a function of turning angle 
and chord length and, in general, lies between 1.2 and 
17. Figs. 2 and 3 show the effect of solidity on ef- 
ficiency and deviation for some typical turbine blades. 
Deviation is the angle by which the airflow leaving the 
blade does not follow the physical exit angle of the 
blade. These curves were obtained from the Wright 
Aeronautical stationary cascade test rig shown in Figs. 
4and 5. With this rig it is possible to obtain sonic 


velocities in the nozzles. To check the information ob- 
tained from this rig, tests are run on a rotating rig such 
as shown in Fig. 6. The stand incorporates a combus- 
tion chamber for supplying hot gas and steam ejectors 
to simulate altitude exhaust. 

Trailing-edge thickness should not exceed 4 per ceyt 
of the blade pitch. Thick trailing edges result in high 
wake losses, which, in addition to reducing nozzle 
efficiency, may change the effective gas discharge 
angle. 

The leading-edge contour should, in general, be ellip- 
tical. The theory that a blunt nose, though less sensi- 
tive to incidence angle changes, gives a lower peak ef- 
ficiency has been fairly well established for isolated air- 
foils operated at low Mach Numbers. Tests of airfoils 
in cascade indicate that, for a given airfoil, there is an 
optimum nose shape, determinable by trial and error, 
which will provide both high peak efficiency and broad 
incidence angle range. 

Slightly negative incidence angles appear to give 
better efficiency than zero or slightly positive angles. 
Considerable care must be used in the evaluation of in- 
cidence angle effects, since, if incidence angle is changed 
with no change in blade shape, air turning angle is 
also changed. Thus, a change in efficiency, though 
credited to an incidence angle effect, may in reality 
result from the accompanying change in air turning 
angle. 

Stress and vibration factors influence the blade shape 
design especially for the rotor. For a given blade area 
at the tip section, the area at the root is sometimes dic- 
tated by the stresses due to centrifugal force. So also 
the center of gravity of the root section must be in 
nearly a radial line with the center of gravity of the tip 
section. 

The permissible twist of the blade, root to tip, is de- 
pendent upon the twist forces and the stress they pro- 
duce at various sections of the blade. In highly cam- 
bered blades, the choice of incidence angle is influenced 
considerably by the allowable twist. 


(IV) Miscellaneous Design Factors 


Radial clearances should be as low as possible, par- 
ticularly at the outer diameter where fluid density is 
highest and where area per unit of clearance width is 
greatest. This clearance is at present held below 0.020 
in. per in. of blade height. Axial clearance is not be- 
lieved to be important and is usually held within the 
range of !/, to '/2 in. 

The shape of the inner and outer contours of the an- 
nulus should be such that a smooth and gradual in- 
crease of annular area is provided in order to ensure 
against breakaway and excessive boundary-layer 
thickness. Desired contour shape is obtained by 
minor adjustments of the individual nozzle Mach 
Numbers. 
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The turbine machine is designed for a particular alti- the continuity equation W/A V/o. If the correct 


tude but must, in addition, perform satisfactorily at 
other altitudes. If turbine inlet conditions and exhaust 
pressure at any altitude are known, a valid performance 
calculation is possible. Unfortunately, the only 
method of obtaining such a calculation is the involved 
and laborious one of cut and try. Briefly, a weight 
flow is assumed, and the velocity in each successive 
nozzle is then determined by trial and error to satisfy 


weight flow has been assumed, the velocities in the in- 
dividual nozzles will have adjusted themselves to give 
the required discharge pressure at the last nozzle. Ii 
the required discharge pressure does not result from 
the calculation, a new weight flow involving new veloct- 
ties must be assumed and the entire calculation re 
peated. 
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UNSOLVED PROBLEMS 


A close study of the above design procedure indicates 
that there are many assumptions and ‘‘guesstimates”’ 
that still need to be accurately evaluated. One of the 
first of these problems is the analysis of the action of 
the boundary layer at the inside diameter (I.D.) and 
outside diameter (O.D.) of the blade annulus, as well 
as the action of the boundary layer on the blade itself. 
Its effect on the velocity vectors must be considerable 
and may mean that we do not obtain the velocity vector 
diagram that we calculate. When this effect is more 
completely evaluated, it may be possible to improve ef- 
ficiency by revising the vector diagram to take account 
of it. For example, if the boundary layer at the I.D. 
of the stator blades reduces the weight flow of air 
through the sections near it, it must tend to increase 
weight flow and velocities at other sections. The ques- 
tion of how far up the blade this can influence the vector 
diagram is, as yet, unanswered. 


The problem of obtaining high efficiency at high turn- 
ing angles is a paramount one. Calculations indicate 
that wall-friction and trailing-edge wake losses may, 
with careful design, be kept reasonably low, and, thus, 
if a blade shape can be drawn which prevents break- 
away and local shock, it would seem reasonable to 
expect to obtain high efficiencies with high turning 
angles, 


In the selection of blade chord length and shape we 
encounter the prgblem of Reynolds Number or turbu- 
lence. In the case of the isolated airfoil as in an air- 
craft wing, the turbulence is stated as a Reynolds 
Number that is a function of chord length and free- 
stream velocity ahead of the airfoil. In the turbine 
case, we encounter a series of airfoils in cascade having 
a much greater velocity difference, entering to leaving, 
and a much more curved airfoil. In the design of such 
blades we are interested in obtaining a turbulence that 
will give us a low drag coefficient that will not vary 
appreciably under a great number of operating con- 
ditions. 


The question of optimum incidence angle in a cascade 
of turbine blades is still unsettled. This angle is, of 
course, tied in closely with turning angle, blade shape, 
and passage area. For instance, a change of shape in 
more than the first 10 per cent of the chord of the 
blade in order to obtain a more negative angle of in- 
cidence may bring about a divergence in the passage 
area. 


We encounter the problem of selecting the optimum 
leaving velocity from the blade for a given relative inlet 
velocity or Mach Number. In other words, is it bet- 
ter to accelerate more and increase the favorable pres- 
sure gradient at the expense of the possible increase 
in drag due to the increase in friction at the higher 
velocities? A corollary to this is the question of how 
much of the loss in a blade is due to the turning 
or the work done and how much is due to inlet veloc- 
ity. 

When designing a blade and passage we are faced 
with the problem of deviation of the air from the actual 
blade angles especially at the exit. We need to know 
more about the flow pattern that dictates this devia- 
tion, but, still further, we must apply corrections such 
as overturning of the blade for this deviation in order 
to be sure to get the work done. However, overturning 
actually reduces the physical throat area, and, if throat 
Mach Numbers are already close to one, a reduction in 
weight flow will be involved unless the blade height is 
increased. This blade height change means a lower 
axial velocity after the blade and thus affects our veloc- 
ity diagram. 

More must be learned of blade vibration and damping 
resulting from aerodynamic forces, *because the vibra- 
tion stress adds to the centrifugal stresses and helps dic- 
tate the relation of blade thickness at the root com- 
pared with that at the tip section. This, in turn, pre- 
scribes a relation of the blade shape at the root to the 
blade shape at the tip. Solutions are desired to stress 
problems involving the allowable relative twist of the 
blade root to tip, since such twist helps govern the 
angle of incidence of the blade at the root for a given 
angle of incidence at the tip. 
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Appendix I 


HyPoTHETICAL 4,000-LB. THrRusT TURBOJET ENGINE. 
DETAILED CALCULATIONS FOR VECTOR DIAGRAM VE- 
LOCITIES AND NOZZLE AREAS 


Initial conditions, sea-level operation, 4,000-lb. 


thrust turbojet engine: 


Airplane speed 0 mph. (take-off condi- 


tions) 

Turbine inlet temperature = 1,500°F. 
Turbine inlet pressure 

(total) = 179.3 in. Hg 
Turbine weight flow = 70 lbs. per sec. 
Turbine r.p.m. = 7,800 
Assumed nozzle efficiency = 92 per cent 
Required turbine output = 11,720 hp. 
Maximum allowable turbine 

diameter = 32 in. 
Nomenclature 

V = velocity, ft. per sec. 


ti = temperature, °R. 
HH = actual enthalpy, B.t.u. per Ib. 
H’ == theoretical enthalpy, B.t.u. per Ib. 
AH = actual enthalpy change, B.t.u. per Ib. 
AH’ = theoretical enthalpy change, B.t.u. per lb. 
v = specific volume, cu.ft. per Ib. 
a = pressure, in. of mercury absolute 
PR’ = relative pressure, as used in Keenan and 
Kaye’s Thermodynamic Properties of Air 
W = weight flow, lbs. per sec. 
-A = area, sq.ft. 
a = area, sq.in. 
U = rotor speed, ft. per sec. 
AU = vector sum of tangential components of 


rotor entering and leaving velocities 


Determination of Root and Tip Design Radii 


Choose a rotor tip speed of 1,000 ft. per sec. Then, 


tip radius 


- 60 ) = 


1.222 ft. = 14.70 in. 


Assuming critical flow conditions at first stator noz- 
zle throat, determine required area from standard 


nozzle flow formula: 

W VTintet 

= ——— (3.82) = 
P. inlet 


(70) 


—————— (3.82) = 66 sq.in. 


‘179.3 
Assume an average nozzle angle of 18°. Then 
throat area 66 
‘ 1 = = = Zia. i . 
annular area sin 18° 0.3090 5 sq.in 


A porn 213.5 
root radius = _ 213.5 
12.21 in, 
12.21 
Uso = 147 (1,000) = 830 ft. per sec. 
From the formula 
UAU = HP(550g)/W 
determine the quantity UAU. Thus, 
11,720 X 550 X 32.17 
UAU, total root = 2,960,000 


70 


Assume a work distribution of 60 per cent first stage 
and 40 per cent second stage. Then, 


UAU first stage, root = 0.60 & 2,960,000 = 1,776,000 
° AU ip = 1,776,000/1,000 = 1,776 


“Assuming equal work root and tip, 
AU oot = 1,776,000/830 = 2,140 
Then, 
Uroot + AU soot = 830 + 2,140 = 2,970 


The total base length of the vector diagram, which is 
always equal to U + AU, is thus determined for the 
first stage. 

In order to provide for decreasing Mach Number 
through the machine, the diagram base is so divided 
that the right-hand side, which is actually the tangential 
component of stator discharge velocity, is slightly 
greater than the left-hand side. A right-hand side 
length equal to 54 per cent of the total base length is 
ordinarily used. Thus, 


Vi root tan = 0.54 X 2,970 = 1,601 ft. per sec. 
Try Vaxiai = 500 ft. per sec. Then, 
tan Groot = 500/1,601 = 0.312, 


@root = 17.3° 


Since 18° was chosen for an approximate average 
nozzle angle and since root angle is always smaller 
than average angle, the values of 500 ft. per sec. and 
17.3° Groot are satisfactory. Then, 


V1 coor = W1601? + 500? = 1,677 ft. per sec. 


Having determined the first stator root absolute veloc- 
ity and the axial and tangential components thereof, 
the remaining vector diagram calculations are set down 
in tabular form. 


Turbine Entrance Conditions 


Turbine entering enthalpy = H; = 398.47 root and 
tip; corresponding relative pressure, PR; = 3322 
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from air tables and the given inlet temperature of 


1,500°F.). 


Diagram Velocities and Enthalpy Drops 


(A) First Stage. 


1% 
Virel = 


Root—First Stator: 


1,677 ft. per sec. (as determined previously, 
M, = 0.835) 


= (1,677/223.7)? = 56.17 
= 56.17/0.92 = 61.10 B.t.u. per lb. 


398.47 — 56.17 = 342.30 (Ty = 1,755°R., 


PR, = 213.60) 
398.47 — 61.1 = 337.37 (PR,’ = 204.99) 
(204.99/332.2)179.3 = 110.8 in. Hg 


0.755(1,755/110.8) = 11.98 cu.ft. per Ib. 
V (1,601 — 830)? + 5002 = 919 


Root—First Rotor: 


Vexial = 550 ft. per sec. 
. Vereltan = U + AU — YV, abs tan = 2,970 — 

1,601 = 1,369 

V> rel = V1,369? + 550? = 1,475 ft. per sec. 

AH = (1,475/223.7)? = 43.5 

Hs = 342.3 + (919/223.7)? — 43.5 = 359.13 
— 43.5 = 315.63 

- = 1,656°R. (PR; = 169.98) 

AH’ = 43.5/0.92 = 47.3 

Hs; = 359.13 — 47.3 =, 311.83 (PR,’ = 
164.40) 

Ps = (164.4/213.6) 110.8 = 85.3 in. Hg 

is = 0.755 (1,656/85.3) = 14.68 cu.ft. per 
Ib. 

Veabs = V (1,369 — 830)? + 5502 = 770 


Tip—First Stator: 


Vitan = tanoor = 
tip 
Ge 1221) 
1,601 = 1,329 
14.71 
= V 1,329? + 5002 = 1,420 ft. per sec. 
AH = (1,420/223.7)? = 40.30 
A, = 398.47 — 40.3 = 358.17 (JT, = 1,813, 
PR, = 243.2 25 >) 
AH’ = 40.3/0.92. = 43.8 
Hy! = 398.47 — 43.8 = 354.67 (PR,’ = 
236.45) 
P, = (236.45/332.2) 179.3 = 127.7 in. Hg 
= 0.755 (1,813/127.7) = 10.74 cu.ft 
per Ib. 
Virel = = V (1,329 — 1,000)? + 5002 = 599 
Tip—First Rotor: 
Vaxiat tip Vaxiat = 550 ft. per sec. 
Vareltan = U +- AU — V, abs tan = 1,000 + 
1,776 — 1, nad = 1,447 
V2 rel = V1,447 4472 + 5502 = 1,548 ft. per sec. 


AH = (1,548/223.7)? = 47.9 

Hs = 358.17 + (599/223.7)? -— 47.9 = 
365.34 — 47.9 = 317. 44 

Ts = 1,663°R. (PR; = 172.72) 

AH’ = 47.9/0.92 = 52.0 

H;! = 365.34 — 52.0 = 313.34 (PR = 
166.62) 

Ps = (166.62/243.25) 127.7 = 87.4 

Us = 0.755 (1,663/87.4) = 14.38 cu.ft. per lb. 

Veabs = V(1,447 — 1,000)? + 550? = 709 


(B) Second Stage. Root—Second Stator: 


UAU = 2,960,000 — 1,776,000 = 1,184,000 

AU = 1,184,000/830 = 1,428 

U+ AU = 830 + 1,428 = 2,258 

right-hand side = 0.54 X 2,258 = 1,219 

V axial = 600 ft. per sec. 

Vi = V1,2192 + 6002 = 1,359 

AH = (1,359/223.7)? = 36.9 

AH’ = 36.9/0.92 = 40.1 

Hs = 315.63 + (777/223.7)? — 369 = 
327.68 — 36.9 = 290.78 

Ts = 1,564° R. (PRs aad 135.7 5) 

H,’ = 327.68 = 40.1 = 287. 58 (PR,’ = 
131.70) 

Pe = (131.70/169. 66.0 

% = 0.755(1564/66.0) = 17.9 

Virel = V(1,219 — 830)? 6002 = 715 


Root—Second Rotor: 


4 = 620 ft. per sec. 


axial 

V2 rel tan = 2,258 — 1,219 = 1,039 

Verel = V1,0392 + 6202 = 1,210 

AH = 1,207/223.7 = 29.1 

H; = 290.78 + (715/223.7)? — 29.1 = 
301.01 — 29.1 = 271.91 

ti = 1,493°R. (PR; = 113.40) 

AH’ = 29.1/0.92 = 31.7 

H,’ = 301.01 — 31.7 = 269.31 (PR; = 
110.50) 

P; = (110.50/135.75)66.0 = 53.8 

V7 = 0.755(1,493/53.8) = 21.00 cu.ft. per Ib. 

V2abs. = V (1,039 — 830)? + 6202 = 654 


Tip—Second Stator: 


V; tan = (12.21/14.71) 1,219 = 1,010 

Vi = V1,010? + 600? = 1,175 

AH = (1,175/223.7)? = 27.55 

He = 317.44 + (709/223.7)? — 27.55 = 
327.51 — 27.55 = 299.96 

Ts = 1,598°R. (PRs = 147.75) 

AIT’ = 27.55/0.92 = 29.9 

HI,’ = 327.51 — 29.9 = 297.61 (PR,’ = 
144.60) 

Ps = (144.60/172.72)87.4 = 73.1 
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rel 


IL 


0.755 (1,598/73.1) = 16.5 
= V(1,010 — 1,000)? + 6002 = 600.1 


Tip—sSecond Rotor: 


Vasial = 620 


%reltan = U + AU — JV, abs tan = 1,000 + 
(1,184,000/ 1,000) — 1010 = 1,174 

Verel W1,174? + 6202 = 1,327 

AH = (1,327/223.7)? = 35.2 

HH; = 299.96 + (600.1/223.7)? — 35.2 = 
307.15 — 35.2 = 271.95 

r = 1,493°R. (PR; = 113.40) 

AH’ = 35.2/0.92 = 38.3 

HH,’ = 307.15 — 38.3 = 268.85 (PR,’ = 
110.01) 

P; = (110.01/147.75) 73.1 = 54.5 . 

07 = 0.755(1,493/54.5) = 20.7 cu.ft. per Ib. 

V2abs = V (1,174 — 1,000)? + 620? = 644 

Overall Efficiency. Root: 

actual overall pressure ratio = 179.3/53.8 = 3.33 

PR,’ overall = PR;/3.33 = 332.2/3.33 = 99.8 


H’; overall = 259.13 
AH’ overall = H; — Hy,’ overall = 398.47 — 259.13 


= 139.34 
UAU = 2,960,000 
QUAU/2gJ 
AH’ — (V,y/223.7)? 
118.30 
139.34 — (654/223.7)2 
118.3/130.79 = 0.91 
Tip: 


actual overall pressure ratio = 179.3/54.9 = 3.261 
PR,’ overall = PR3/3.261 = 332.2/3.261 = 101.8 
H,’ overall = 261.0 

AH’ overall = H; — H,’ overall = 398.47 — 261.0 


= 137.47 

UAU = 2,960,000 
E 2UAU/2gJ 


118.3 
137.47 — (644/223.7)? 
118.3/129.8 = 0.915 
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CALCULATION OF ANNULAR AREA 


To calculate the annular area required to pass the 
design weight flow, the basic continuity equation jg 
used. Thisis: Q = AV, Wr =Al,or WA = Vy 

In averaging the IV/A for the entire nozzle, the values 
for the two sections considered are weighted in direct 
proportion to the magnitude of the radii at which the 
sections are located. This is done because of a fixed 
annular height; the area of an annulus increases dj. 
rectly with the radius. Thus, 


(Vz/ve)Re + (Vr/v)R, 
A average Re + R, 


W(Re + R:) 
(Vr/vr)Re + 


For first stator: 


70(12.21 + 14.7) X 144 
(500/11.98)12.21 + (500/10.74)14.7 
227.0 sq.in. 


@ annular = 


For first rotor: 


70(12.21 + 14.7) x 144 7 
(550/14.68)12.21 + (550/14.38)14.7 
265.9 sq. in. 


a@ annular = 


For second stator: 


70(12.21 + 14.7) 144 
(600/17.9)12.21 + (600/16.5)14.7 
287.0 sq.in. 


@ annular = 


For second rotor: 
70(12.21 + 14.7) & 144 


(620/21.0)12.21 + (620/20.7)14.7 
338.0 sq.in. 


@ annular = 


Having established the required annular areas, the 
blade heights required are easily determined. Blade 
entrance and exit angles are obtained from the vector 
diagram, using proper correction factors for deviation 
of air flow from the physical angles. The amount of 
correction depends upon the characteristics of the par- 
ticular blade shape involved. 
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The Aerodynamic Forces on an Oscillating 
Flap at Supersonic Speeds 


JOHN W. MILES* 
University of California at Los Angeles 


SUMMARY 


Using the general results developed in an earlier paper, the 
complex flutter coefficients A,, (hinge moment due to plunging), 
Arq (hinge moment due to pitching), A.g (lift due to flapping), 
A,g (pitching moment coefficient due to flapping), and A.g 
(hinge moment coefficient due to flapping) are found. All of 
these coefficients are expressed in terms of the previously known 
coefficients Aga, Aak, Aca, and Ay. Numerical results are given 
in the form of curves for M = 1.2, 1.4, 1.6, and 2.0, the real and 
imaginary parts of the coefficients being plotted vs. & for a range 
of k = Oto 1.4. In the case of Agg, Asg, and A.g, the results 
are applicable to any flap, while the numerical results for Aq and 
Apa are given for a 50 per cent chord flap. 


NOTATION 
a, = amplitude of sth degree of freedom 
b = semichord of airfoil 
C = fraction of semichord of the flap hinge aft of the mid- 
chord 
h = plunging motion 
j = (—1)!/2 = imaginary unit 


k = (wb/U) = reduced frequency 


p = pressure (gage) 

t = time 

x,2 = airfoil coordinates 

A,;,s = dimensionless flutter coefficient giving force in rth 
degree of freedom due to motion in sth degree of 
freedom 

G(v) = Green’s function; cf. Eq. (10) 

Jo = Bessel function, first kind, zero order 

L = iift 


M == Mach Number; also mid-chord stalling moment 

U free-stream velocity 

a = angle of attack 

8 = (M? — 1)!/2 in. Egs. (4), (5), (6), (7), (10); 
denotes flap deflection 

angular frequency of harmonic motion 


otherwise 


INTRODUCTION 


I AN EARLIER PAPER! the Green’s function giving the 

pressure on an airfoil in a supersonic stream due to 
an arbitrary harmonic velocity distribution was de- 
rived. The result was utilized to find the lift and mid- 
chord pitching moment due to plunging and pitching 
motions and, later,” to find the response of an airfoil 
entering a sharp edged gust or to a sudden change in 
angle of attack or flap deflection. The present paper 
will extend the results to include the forces on an oscil- 
lating flap, because of both its own motion and the mo- 
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tions of the airfoil as a whole. The notation and as- 
sumptions, namely, small perturbations, the specifi- 
cation of the velocity at the mean location of the airfoil, 
and two-dimensional flow are the same as in reference 1. 


VELOCITY DISTRIBUTIONS 


The motions of the airfoil due to plunging (4) and 
pitching (a) about the mid-chord are defined by Eqs. 
(17) and (18) of reference 1; to these a third motion, 
elevon flapping (8), will be added. The displacements 
in these three modes of motion are described by 


t) = (1) 

t) = —a,bxe™ (2) 

za(x, t) = —agb(x — cle, x (3) 
=, 


The positive displacements in these modes are illus- 
trated in Fig. 1. 


FLUTTER COEFFICIENTS 


The lift, mid-chord stalling moment, and hinge 
moment (positive when acting to increase 8) due to the 
motion 2,(¢) are given by 


1 
L(t) = (= Jae. = fou thdx (4) 


1 
M(t) = (= ase Ay =2 xp.(x, (5) 


1 
N,(t) = wi Jase Ay = (x — c)p,(x,t)dx (6) 


and are nondimensionally described by the flutter coef- 
ficients A,,. 

In deriving the flutter derivatives it is convenient to 
express the pressure distribution in the form 


X=C 
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where /, is a nondimensional distribution in the vari- 
able k(x +.1). From Eqs. (12), (20), and (21) of 
reference 1 it is found that 


Pr(ku) = G(ku) + So™ G(v)dv (8) 
ku 


Palku) = ( 1+ prlku) — pr(v)dv (9) 
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G(v) = (M/8?)0] (10) 


where (x + 1) has been replaced by wu. 
Introducing Eq. (7) in Eqs. (4)-(6), the flutter 
coefficients may be rewritten 


ps(ku)du (11) 
mk Jo 
a3 

= (u — 1)p,(ku)du (12) 
rk Jo 


A,,(k, M) = 


Ags(k, M) 
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k’ = [(1 + 0)/2]k (16) | r 4 
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l+c 
Ac (A k, u)| (17) Now from Eqs. (9), (11), and (12), it is found that 


I To evaluate A'ya, Eq. (9) may be substituted in i) 
Eq. (15) to obtain Agalk, M) + Acalk, M) = (: + b [Ag(k, M) + 


2) 1 2 ku 7 
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'w (19) 
(w Pr(v)dv (18) Comparing Eqs. (17), (18), and (19) yields the result: 


M, = (755) arate, at, 


0247 E k, u) + (2+ u) (20) « 
2 
0 Substituting Eqs. (18) and (20) in Eq. (14) yields 
the desired coefficients: 
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ForRCE AND MOMENT DUE TO OSCILLATING FLAP 


In computing the forces due to the motion of an os- 
cillating flap, as described by Eq. (3), it may be con- 
cluded immediately (since the flow is supersonic) that 
the pressure distribution due to the motion van- 
ishes except over the flap itself. It follows that 
the forces due to an oscillating flap are analogous to 
those due to an oscillating (in pitch) airfoil of chord 
[(1 — c)/2]26, differing only in the fact that the rotation 
of the flap and the desired hinge moment are about the 
leading edge, whereas these quantities were computed 
with respect to the mid-chord in the case of the pitching 
airfoil. It therefore follows directly that 

2 


3 
Avg (k, M, c) k, u) 
1-¢ 2 


1 


2 4 l-—c 
(+ -) M,c) = Aag k, u) 
2 $ 
Auf k, Mt) (=) A.g(k, M,c) (24) 


where the factor of [(1 — c)/2]* enters in Eq. (23) 
because b? enters the nondimensionalizing coefficient in 
Eq. (4), where [(1 — c)/2]* enters in Eq. (24) because 
of 6‘ in Eq. (5), and where [(1 — c)/2]k is used because 
k is defined in terms of b. 
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The remaining coefficient A,g, giving the pitching 
moment due to flapping, is given by 


Aas(k, M, c) = Anglk, M, c) — cAcs(k, M,c) (25) 


as may be seen from Eqs. (11)—(13). 


NUMERICAL RESULTS 


The coefficients A, are conveniently broken into 
their real and imaginary parts, viz: 


An = Rr + jln (26) 


Using the results of reference 1, the real and imag- 
inary parts of the above flutter coefficients are plotted 
in Figs. 2-9. The coefficients Rys, Rp, and are 
plotted so that the results are valid for all values of 
c, while those for the remaining coefficients are plotted 
only for c = 0, corresponding to a 50 per cent chord 
flap. The.results for other values of ¢ are easily ob- 
tained from the results of references 1 or 3 and from 
Eqs. (21), (22), and (25) above, the last being par- 
ticularly simple. 
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